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The  variational  effect  of  wind  shears,  convection,  and  temperature 
gradients  on  atmospheric  density  distributions  characterize  atmospheric 
turbulence.  The  temperature  structure  parameter,  CT2,  and  the  index  of 

refraction  structure  parameter,  Cn2,  quantify  two  density  distributions 

that  arise  in  acoustic  and  electromagnetic  propagation.  [Ref.  1] 

The  density  distributions  quantified  by  Cn2  and  C,2  induce  phase 

fluctuations  in  propagating  electromagnetic  waves.  The  density 
distributions  of  atmospheric  turbulence  degrade  the  spatial  and  temporal 
coherences  of  initially  coherent  electromagnetic  radiation.  Similar,  more 
pronounced  effects  occur  in  the  propagation  of  acoustic  waves.  [Ref.  1] 

A  number  of  techniques  for  correcting  these  problems  have  been 
proposed,  particularly  with  respect  to  electromagnetic  propagation.  These 
corrective  techniques  include  the  use  of  modular  mirrors  [Ref.  2],  optical 
phase  conjugation  [Ref's.  2  and  3],  optical  and  digital  signal  processing 
[Ref.  A],  and  simple  avoidance  of  turbulent  air  masses.  The  appropriate 
timing  of  transmissions  or  suitable  location  of  transmission  and  reception 


facilities  accomplishes  avoidance  [Ref’s.  5  and  6],  These  techniques  can 
be  applied  alone  or  in  combinations. 

A  knowledge  of  the  history  of  the  density  or  refractive  structure 
along  a  proposed  transmission  path  is  required  to  effectively  employ  a 
number  of  these  compensatory  techniques.  Time-averaged  CT2  or  Cn2 

profiles  charactize  this  history  [Ref's.  1  and  7], 

Wroblewski  and  Weingartner  developed  and  Moxcey  refined  a  high 
resolution,  computer  controlled  acoustic  sounder,  or  echosounder,  system. 
Their  computer  controlled  echosounder  system  uses  a  transmitted  pulse  of 
acoustic  energy  as  a  probe  of  atmospheric  structure.  The  interaction  of 
the  acoustic  wave  packet  with  atmospheric  turbulence  is  used  to  develop  a 
time-averaged  C,2  profile  for  a  short  spa  dal  range  [Refs.  5,  6  and  81.  The 

acoustic  sounder  may  also  be  used  to  verify  certain  long  range  optical 
measurements  of  turbulence  [Ref.  6], 

Development  of  the  CT2  profile  uses  a  CT2  expression  developed  from 

the  echosonde  equation  summarized  by  Neff  [Ref.  7]  and  the  empirical 
acoustic  backscatter  cross  section  per  unit  volume  expression  of  Tatarski 
[Ref.  9]  iRef's.  5  and  6]  The  accuracy  of  the  absolute  values  of  this  CT2 

profile  depends  on  the  calibration  of  the  acoustic  antenna  parameters 


2 


contained  in  the  CT2  expression.  The  antenna  parameters  of  interest  are 
Et,  the  efficiency  of  conversion  of  electrical  to  acoustic  power,  and  Er,  the 

efficiency  of  conversion  of  acoustic  power  to  electrical  power. 

In  this  thesis  acoustic  energy  Lackscattered  off  acoustically  hard 
spheres  was  used  to  determine  Er  and  Et  Acoustic  range  equations  were 

software  implemented  to  perform  the  necessary  calculations.  Data  and 
calculational  results  provided  first  order  verification  of  the  calibration 
process. 
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A.  CT2  REVIEW 

Wroblewski  and  Weingartner  applied  the  echosonde  equation, 

Pr=  Er[P(Et][e— ][^^,0][Ae]  (1) 

summarized  by  Neff  [Ref.  7]  and  the  acoustic  backscattering  cross 
section  per  unit  volume  expression  given  by  Tatarski  [Ref.  9]  as 

V  C  2 

=  0.0039  k  3  -1_  (2) 

to  develop  the  volume  averaged  expression 

ct2=  J_  Iii-LJ-SLR*,**  (3, 

00039  ErE,  k's  CT  AG  P, 

where 

•  Pr  is  the  electrical  power  of  the  reflected  signal  returned  to  the 
acoustical  antenna  array, 

•  Er  is  the  efficiency  of  conversion  of  the  returned  acoustic  power  to 

Pr 

•  Pt  is  the  electrical  power  supplied  to  the  acoustical  array, 

«  Et  is  the  efficiency  of  conversion  of  Pt  to  acoustic  power, 

•  a  is  the  average  attenuation  per  unit  distance  along  the  transmission 
path, 

•  R  is  the  range  from  the  array  to  the  target  or  backscattering  volume. 


•  c  is  the  average  speed  of  sound  along  the  transmission  path, 

•  t  is  the  length  in  time  of  the  transmitted  acoustic  pulse, 

•  A  is  the  antenna's  aperture  area, 

•  G  is  the  antenna's  effective  aperture  factor, 

•  k  is  the  wavenumber  of  the  incident  acoustic  energy,  and 

•  T0  is  the  average  temperature  (in  degrees  Kelvin)  along  the 

transmission  path. 

[Ref's.  5  and  6] 

The  conditions  at  the  time  and  location  of  each  application  determine 
the  value  of  the  terms  T0,  k,  c,  x,  Pr,  Pt,  R  and  a  in  the  CT2  expression. 

The  terms  A  and  G  depend  on  the  antenna  design. 

Adverse  field  effects  such  as  dust,  debris  and  vibration  loosened 
connections  cause  the  remaining  two  terms,  Er  and  Et,  to  change  over  the 

life  of  the  system.  This  chapter  develops  the  theoretical  expressions 
supporting  a  determination  of  Er  and  Et. 


B.  ACOUSTIC  RANGE  EQUATIONS  AND  EFFICIENCIES 

In  the  following  two  sections  the  acoustic  range  equation  of  Neff 
[Ref.  7]  is  re-developed  for  targets  of  small  area  relative  to  the 
cross-section  of  the  ensonifying  field.  Neff  [Ref.  7],  Skolnik  [Ref.  10],  and 
Probert-Jones  [Ref.  1 1  ]  were  used  for  guidance. 
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As  defined  previously,  the  electrical  power  supplied  to  the 
acoustic  array  is  Pt.  it  is  converted  at  an  efficiency  Et  to  an  emitted 

acoustic  power, 

Pa  =  Pt  Et .  (4) 

Pfl  is  spread  over  a  4tt  solid  angle. 

A  far  field  range  R  is  defined  by  the  condition 

A/2  »  ILI  -  IRI  (5) 

for  the  geometry  in  Figure  1.  A  is  the  wavelength  of  the  emitted  acoustic 
wave. 


Array 


Figure  1:  Calibration  (Far  Field)  Geometry 
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The  acoustic  intensity, 


l8  -  P8  /  4tTR2,  (6) 

describes  the  emitted  acoustic  power  density  of  an  isotropic  emitter  at  a 
far  field  range  R.  Ifl  is  in  units  of  [W/m2]. 

If  the  emitted  power  is  directional,  the  intensity  la  is  modulated 
by  a  geometrical  gain, 

G(Q)  =  GolF(O)|2  =  GolF(0,e)l2  (7) 

Note  this  is  not  the  same  G  found  in  equation  ( 1 ).  Now, 

lF(0,0)j 2  =  1,  (8) 

and  Go  is  thefore  the  centerline  gain  of  the  acoustic  array.  [Ref.  11]  The 

centerline  direction  (0,0)  is  taken  to  be  normal  to  the  face  of  the  acoustic 
array.  The  acoustic  intensity  lfl  is  now  written 

I,  =  J±__  G (Q)  (9) 

In  a  real  propagating  medium  intensity  attenuation  occurs  along 
the  transmission  path.  The  radial  term  e"0^  accounts  for  this 
attenuation  when  the  atmospheric  path  is  assumed  to  be  reasonably 
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homogeneous.  The  acoustic  intensity  is  then 


4iTR 


G(O)  e 


-aR 


00) 


A  target  of  cross  sectional  area  At  t  intercepts  a  power 


PtgL  =  *S  Atgt 


(11) 


Expanding  this  expression  with  equation  (10)  yields 

P  =[5l 
W  L  4tT 


(12) 


For  a  target  of  cross  sectional  area  Atgt  with  a  sufficiently 
small  solid  angle  Algl  /R2  little  change  in  intensity  occurs  about  the 

direction  0.  Equation  ( 1 2)  is  now  arranged  as  a  one  way  range 
equation  along  the  centerline  direction  (0,0), 


A*„4  t  2 


(13) 


Isolating  EtG0  in  the  above  expression  yields 


E.G,  =  e“R 


(14) 


V  * 
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Electrical  power  is  expressed  as 


P  = 


v  2 

vrms 

Z 


V  ms 

Z  ' 


(15) 


where 


•  Vpms  is  the  empirical  root  mean  square  voltage, 

•  V  is  V  2,  and 

ms  rms  * 

•  1  is  the  electrical  impedance. 


EtGQ  is  now  written 


EtG0 


=  4ff  R 


Z  e  CtR  tqt 

V 


(16) 


The  ratio  Ptgt/Atgl  is  defined  as  the  acoustic  intensity  lr 


received  by  the  target.  Therefore, 


E,60  =  4«R2 


eKR  |. 


(17) 


2.  Two  Way  Range  Equation  and  Er 

For  an  acoustically  hard  target  the  total  power  reflected  is  Ptgt. 

The  power  reflected  from  the  target  surface  has  a  normalized, 
directional  intensity  distribution  described  by  the  differential  scattering 
cross  section,  <r 
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The  direction  from  the  target  to  the  array  is  the  backscattered 
direction.  The  differential  scattering  cross  section  in  the  backscattered 
direction  is  calied  the  normalized  backscattered  cross  section,  <rb.  <rb  is 


defined  as  the  power  reflected  toward  the  source  per  unit  solid  angle, 
normalized  by  the  incident  power  density  over  a  4it  solid  angle. 


The  backscattered  intensity  from  the  target  is 


(10) 


An  array  of  aperture  area  A  at  range  R  intercepts  an  attenuated,  received 
acoustic  power, 

Pr.  ’  V  A  -  C9) 

Substituting  equation  ( 1 8)  into  equation  ( 1 9)  yields 

P  =  „  e"“R  A  (20) 

r”  4ffR2  b 


Pra  is  converted  to  returned  electrical  power,  P  at  a  return 


efficiency,  Er; 


Since  electrical  power  is  expressed  as 


Vrms  _  V  ms 
Z  Z  ' 


(15) 


electrical  power  representing  the  returned  acoustic  signal  from  a 

target  placed  along  an  emitter's  direction  (0,0)  is  written  as 

(Vms)r  .  (Vms)t[EtGo  e-«R  AtqfirErgh  e~«R  A  ~j  (22) 
Zr  "  Zt  t  471  r2  JL  4TJ  r2J  . 


For  a  passive,  linear  acoustic  array  circuit  the  concept  of 
impedance  reciprocity  is  invoked  [Ref.  12]  to  yield  a  two  way  range 
equation  for  scattering  from  a  small  target; 

i 


p  -  f  (  Vms)t 
"L(Vms)r 


EtEr 


e"20CR 
2  c 

(4*r) 


(23) 


This  is  written  to  isolate  Er; 

E  .  (VmsV  (4TIR2)2  e2«*  J_  ,, 
r  '  (Vms),  E,G0  A  Aw  «b  - 

A  part  of  the  return  signal's  power  is  unwanted  noise. 
Subtracting  noise  power  from  the  returned  acoustic  power  yields 


2^ 


r  =  (V  ms  )r~(V  ms)n  (4TfR^) 
(Vms)t  Et60 


20CR 


A  A 


tgi 


(25) 


As  return  voltages  are  quite  small  they  are  analyzed  after 


electrical  amplification.  Calling  the  electrical  amplification  a  gain 


allows  introduction  of  the  final  modification  to  E„: 

6  r 

F  _  (Vms)r-(Vms)n  (4WR2)2  e2KR  I 
r  ’  (Vms>,  Ge2  E,G0  A  Atgt  a„ 


C.  ATTENUATION  AND  RANGE  DETERMINATION 


(26) 


Employing  equations  ( 1 7)  and  (26)  in  the  calibration  process  requires 
calculating  the  attenuation  a  and  the  range  R  from  the  array  to  the 
target. 

1.  Attenuation 


The  calculated  attenuation  cc  for  an  atmospheric  propagation 
path  is  assumed  to  be  the  combined  result  of  molecular  and  classical 
absorption  in  the  atmosphere.  The  classical  absorption  otc,  is  attributed 

to  heat  conduction  and  viscous  effects  while  the  larger  molecular 
absorption  amo]  is  attributed  to  the  excitation  of  internal  energy  modes  of 

atmospheric  gases  by  the  propagating  sound  energy.  [Ref.  7] 

Information  from  Neff  [Ref.  7],  Businger  [Ref.  13],  Neiburger 
[Ref.  14],  and  Fuller  [Ref.  15]  is  applied  to  make  a  reasonable  determination 
of  a.  The  determination  is  made  using  the  average  measured  temperature 
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Tc  in  degrees  Celsius  along  the  propagation  path,  the  measured 


atmospheric  pressure  P  in  millibars,  and  the  measured  relative  humidity 
Rh  in  percent. 

The  molecular  attenuation  coefficient,  a  .,  is  expressed  by  the 


empirical  relationship 


’ratio 


amoi is  in  units  of  [_dB/m], 

Now, 

<xmax  =  0.0078  fm  (T*r2-5  e7  77(,-,/T*)  (28) 

is  the  maximum  absorption  at 

(10  +  6600  (100f-)  +  44,400  (100f-)2)  P*  (2g) 

f  is  in  units  of  [Hz], 

m 


Also, 


L  =  Es Rh  (30) 

P  P  "  E3(  1  ♦  Rh) 


is  the  mole  ratio  of  water  vapor  in  the  atmosphere.  The  water  vapor 
pressure  e  is  expressed  in  millibars. 
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(31) 


Es  =  io^-2353/To> 
expresses  the  saturation  vapor  pressure  Es.  TQ  is  the  average  absolute 

temperature  along  the  transmission  path,  as  defined  in  Section  A. 

T*  =  ( 1 .8  Tc  ♦  492)  /  5 1 9  ,  (32) 

and 


P*  =  P/  1014. 


(33) 


Finally, 


^ratio  "  ^  » 


m 


(34) 


where  f  is  the  operating  frequency  of  the  echosounder  in  [Hz],  [Refs.  7, 
13,  1 4  and  1 5] 

The  classical  attenuation  coefficient,  a  .,  is  approximated  by 


<xcl  =  1.74*  10_,0(f)2 


(35) 


ac|  is  in  units  of  l-dB/ml  [Ref.  7] 


For  a  total  attenuation 


«t  =  aci +  amol 


(36) 


and  the  relationship 
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-  «xtR  =  10  log  e-ctR  (37) 

the  average  attenuation  per  meter  is  finally  written  as 

«  =  )n)0  <38) 

Equation  (37)  is  derived  from  unit  definitions  and  equation  one  of  Neff 
[Ref.  71 

2.  Range 

The  coordinate  system  for  the  backscattered  cross  section 
determination  uses  the  center  of  the  target  spheres  as  the  coordinate 
system  origin.  Therefore,  the  radius  a  of  a  target  sphere  and  the  time  of 
return  tr  of  an  acoustic  signal  from  the  acoustic  echosounder  are  used  to 


determine  the  range  R  from  the  target  to  the  echosounder; 


(39) 


c  is  the  average  speed  of  sound  aung  the  transmission  path,  and 

V  “  ^receive  "  ^transmit  ■ 

For  the  range  of  temperatures  expected  in  normal  echosounder  operation 
the  speed  of  sound  in  meter  per  seconds  in  dry  air  is 

-  20.05  VT7  .  (41) 
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For  moist  air  c  becomes 


'  ■  Sri*  -  Sry<*  *0.I4»/P).  (42) 

where  el P  is  defined  in  equation  (30).  [Ref.  1] 


D.  BACKSCATTERED  CR05S  SECTION  DETERMINATION 

Employing  equations  (17)  and  (26)  in  the  calibration  process  also 
requires  calculating  <rb,  the  backscattered  cross  section.  Recall  from 

Paragraph  B.2  that  the  backscattered  cross  section  is  the  differential 
scattering  cross  section  in  the  direction  of  the  energy  source. 

1.  Initial  Conditions  and  Assumptions 

In  determining  the  differential  scattering  cross  section  <r  it  is 
assumed  the  incident  wave  state  is  unchanged  after  scattering.  That  is, 
for  a  stationary,  acoustically  hard  target  it  is  assumed 

^  ~  ^incident  ~  ^scattered  ’ 

and  that  incident  k  has  good  definition,  i.e.  the  incident  wave  can  be 
considered  monochromatic  [Ref.  16], 

Selection  of  acoustically  hard  spheres  as  targets  further 
simplifies  the  determination  of  the  differential  scattering  cross  section 
solution.  This  target  selection  sets  the  boundary  condition  at  the  target 
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surface  as  the  Neumann  boundary  condition 


(44) 


V  is  the  velocity  potential  defined  by  the  acoustic  propagation  velocity  c 
as 


c  *  7  V  .  (45) 

n  lies  along  the  direction  normal  to  the  target  surface.  Also, 

P  =  -Po-ff  (46) 

for  the  pressure  p  and  rest  density  po  of  the  propagating  medium.  [Ref.  16]. 


Packet  of  N  Waves ;  Coordi nate 

W,NA»a.  System. 

Figure  2.  Incident  Wave  and  Target  Description 


The  target  sphere  is  aligned  along  the  centroid  of  the  main 
acoustic  lobe  of  the  acoustic  arrray  echosounder,  that  is,  along  the 
direction  (0,0),  normal  to  the  face  of  the  array.  The  target  sphere  is 
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selected  to  be  sufficiently  small  with  respect  to  incident  beam 
dimensions  so  that  an  essentially  uniform  intensity  is  intercepted  by 
the  sphere  [Ref.  1 7j.  Since  the  target  intercepts  a  very  small  proportion 
of  the  centroid  of  the  transmitted  sounding  lobe  no  accommodation  for 
sounding  lobe  geometry  over  the  target's  angular  area  is  made. 

Acoustic  gain  is  simply  taken  as  G0. 

Target  placement  is  also  at  a  far  field  range.  This  allows 
spherical  incident  waves  to  be  treated  as  planar  incident  waves  and 
places  the  target  outside  the  range  associated  with  the  ring  time  of  the 
echosounder  [Ref.  8],  Placement  of  the  target  is  also  close  enough  to 
the  array  to  support  an  assumption  of  transmission  medium  homogeneity 
along  the  transmission  path. 

Finally,  the  potential  describing  the  target  falls  off  faster  than 
1/r.  The  describing  potential  for  an  acoustically  hard  sphere  with 
radius  a  is 

r  r  i  a 

V(r) =  (47) 

-  0  r  i  a 

Thus,  for  r  >  a  V(r)  falls  off  faster  than  1/r. 
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Note  that  with  the  exception  of  the  specific  surface  boundary 
condition  for  acoustic  waves  the  above  conditions  and  assumptions  closely 
correspond  to  those  for  plane  electromagnetic  waves  incident  on  a 
perfectly  conducting  sphere,  and  for  perfectly  elastic  nuclear  scattering 
in  a  center  of  mass  coordinate  system.  [Refs.  16,  18,  19  and  20] 

2.  The  Incident  Wave 

The  asymptotic  form  of  an  incident  plane  wave  approaching  a 
scattering  target  but  beyond  the  influence  of  the  target's  potential  V(r)  is 
the  undisturbed,  or  free,  plane  wave.  In  the  spherical  polar  coordinates 
(r,0,e)  illustrated  in  Figure  3  a  free,  monochromatic  plane  wave  of 
unit  amplitude  propagating  in  the  direction  Q0  given  by  (0o.eo)  is 

f(r,t)  =  e  »kr[cos(0o)cos(  0)  +  sin(0o)cos(  e  -  e0)]  e-wot  (4Q) 

[Ref.  16] 


Figure  3:  Spherical  Polar  Coordinates 
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For  a  wave  of  unit  amplitude  incident  along  the  z  axis  ,  Oo  »  0 


and  the  free  wave  f(r,t)  is  presented  as 

f(r,t)  =  e  ikr  e"10)t  =  e  ikz  e_kDt .  (49) 

Note  this  is  a  solution  with  axial  symmetry  for  the  free  wave  equation 
[Ref.  181 

The  wave  equation  describing  the  behavior  of  the  undisturbed  or 
free  monochromatic  wave  J(r\t)  is 

( V2--L-^,)f(r,t)  =  0.  (50) 

at 


c  is  the  speed  of  propagation  of  the  free  wave.  If  the  wave 
function  f(r,t)  has  a  simple  harmonic  time  dependence,  i.e., 

j(r.t)  *  e'“>*  f(r) ,  (51) 

then  for  a  wave  number 

k  h  oj/c  =  2t 1/A  (52) 


the  free  wave  equation  becomes  the  Helmholtz  equation 

(v2  +  k2]f(r,t)  =  0 


[Ref.  211 


(53) 
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The  stationary  solution  f(r)  of  the  free  wave  equation  can  be 
written  as  the  product 

f(r)  =  1  ,(Qi >  J2(q2) f3(q3)  (54) 

for  orthogonal  curvilinear  coordinates  q,,  q2,  q3  [Ref.  161.  Recall  an 

undisturbed  plane  wave  has  a  cylindrical  symmetry  about  the  axis  of 
propagation  [Ref.  22].  Using  the  spherical  polar  coordinates 

d]  =  r,  q2=0,  q3=  e,  (55) 


the  stationary  solution  describing  an  undisturbed  plane  wave  propagating 
along  the  z  axis  is 


f(D«  ft<r) j2(0)  . 


(56) 


Substituting  f(r)  into  the  free  wave  equation  expressed  in 
spherical  coordinates  and  separating  the  coordinate  dependencies  yields 
the  separated  free  wave  equation, 

(kr)2’f  &r2&<t  =  ‘  f— 1 —  f-  <57> 

^  ur  ar  •  j2smb  der  60  <- 

Using  a  separation  constant  of  1(1+ 1 )  the  equation  for  the  angular 


dependency  of  the  free  wave  becomes 


7— 1 - -4-  <sin0-^-f_)  ♦  l(l  ♦  1)  =  0  (58) 

)2  sin  0  d#  dff  £ 
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Letting 


cos(0)  =  x  (59) 

in  the  angular  equation  yields 

[(1  -x2)t4  -  2x  +  1(1  +  1)]  f2  =  0  .  (60) 

1  dx2  dx  c 

The  previous  equation  is  known  as  the  Legendre  differential 
equation  [Ref.  21).  The  most  useful  solutions  to  the  Legendre  differential 
equation  are  the  Legendre  polynomials,  designated  by  P,(x).  They  can  be 

generated  by  the  recursion  relationships 


P0(x)=  1, 

(61a) 

Pt(x)  =  x 

(61b) 

and 

P,(x)  = 

[2 -1/11  xPM(x)- [1-1/11  P,_ 

,2(x) 

(61c) 

[Ref.  21], 

Therefore, 

f2>1(0)  =  P,(COS  0)  . 

(62) 

From  equation  (57),  the  equation  describing  the  radial 

dependency  of  the  free 

i  wave  is 

{ 

d2  2  d  .  [fc2_  10*1) 

dr2  r  dr  r2 

(63) 
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Changing  scale  to  the  variable  s  -  kr  yields  the  differential  equation 


1  ds2  s  ds 


-  [1 


10+1) 

S2 


3}ti  =  o 


(64) 


[Ref.  17]. 


This  differential  equation  has  the  spherical  Bessel  function  j,(s) 
with  amplitude  A,  as  a  solution; 

J,  ,(s)  »  A,  j,(s)  (65) 

[Ref.  18). 


The  stationary  solution  of  the  free  wave  equation  in  spherical 
coordinates  is  now  written  as  the  partial  wave 

J,(r)  =  P,(cos  0)  [A,  j,(kr)] .  (66) 

The  most  general  form  of  this  free  wave  equation  solution  is  a  Fourier- 
Bessel  series, 

oo 

f(r,0)  =  2qA,  j,(kr)  P,(cos  0)  (67) 


[Refs.  17  and  21], 

A  plane,  free  wave  of  unit  amplitude  propagating  in  the  direction 
00=  0  is  then 

OO 

eiKZ  =  2^,  j,(kr)  P,(cos  0)  .  (68) 
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The  norm  of  the  Legendre  polynomials,  the  asymptotic  expression 

(69) 


sin  (kr-  $  1) 
),<*r> - jH- 


and  the  process  detailed  in  Elton  [Ref.  18]  are  used  to  determine  A,, 

Aj  =  i1  [21  +!]  .  (70) 


Therefore, 

e,kz  -  2M1  [21  ♦ 1  ]  jj(kr)  P^cos  0)  (71) 

[Ref's.  1 8  and  2 1  ]. 


3.  Detection  fonriitinns  and  Assumptions 

As  the  dimensions  of  the  incident  acoustic  wave  packet  are 
greater  than  the  size  of  the  target  it  is  expected  that  a  portion  of  the 
incident  wave  packet  continues  to  propagate  past  the  target.  The  total  or 
disturbed  wave  description  is  then  a  superposition  of  the  incident  and 
scattered  wave  descriptions.  [Ref.  1 7] 

The  acoustic  array  serves  as  both  transmitter  and  receiver,  or 
detector.  This  keeps  the  detector  outside  the  path  of  continued 
propagation  of  the  incident  wave  packet.  [Ref.  17]  Therefore,  only  the 
scattered  wave  description  is  of  interest  in  the  calibration  of  the  array. 
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The  placement  of  the  target  at  a  far  field  range  simultaneously 
places  the  detector  in  the  region  of  asymptotic  behavior  of  the  scattered 
wave. 

4  The  Disturbed  Wave 

Once  an  incident  wave  described  by  the  velocity  potential  l/(r) 
enters  the  region  of  influence  of  a  scattering  target's  spherical  potential 
V(r)  the  wave  equation  becomes  the  disturbed  wave  equation 


{  V2  +  [k2-  V(r)]}  K(r)  =  0  (72) 


[Ref.  1 7],  Comparison  of  the  form  of  the  disturbed  wave  equation  to  the 
free  wave  equation  gives  an  expectation  of  a  solution  of  the  form 

l/(r)  *  K,(r)  K2(b)  .  (73) 

Substituting  K(r)  into  the  disturbed  wave  equation  and 
separating  variables  yields 


1  d  r2d 


l^sin  0  d0 


d0 


The  angular  equation  for  the  disturbed  wave  is  identical  in  form 
to  the  angular  equation  for  the  free  wave  and  independent  of  the  Neumann 
boundary  condition  d  V/br  =  0.  Therefore, 


2  =  P](C0S  0) 


(75) 
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Using  the  previous  substitution  of  s  =  kr  yields  the  now  modified 


Bessel  equation  as  the  radial  equation  for  the  disturbed  wave; 


2_.d-4  [  ]  -  10+  =  o  (76) 

s  ds  s2  J  1 


Following  the  solution  development  in  Elton  [Ref.  18]  and 
applying  the  Neumann  boundary  condition  at  r  =  a  leads  to  a  solution 
l/]  ,(r)  =  i1  [21  +  1  ]  [jj(kr)  -  a,’  h/ 1  >(kr)]  .  (77) 

The  relationship 

*  j>a) 

a,'  =  i  e',fii  sin  5,  =  — — n —  (78) 

h  (ka) 

is  defined  from  Morse  [Ref.  17]  and  Bowman  [Ref.  16].  Additionally, 

•  the  j,(kr)  are  the  spherical  Bessel  functions  of  the  first  kind, 
argument  kr,  and  order  1, 

•  h,(1Hkr)  are  the  Hankel  functions  of  argument  kr  and  order  I, 

•  the  fi,  are  the  phase  shifts  of  the  usual  partial  wave  description  for 
spherical  scattering  solutions  [Refs.  17  and  20], 

•  the  j,'(ka)  are  the  first  derivatives  of  the  j,(kr),  taken  with  respect 
to  the  argument  kr  and  evaluated  at  r  =  a, 

•  h,(1)  (ka) »  j,'(ka)  +  i  y,'(ka),  and  (79) 

•  the  y,*(ka)  are  the  first  derivatives  of  the  y,(kr),  the  spherical  Bessel 

functions  of  the  second  kind  with  argument  kr  and  order  1,  taken  with 
respect  to  the  argument  kr  and  evaluated  at  r  =  a. 
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The  general  series  solution  for  the  total  disturbed  wave  Kr)  is  then 

r)=  jzQil  [21  +  l]  P^cos  0)  [jj(kr)  -  h,(1)(kr)]  .  (80) 

5.  The  Scattered  Wave  and  its  Amplitude 

The  general  asymptotic  form  for  a  disturbed  wave  Kr, 0) 

resulting  from  a  plane  wave  1/^,0)  incident  along  00  =  0  and  scattered  by 
a  spherical  potential  V(r)  is 

l/(r,0)  =  elkz  +  s(0)  eikr/r  (81) 

[Ref.  1 8],  This  description  of  the  asymptotic  disturbed  wave  is  a 
superposition  of  the  incident  plane  wave  eikz  and  a  spherical,  scattered 
wave  s(0)eikr/r.  The  term  s(0)  is  called  the  scattering  amplitude 
[Ref’s.  18,  19  and  221. 

Now,  recalling  the  series  equivalence  for  the  plane  wave  incident 
along  O0  in  equation  (71 ),  and  comparing  the  general  series  solution  for 

the  disturbed  wave  in  equation  (80)  to  its  asymptotic  form  in  equation 
(8 1 )  gives 

s(Q)eikr/r  =  -  2  i1  [21  +  1  ]  P,(cos  0)  a,’  h S 1  )(kr) .  (82) 

1=0  1  11 
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This  is  a  series  expression  for  the  scattered  wave.  Therefore, 

s(O)  =  — —  i  2  [2!  ♦  i]  Pi(cos  0)  a,’  (83) 

k  i=o  1  1 

[Ref's.  1 6  and  1 8]. 

6.  The  Backscattered  Cross  Section 

The  scattering  amplitude  s(O)  is  related  to  the  differential 
scattering  cross  section  <r(Q)  =  <rAfgt  by 

<r(0)  =  |s(0)|2  (84) 

[Ref's.  18,  19  and  22], 

Applying  the  series  expression  for  the  scattering  amplitude  in 
equation  (83)  to  the  previous  equation  yields 

12  i  [21  +  l]  a,'  P^cos  o)|2  .  (85) 

1=0  1  1 

0  =  Tl  this  gives  a  backscattered  cross 

section  of 

fftAgt  "  ffB  =  TT  [21  +  1  ]  ai' |2  <86a> 

[Ref's.  16  and  17],  For  the  spherical  target  Algt  =  Tta2,  so 

i  ln[2l*  t]a'l2.  (86b) 

6  (ko)2  '=0  ' 


4xf 


'V=i(°)*  — 


For  the  backscattered  direction 
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It  should  be  noted  a  plane  wave  of  unit  amplitude  incident  from 
the  positive  z  direction  and  backscattered  along  0  =  0  would  introduce  a 
factor  of  (-1 )'  into  the  previous  sum  [Ref.  16]. 

E.  SPHERICAL  BESSEL  FUNCTIONS  AND  THEIR  DERIVATIVES 

Implementation  of  the  previous  equation  necessitates  calculation  of 
the  a,’.  The  recursion  relations 

j0‘(ka)  «  [-j0(ka)  +  cos(ka)]/ka ,  (87a) 

jj’fka)  *  {-[l  ♦  l][j,(ka)]/ka}  +  jM(ka) ,  (87b) 

y0'(ka)  =  t-y0(ka)  +  sin(ka)]/ka  (87c) 

and 

y,'(ka)  =  {-[l  +  1  ][y  j(ka)]/ka}  +  y ,(ka)  (87d) 

determine  the  j,'(ka)  and  y,*  (ka)  needed  for  the  calculation  of  a(‘  in 
equation  (78)  [Refs.  21  and  23], 

Determination  of  the  regular  and  irregular  spherical  Bessel  functions 
at  r  =  a  supports  the  use  of  the  recursion  relations  in  equations  (87).  As 
ka  is  real  the  irregular  spherical  Bessel  functions  are  stable  enough  to  use 
the  recursion  relations 
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y0(ka)  =  -[cos(ka)]/ka , 
y  ,(ka)  =  [y0(ka)  -  sin(ka)]/ka 


(88a) 


(88b) 

and 

y,(ka)  =  {[21  -  1  ][y ,(ka)]/ka }  -  y  ,_2(ka)  (88c) 

in  their  determination  [Refs.  23  and  241 

Because  the  size  of  the  argument  of  the  regular  spherical  Bessel 
functions  may  be  very  small  the  method  of  negative  continued  fractions  is 
appropriate  to  the  fast  and  accurate  determination  of  a  converging 
solution  to  these  functions  [Ref.24],  The  regular  spherical  Bessel  function 
reciprocals  are  therefore  determined  by  the  negative  continued  fraction 
calculation 

_ =  1  [~3[bt  tb1  2,b1  jl..  [b, j3  (89) 

jj(ka)  j0(ka)  [bj  21 . Ib]  2]  [bj  jJbj  2]... 

Now, 

b)  m  =  2(1  +  m  -  l  )/ka  (90) 
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and 


_J _ 

bl,ttt-1 


1 


=  b 


_ L 

bl,m-1 


1 _ 

bl,m-2  1 - 

bl,m-3 


(90 


(Ref.  25] 

While  the  regular  spherical  Bessel  function  reciprocal  calculation 
of  equation  (89)  would  appear  to  involve  an  infinite  number  of  terms,  the 
calculation  can  in  fact  be  terminated  when  the  m1*  numerator  term  equals 
the  m01  denominator  term  to  the  desired  precision  for  the  spherical  Bessel 
function's  solution.  [Ref's.  24  and  25] 
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Small  quantities  of  data  were  collected  in  the  anechoic  chamber  at 
the  Naval  Postgraduate  School.  This  data  aided  first  order  verification  of 
the  calibration  process  and  its  supporting  software.  A  brief  description 
of  the  anechoic  chamber  is  provided  in  Appendix  A. 

The  calibration  process  and  data  development  occured  in  two  stages. 
The  first  stage  used  equation  (26)  and  two  way  propagation  path 
measurements  to  determine  <rbErEtG0.  The  second  stage  determined  EtG0 
using  equation  ( 1 7)  and  one  way  propagation  path  measurements. 

Er  was  then  determined  from  <rbErEtG0,  EtG0  and  the  calculated  <rb  of 

equation  (86b). 

Because  reciprocity  was  assumed  for  Er  and  Et 

ErEt  =  Er2  (92) 

provided  the  desired  calibration  information. 

Supporting  software  was  written  to  perform  calculations  required  for 
calibration.  The  next  section  describes  the  development  of  this  software. 


A.  CALIBRATION  SOFTWARE  DEVELOPMENT 

Calibration  software  was  written  in  BASIC.  It  was  used  with  a 
BASIC  4.0  compiler  on  a  Hewlett  Packard  3 10  computer.  Computer 
performance  was  enhanced  by  an  Infotek  floating  point  processor. 

Notation  in  the  programs  sometimes  differs  from  that  in  the  thesis 
text  as  a  result  of  H.P.  keyboard  limitations.  Variable  definitions  are 
located  at  the  beginning  of  subprograms  or  modules.  These  definitions 
are  provided  throughout  a  program  as  an  aid  in  determining  program  to  text 
variable  equivalence. 

The  main  effort  of  software  development  was  calculation  of 
backscattered  cross  sections.  This  was  done  with  equation  (86b). 
Information  available  for  summation  loop  verification  was  differential 
scattering  cross  section  information.  Differential  scattering  section 
software  was  developed  first. 

I.  Differential  Scattering  Cross  Section  Software  Development 
The  differential  scattering  cross  sections  <r  were  calculated 
using  equation  (85)  and  equation  (78). 

To  support  equation  (78)  employment  the  subprograms 
calculating  the  values  of  the  Legendre  polynomials  and  the  values  of  the 
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regular  and  irregular  spherical  Bessel  functions  were  developed  first. 

Each  subprogram  was  developed  and  tested  separately  in  test  programs. 
Each  test  program  was  designed  to  provide  a  tabular  output  for 
verification  against  values  found  in  Abramowitz  and  Stegun  [Ref.  23]. 

The  Legendre  values  were  calculated  using  equations  (61 ).  The 
irregular  spherical  Bessel  functions  were  calculated  using  equations  (88). 
The  regular  spherical  Bessel  functions  were  calculated  using  equations 
(89),  (90)  and  (91). 

Following  verification  these  subprograms  were  assembled  with 
subprograms  calculating  the  derivatives  of  the  regular  and  irregular 
spherical  Bessel  functions,  the  normalized  differential  scattering  cross 
section  <r,  the  scattering  modulus  of  Bowman  [Ref.  16],  and  with  three 
graphics  subprograms.  This  assembly  created  the  differential  scattering 
cross  section  program  found  in  Appendix  B. 

The  values  of  the  regular  and  irregular  spherical  Bessel  function 
derivatives  at  r=a  were  calculated  using  equations  (87). 

Differential  scattering  cross  sections  were  calculated  using 
equation  (85).  The  number  of  orders  1  required  for  the  summation 
calculation  in  equation  (85)  was  determined  by  generating  tabular  outputs 
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and  determining  the  number  of  orders  required  for  numerical  stability  to 
occur  to  12  decimal  places.  This  stability  process  was  based  on  the 
concept  of  the  phase  shift  8,  of  each  outgoing  partial  scattered  wave 

becoming  progressively  smaller  with  increasing  1  until  a,'  (equation  (78)) 

became  negligible  for  the  precision  desired  [Ref's.  18,  19,  20  and  22]. 

Tabular  outputs  covered  a  ka  value  range  of  one  to  thirty  at 
varying  scattering  angles,  including  0  =  rr.  This  ka  value  range  selection 
was  based  on  the  expected  maximum  target  size  to  be  encountered  during 
actual  echosounder  employment. 

Twenty-one  orders  greater  than  the  integer  value  of  the 
argument  ka  were  determined  to  be  required  for  equation  (85)’s 
summation.  The  number  of  orders  required  was  greater  than  the  usually 
cited  relationship 

ka  ~  /1(1  +  1)  (93) 

gives.  It  should  be  noted  that  the  precision  criteria  of  <r  was  greater  than 
that  normally  encountered.  [Ref's.  18,  19,  20  and  22] 

The  graphics  subprograms  of  the  differential  scattering 
cross  section  program  provided  outputs  for  comparisons  to  Figure  80  in 
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Morse  [Ref.  17]  and  Figure  10.10  of  Bowman  [Ref.  16].  These  comparisons 
were  used  to  verify  the  functioning  and  output  of  the  differential 
scattering  cross  section  program. 

The  graphics  subprograms  used  to  generate  Figure  4  provide 
linear  polar  plots  of  the  differential  scattering  cross  section  <r(0) 
normalized  by  the  cross  sectional  area  of  the  target  sphere.  These 

linear  polar  plots  show  the  normalized  differential  scattering  cross 
section  of  acoustically  hard  spheres  ensonified  by  plane  acoustic  waves 
incident  from  the  left,  i.e.  from  0  =  rr. 

The  first  view  of  each  sheet  in  Figure  4  shows  the  entire 
normalized  differential  scattering  cross  section  pattern,  including  the 
forward  lobe,  or  “shadow  zone",  along  0  =  0.  The  second  view  of  each 
sheet  is  a  magnified  plot  showing  the  scattering  detail  at  a  scale 
comparable  to  the  normalized  backscattered  cross  section  value.  Note 
the  two  views  are  identical  on  Sheet  1.  The  need  for  magnification  with 
increasing  ka  becomes  apparent  on  subsequent  sheets. 
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Figure  A,  Sheet  1  of  3.  Differential  Scattering  Cross  Section 
Program  Output  for  Comparison  to  Figure  80  of  Morse  [Ref.  1 7], 
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Figure  4,  Sheet  2  of  3:  Differential  Scattering  Cross  Section 
Program  Output  for  Comparison  to  Figure  80  of  Morse  [Ref.  17]. 
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Figure  4,  Sheet  3  of  3:  Differential  Scattering  Cross  Section 
Program  Output  for  Comparison  to  Figure  80  of  Morse  [Ref.  17], 
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Figure  4  was  used  for  comparison  to  Figure  80  of  Morse  [Ref.  1 7] 
The  radial  scale  of  Morse  was  not  graduated  so  verification  consisted  of 
a  check  for  the  correct  number  and  approximate  angular  location  of 
minimums  and  maximums  in  the  differential  scattering  cross  section 
patterns. 

The  graphics  subprograms  used  to  generate  Figure  5  provide 
semi -logarithmic  plots  of  the  scattering  modulus  5  of  Bowman  [Ref.  16]. 

5  =  k  s(a)  /  /4 tt  (94) 

is  the  scattering  modulus  5  of  Bowman  [Ref.  1 6].  s(a)  is  the  scattering 
amplitude  s(0),  as  defined  by  equation  (83),  with  a  factor  of  (-1)' 
introduced  into  the  sum,  that  is, 

s(e)  =  ~-i  2  (-1)1  (21  +  1)  P,(cos  0)  a  '  (95) 

k  i=o  11 

s(0)  is  normally  used  to  calculate  <r(Q),  per  equation  (84).  The  scattering 
modulus  5  of  Bowman  [Ref.  16]  could  also  be  used  to  calculate  <r(0).  Note 
the  introduction  of  (- 1 )'  into  the  summation  corresponds  to  the  case  of  an 
acoustically  hard  sphere  ensonified  by  a  plane  acoustic  wave  incident  from 
0  =  0. 
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Comparison  of  Figure  5  to  Figure  10.10  of  Bowman  [Ref.  16] 


shows  a  match  in  both  0  and  5.  This  graphical  match  verifies  the 
functioning  and  output  of  the  differential  scattering  cross  section 
program  of  Appendix  B. 


C  =  343  M/S. 

F  =  5000  HZ. 

A  =  .0109180290961  M. 

K*A  =  1  . 

NORMALIZED  BCKSCTR  =  .880073464081 

Figure  5,  Sheet  1  of  3:  Differential  Scattering  Cross  Section 
Program  Output  for  Comparison  to  Figure  10.10  of  Bowman  [Ref.  16]. 


A 1 


CTRG  MODULUS 


LQ 


m 


0  30  60  20  120  150  180 


PHI  c DEGREES) 


C  =  343  M/S. 

F  =  5000  HZ. 

A  =  .0545901454805  M . 

K*A  =  5  . 

NORMALIZED  BCKSCTR  =  1.08120386888 


Figure  5,  Sheet  2  of  3:  Differential  Scattering  Cross  Section 
Program  Output  for  Comparison  to  Figure  10. 10  of  Bowman  [Ref.  16], 
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2.  Backscattered  Cross  Section  Software  Development 

The  differential  scattering  cross  section  program  of  Appendix  B 
was  modified  to  yield  the  backscattered  cross  section  program  of 
Appendix  C.  With  the  modification  the  normalized  differential  scattering 
cross  section  is  calculated  for  the  backscattered  direction  only.  Equation 
(86b)  is  used  for  this  calculation. 

Additionally,  the  graphics  subprograms  of  the  differential 
scattering  cross  section  program  were  replaced  by  one  graphics 
subprogram.  The  new  graphics  subprogram  was  designed  to  provide  a 
rectangular  graph  of  the  normalized  backscattered  cross  section,  <rb. 

<rb  is  plotted  as  a  function  of  ka  for  a  specifiable  range  of  values  of  ka. 

Figure  6  was  generated  by  the  backscattered  cross  section 
program.  It  was  used  for  comparison  to  Figure  10. 1 1  of  Bowman  [Ref.  16]. 
This  comparison  served  as  a  verification  of  the  calculated,  normalized 
backscattered  cross  section  values.  Note  the  ka  axis  is  labeled  K*A. 

Figure  7  was  generated  to  verify  the  normalized  backscattered 
cross  section  profile  was  approaching  a  value  of  one  as  ka  left  the  Mie 
scattering  regime  and  entered  the  "geometric  optics"  region.  Note  the  ka 
axis  is  labeled  K*A. 
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Figure  6:  Backscattered  Cross  Section  Program  Output 
for  Verification  of  Calculated  Normalized  Backscattered 
Cross  Section  Values  by  Comparison  to  Figure  10. 1 1  of 
Bowman  [Ref.  16], 
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Figure  7:  Backscattered  Cross  Section  Program  Output 
for  Verification  of  the  Normalized  Backscattered  Cross 
Section  Asymptotic  Behavior  with  Increasing  ka. 


3.  Acoustic  j^rray  Efficiency  Software  Development 

The  acoustic  array  efficiency  program  in  Appendix  D  determined 
Er  Et.  In  general,  each  major  section  of  Chapter  1 1  of  this  thesis  and  each 

set  of  analyzed  data  is  represented  by  a  subprogram  of  the  acoustic  array 
efficiency  program. 

The  acoustic  array  efficiency  program  of  Appendix  D  includes  the 
non-graphics  subprograms  present  in  the  backscattered  cross  section 
program  of  Appendix  C.  It  also  includes  subprograms  designed  to  furnish 
stored  data,  calculate  attenuation  and  range  using  equations  (27)  through 
(42),  and  use  developed  information  to  calculate  EtG0  with  equation  (17). 

Equation  (86b)  was  used  to  calculate  <rb. 

The  appropriate  subprogram  and  main  program  calculational 
results  were  then  used  in  equation  (26)  to  determine  Ep.  Er  was  used  in 

equation  (92)  to  determine  Er  Et 

B.  DATA  ACQUISITION  AND  SUPPORTING  EQUIPMENT 

Data  collection  occured  28  November  1987,  30  November  1987,  and 
7  December  1987. 
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Data  collected  28  November  and  30  November  provided  information 
needed  for  calculation  of  <r.E„E.Gn.  Data  collected  7  December  allowed 

brio 

calculation  of  E,Gn  and  the  calculation  of  E 's. 

to  r 

Collected  data  is  presented  in  Appendix  E. 

i.  Two  Way  Propagation  Path  Measurements 

Two  way  propagation  path  data  were  collected  28  and  30 
November.  Measurements  were  used  to  calculate  the  product  <rbErEtG0. 

Figure  8  provides  a  connection  schematic  for  the  data  acquisiton 
equipment. 


Figure  8:  Data  Acquisition  Equipment  Set-up 
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Four  target  spheres  were  available  and  measured  for  size.  They 
included  one  hollow  aluminum  sphere  (sphere  number  one)  and  three  solid 
aluminum  and  brass  spheres.  Table  ( 1 )  lists  the  data  and  the  calculated 
values  of  Algl  =  na2,  ka,  and  <rb.  <rb  was  calculated  from  equation  (86b). 

Table  1:  Target  Sphere  Dimensional  Data  and  Normalized 
Backscattered  Cross  Section  Calculationa!  Results  for 
Frequency  of  5000  [Hz]. 


Target 

No. 

Sphere 

Diameter  [m] 

Sphere 
Radius  [mj 

V  i™2! 

ka 

1 

.2546  +.0006 

.1273 

.05091 

1 1.62 

.9142 

2 

.1013  +.0001 

.05065 

.008060 

4.619 

1.000 

3 

.076341.00003 

.03817 

.004577 

3.479 

1.051 

4 

.06240+00008 

.03120 

.003058 

2.843 

.6836 

Figure  9  shows  the  target  spheres’  theoretical  normalized 
backscattered  cross  section  distribution  plotted  on  a  graph  generated  by 
the  backscattered  cross  section  program  of  Appendix  C. 
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Figure  9:  Target  Spheres'  Normalized  Backscattered  Cross 
Section  Distribution  as  Calculated  from  Measured  Diameters, 
Using  an  Assumption  of  Perfect  Acoustic  Hardness. 


The  target  spheres  were  suspended  in  turn  in  the  anechoic 
chamber.  Placement  was  made  using  the  criteria  in  Section  D.  1  of 
Chapter  II.  The  range  R  from  the  array  to  a  target  was  determined  to  be 
5.30  [ml. 

Temperature  Tc  and  relative  humidity  Rh  along  the  propagation 
path  were  determined  with  a  Weathermeasure  Temperature  and  Relative 
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Humidity  Meter.  T.  was  usually  about  21  [°c]  with  a  relative  humidity  of 
approximately  50%.  Atmospheric  pressure  P  was  read  from  a  barogram 

supplied  by  the  Naval  Postgraduate  School  meteorological  station. 

Specific  Tc,  Rh,  and  P  data  are  furnished  in  Appendix  E. 

The  unshrouded,  close  packed,  hexagonal  acoustic  array  of 
Moxcey  [Ref.  8]  was  used  to  ensonify  each  target. 

A  =  Ns  it  rs2  (96) 

determined  the  aperture  area  A  of  the  array.  Ns  =  1 9  was  the  number  of 
speakers  in  the  array  and  rs  =  .0381  [m]  was  the  average  speaker  horn 
radius  [Ref.  8]. 

Moxcey  [Ref.  8]  determined  the  half  width  of  the  array's  main 
lobe  to  be  12°.  The  angular  half  widths  of  targets  two,  three  and  four 
were  all  less  than  5%  of  the  main  lobe’s  half  width.  Target  one's  angular 
half  width  was  1 1.5%  of  the  main  lobe's  half  width. 

A  sinusoidal  input  from  a  Hewlett  Packard  33 1 4  A  function 
generator  drove  the  array  via  a  pre-amplifier.  Signal  input  from  the 
HP33 1 4A  function  generator  to  the  pre-amplifier  was  set  for  20  cycles  at 
5000  [Hz].  The  voltage  input  to  the  pre-amplifier's  power  amplifier 
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terminal  (see  Figure  8)  was  3.5  [V]  for  spheres  two,  three  and  four  and  3.0 
[V]  for  sphere  one. 

Twenty  cycles  is  considerably  less  than  the  number  of  cycles 
expected  in  actual  echosounder  application.  It  was  necessary  to  keep  the 
pulse  packet  this  small  because  the  anechoic  chamber's  corner  was  only 
about  1 .5  [m]  from  the  target  mounting  location.  Voltage  input  for  sphere 
one  was  reduced  to  3.0  [V]  to  avoid  clipping  of  its  greater  return  signal 
voltage. 

Actual  transmission  voltage  inputs  Vt  to  the  array  from  the 

pre-amplifier  were  measured  with  a  Nicolet  3091  oscilloscope  at  the 
pre-amplifier's  array  terminal.  The  array  was  disconnected  from  the 
pre-amplifier  for  Vt  measurements.  Vt  values  are  provided  in  Appendix  E. 

They  are  also  stored  in  the  Trans30  and  Trans35  subprograms  of  the 
acoustic  array  efficiency  program  in  Appendix  D. 

The  Nicolet  oscilloscope  was  connected  to  a  Hewlett  Packard 
7090  A  measurement  plotting  system.  This  connection  provided  a  hard 
copy  plot  of  voltage  vs.  time  traces. 

The  returned  signal  was  amplified  and  filtered. 
Pre-amplifier  gain  Ge  was  determined  asl  1,094  by  Moxcey  [Ref.  8], 
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Filtering  was  done  with  a  Wavetek  Brickwall  filter  model  753  A  to  remove 
the  low  frequency  normal  modes  of  oscillation  of  the  speaker  drivers  from 
the  return  signal.  Filtering  was  accomplished  at  0  [dB]  gain  and  5000  [Hz] 
for  both  low  and  high  pass  settings. 

The  processed  return  signal  voltages  Vr  were  displayed  and 

measured  on  the  Nicolet  3091  oscilloscope.  Vp  measurements  were  made 

at  the  Out  2  terminal  of  the  Wavetek  Brickwall  filter  (see  Figure  8).  They 
were  taken  at  10  [pS]  intervals  for  two  full  cycles  along  the  asymptotic 
region  of  the  return  signal  trace. 

The  Nicolet  3091  oscilloscope  was  connected  to  the  Hewlett 
Packard  7090  A  measurement  plotting  system.  Figure  10  shows  the 
Hewlett  Packard  plotter's  traces  for  target  sphere  two's  return  signal 
voltages  as  printed  and  subsequently  labeled.  The  figure  is  representative 
of  the  oscilloscope  traces  for  each  target  return.  Data  for  all  targets  are 
provided  in  Appendix  E. 

The  time  of  return  t  was  determined  as  the  time  difference 

r 

between  the  first  positive  maxima  of  the  transmitted  pulse  packet  and  the 
first  positive  maxima  of  the  target's  return  pulse  packet. 
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Figure  10,  Sheet  1  of  2:  Target  Sphere 
Two's  5000  [Hz]  Return  Signal  Trace 
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Figure  10,  Sheet  2  of  2:  Target  Sphere  Two's 
5000  [Hz]  Return  Signal  Trace,  Scaled  Down 
For  Oscilloscope  Measurement. 
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Noise  voltage  Vn  measurements  were  made  with  the  same 
measuring  equipment  setup  used  for  Vp  measurements.  The  target  spheres 
were  removed.  V  measurements  were  taken  at  a  time  of  return  t 

n  r 

comparable  to  the  time  of  measurement  of  V 

Vn  measurements  were  actually  made  on  occasions  well 
separated  chronologically  from  the  Vr  measurements.  Therefore,  they 

were  proportionally  corrected  in  line  2530  of  the  acoustic  array 
efficiency  program.  The  corrections  account  for  the  slight  variations  in 
input  voltages  and  propagation  path  attenuations  that  occurred  between 
the  two  measurement  occasions. 

Figure  1 1  shows  a  representative  noise  trace.  Data  from  the 
noise  measurements  is  presented  in  Appendix  E  and  stored  in  the  Noise30 
and  Noise35  subprograms  of  the  acoustic  array  efficiency  program  in 
Appendix  D. 

Unless  otherwise  specified  above  the  remaining  data  from  the 
two  way  propagation  path  measurements  are  stored  in  the  subprograms 
titled  Sphere  1,  5phere2,  5phere3  and  Sphered.  All  remaining  data  are  also 
presented  in  Appendix  E. 
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Figure  1 1:  Representative  5000  [Hz]  Noise  Trace. 
3.0  [V]  Supplied  by  the  HP3314A  Function  Generator. 


2.  One  Way  Propagation  Path  Measurements 

One  way  propagation  path  data  were  collected  7  December. 
Measurements  were  used  to  calculate  EtGo  with  equation  (17). 

Data  from  the  one  way  propagation  path  measurements  are 
provided  in  Appendix  E.  Data  are  also  stored  in  the  Gain30  and  Gain35 
subprograms  of  the  acoustic  array  efficiency  program  of  Appendix  D. 
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For  one  way  path  measurements,  i.e.  for  the  determination  of 
E,G0,  an  Ivie  Electronics  Inc.  IE  30  A  Audio  Analyzer  (Serial  *805  A  426) 

was  used  to  measure  received  target  intensity  lr  in  [dB].  Intensity 

measurements  were  referenced  to  10'12  [W/m2].  The  acoustic  signal 
supplied  for  this  measurement  was  a  continuous  wave  at  5000  [Hz], 
Measurements  were  taken  at  the  suspension  site  of  the  calibrating  targets. 

A  Fluke  8060  A  True  RMS  Multimeter  was  used  to  measure  the 
root-mean-square  voltage  (Vrms)t  supplied  to  the  acoustic  array  at  its 

input.  (Vpms)t  was  squared  to  determine  (Vms)t. 

The  Weathermeasure  temperature  and  relative  humidity  meter 
was  unavailable  for  EtG0  data  collection.  Hygrothermographs  from  the 

Naval  Postgraduate  School's  meteorological  station  were  used  tc  estimate 
the  anechoic  chamber's  relative  humidity  Rh  and  temperature  Tc  on 

7  December.  Estimates  were  based  on  comparisons  to  the 
hygrothermographs  and  measured  Rh  and  Tc  of  28  and  30  November. 
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C.  RE5ULT5 


Data  analysis  was  performed  with  the  aid  of  the  acoustic  array 
efficiency  program  of  Appendix  D.  Data  used  in  calculations  is  furnished 
in  Appendix  E. 

The  product  ffbErE,G0  was  generated  from  two  way  popagation  path 
data  and  application  of  equation  (26).  The  average  of  calculated <rb's 
compared  to  the  average  of  the  experimental  product  <rbErEtG0  gives  a 
normalizing  product  ErEtG0.  The  experimental  product  «-bE  EtG  and  the 
normalizing  product  ErEtG0  are  used  to  calculate  experimental  <rb's. 
These  results  are  supplied  in  Table  2. 


Table  2.  Target  Sphere  Experimental,  Normalized 
Backscattered  Cross  Section  Calculational  Results. 


Target 

No. 

ka  [m2] 

Calc'd 

ffb 

E  E.G 

r  i  o 

Exp’l 

ffb 

1 

1 1.62 

.9142 

55.90 

.947 

2 

4.619 

1.000 

59.88 

1.01 

3 

3.479 

1.051 

61.62 

1.04 

4 

2843 

.6836 

37.98 

.644 

Avg. 

.912 

53.85 

59.0 
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The  experimental  <rb  are  plotted  on  a  graph  generated  by  the 
backscattered  cross  section  program  of  Appendix  C.  This  plot  is  Figure 
12.  It  is  similar  to  Figure  9  and  gives  a-b  as  a  function  of  ka.  Note  ka  is 


labeled  K*A. 


K*fi 


Figure  12:  Target  Spheres'  Experimental,  Normalized 
Backscattered  Cross  Section  Distribution  Compared  to 
a  Calculated  Backscattered  Cross  Section  Distribution 
Curve  Generated  by  the  Backscattered  Cross  Section 
Program  of  Appendix  C. 
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Sphere  1 


As  ka  increases  the  experimental  <rb  distribution  rises  from  just 
below  the  calculated  <rb  curve  to  just  above  it. 

EtGo  was  determined  using  one  way  propagation  path  data  and  equation 
(17).  EtG0,  or  Et>  was  assumed  to  be  the  same  for  spheres  two,  three,  and 

four  because  the  supply  voltage  and  target  locations  were  the  same.  The 
angular  half  widths  of  target  spheres  two,  three  and  four  were  between 
4.6%  and  2.8%  of  the  echosounder's  main  lobe's  half  width. 

Table  3:  Acoustic  Array  Product  EtGo 

Determined  from  One  Way  Propagation 
Path  Measurements  and  Equation  ( 1 7). 


Target 

No. 

HP3314A 

Func'n  Gen'r 
Supply  Volt.  [V] 

EtGo 

1 

3.0 

94.70 

2 

3.5 

95.1 1 

3 

3.5 

95.1 1 

4 

3.5 

95.1  1 

Calculated  <rb,  EtGo  and  two  way  propagation  path  data  were  used  in 
equation  (26)  to  determine  Ef.  «rb  values  were  calculated  from  the  target 
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spheres'  measured  diameters,  2a.  A,gt  was  calculated  as  Tta2.  Calculation 
results  are  furnished  in  Table  4. 

Table  4:  Er  and  EtEr  Er  Determined  From  Two  Way  Propagation 
Path  Data  and  Equation  (26).  EtEr  Determined  From  Assumption 
of  Efficiency  Reciprocity. 


Target 

No. 

V  im2i 

Calc’d 

ffb 

V 

Er 

¥r 

1 

.05091 

.9142 

.04654 

.6456 

.4169 

2 

.008060 

1.000 

.008060 

.6296 

.3964 

3 

.004577 

1.051 

.004810 

.6167 

.3803 

4 

.003058 

.6836 

.002090 

.5841 

.341 1 

The  Er's  determined  for  Moxcey's  unshrouded  hexagonal  array  are  in 

the  range  of  .58  to  .65.  This  is  above  the  value  range  of  0.5  +  0. 1 
determined  by  Weingartner  for  his  square  acoustic  array  [Ref.  6} 

Table  4  results  were  used  to  construct  Figures  13  and  14.  Figure  13 
gives  Er  as  a  function  of  the  calculated,  normalized  backscattered  cross 

section,  <rb  Figure  14  gives  Er  as  a  function  of  the  apparent  target  size, 
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Target  Number 

4  I 


Figure  13:  Acoustical  to  Electrical  Power  Conversion  Efficiencies, 

E  as  a  Function  of  the  Normalized  Backscattered  Cross  Sections, 
r « 

<rb.  Calculated  from  Measured  Sphere  Sizes. 


Target  Number 

4  3  2  1 


o'bNt lrf|21 


Figure  14:  Acoustical  to  Electrical  Power  Conversion  Efficiencies, 
Er  as  a  Function  of  the  Apparent  Target  Sizes,  <rbAlgl. 
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Recalling  equation  (20), 


v[S,]KKRgH [-^i,  <2o> 


and  equation  (12), 


p  =  W 
ra  '  4n'R 


2  e  A  , 


(12) 


serves  as  a  reminder  that  the  received  acoustic  power,  Pra,  depends  on  the 
apparent  target  size,  o^A^. 

Figure  1 4's  trend  and  the  dependency  of  Pra  on  ffbAtgl  were  used  as 
suggestions  for  the  construction  of  "Figure  15.  Figure  15  is  a  plot  of  Er 


Target  Humber 

4  3  2 


1 


Pl*8 


[WJ 


Figure  15.  Acoustical  to  Electrical  Power  Conversion  Efficiencies, 
Ep  as  a  Function  of  the  Received  Acoustic  Powers,  Prg. 
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versus  P 


Pra  was  calculated  using  the  definition  (from  equation  (21 )) 


pra  -  pr  '  Er 


and  the  relationship  (from  equation  (15)) 


Pr  '  <VmsV  /  Z  • 


Note  Figure  15‘s  trend  is  for  Er  to  increase  with  an  increase  in  Pra  for 


the  range  of  target  sizes  investigated. 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  larger  magnitudes  of  Er  determined  for  Moxcey's  hexagonal  array 

suggest  it  may  be  more  efficient  than  Weingartner's  square  array. 

However,  the  drivers  for  each  array  are  identical. 

The  0.08  to  0.15  differences  in  efficiencies  are  more  likely  the  result 
of  the  difference  in  calibration  methods.  Weingartner’s  measurements 
[Ref.  6]  used  steady  state,  continuous  waves.  The  use  of  steady  state, 
continuous  waves  allows  the  development  of  standing  waves.  This 
calibration  investigation  used  pulsed  waves.  The  use  of  pulsed  waves  did 
not  provide  sufficient  time  for  standing  waves  to  completely  develop  in 
the  backscattered  energy.  [Ref.  26]  This  calibration  method  should  provide 
more  accurate  efficiency  values  since  it  more  nearly  approximates  the 
echosounder  system's  actual  field  operating  configuration. 

Additionally,  steel  braid  and  nylon  lines  were  used  to  suspend  the 
target  spheres.  These  supports  probably  reflected  enough  acoustic  energy 
to  have  a  small  influence  on  results.  Other  fixtures  in  the  anechoic 
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chamber,  such  as  support  clamps  and  lights,  may  have  also  contributed 
small  amounts  of  reflected  energy  to  the  returned  signal. 

Calibration  of  the  computer  controlled  echosounder  required 
determination  of  the  product  ErEt.  The  calibration  process  assumed 

Er  ■  E}.  The  calibration  results  in  Tables  3  and  4  for  spheres  two,  three, 
and  four  indicate  Er  varied  with  constant  EtG0,  or  constant  Et.  Such  a 
variation  negates  the  assumption  of  reciprocity  for  Ep  and  Et 

Figure  15  suggests  the  acoustic  array  s  efficiency  Er  increased  in 
response  to  an  increase  in  Pra  and  (Vrms)r. 

The  acoustic  gain  for  spheres  two,  three  and  four  was  probably  not 
constant,  as  assumed.  However,  G0  is  the  maximum  for  G(O),  and  Er  is 

inversely  proportional  to  G0,  per  equation  (26).  An  averaging  process 

accounting  for  the  slight  differences  in  gain  encountered  by  spheres  two, 
three  and  four  would  have  increased  the  data  points  spread  observed  in 
Figures  14  and  15 

Additionally,  the  pre-amplifiers  gain  may  not  have  been  constant  for 
the  range  of  received  voltage  values  encountered.  This  would  have 
introduced  either  an  apparent  increase  in  Er  or  dampened  an  increase  in  Er 
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with  an  increase  in  (Vrms)r.  Such  a  gain  variation  could  have  effectively 

negated  the  assumption  of  efficiency  reciprocity.  However,  this 
pre-amplifier  is  approximately  0.1%  linear  for  less  than  10  [V]  peak 
[Ref.  26],  Received  voltages  in  this  investigation  were  all  less  than  10  [V]. 

B.  RECOMMENDATIONS 

Because  efficiency  reciprocity  is  now  in  question  it  is  necessary  to 
determine  G0  by  calculation  or  measurement.  Such  a  determination  would 

allow  the  desired  calibration  quantity  ErEt  to  be  written  as 

EtG0 

ErE.  =  Er  -if.  <"> 

One  method  of  calculating  G0  would  involve  integrating  the  acoustic 

array’s  calculated  or  measured  intensity  pattern  (see  Moxcey  [Ref.  81)  for 
determination  of  an  equivalent  isotropic  emitter.  To  determine  G0  the 

centroid  intensity  of  the  array's  emitted  intensity  pattern  would  be 
compared  to  that  of  the  equivalent  isotropic  emitter.  This  would  be 
repeated  for  each  transmission  power  level  of  interest. 
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Calculation  of  G0  would  allow  a  determination  of  Et  using  EtG0  .  Such 
a  determination  of  Et  would  allow  a  reciprocity  comparison  to  Er  and 
support  a  more  accurate  evaluation  of  ErEt. 

The  accuracy  of  Er  and  EtG0  could  be  further  refined  as  follows. 

Use  a  calibrating  site  that  allows  the  use  of  a  pulse  train  closer  in 
length  to  actual  field  use. 

Make  target  return  and  noise  voltages  and  site  intensity 
measurements  as  chronologically  close  as  possible. 

Measure  the  barometric  pressure  in  the  propagation  path  directly 
along  with  the  relative  humidity  and  temperature. 

The  received  and  transmitted  electrical  powers  for  Er  and  the 

transmitted  electrical  power  for  EtGQ  were  all  measured  at  the 

pre-amplifier  end  of  the  pre-amplifier/array  transmission  line.  This 
included  transmission  line  power  losses  due  to  impedance  as  efficiency 
decreases.  The  present  calibration  process  includes  the  transmission  line 
impedance  as  part  of  the  array  impedance. 

The  amount  this  impedance  power  loss  actually  affects  Er  could  be 

evaluated  by  determining  the  impedances  seen  in  each  direction  at  the 
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pre-amplifier's  array  terminal.  This  would  be  an  opportune  time  to 
verify  the  validity  of  transmission  line  impedance  reciprocity. 

Use  the  Nicolet  oscilloscope  to  record  a  triggered  readout  from  a 
calibrated  microphone  suspended  at  the  target  site  at  45°  to  the  centroid's 
axis.  This  would  allow  EtG0  to  be  determined  for  a  pulse  packet  rather 

than  for  a  continuous  wave.  Coupled  with  a  proper  calibrating  site 
selection  this  technique  should  minimize  interference  problems  such  as 
anechoic  chamber  corner  reflections. 

A  larger  number  of  apparent  target  sizes  spanning  the  desired  range 
of  apparent  target  sizes  should  be  used  to  redetermine  Er  for  Moxcey's 

array.  The  Er  distribution  should  be  examined  for  a  smooth  and  asymptotic 
behavior.  Comparative  Er  values  for  Moxcey's  array  could  be  determined 

using  the  method  outlined  by  Weingartner  [Ref.  6]. 

Weingartner's  square  array's  Er  could  also  be  determined  with  this 

thesis'  calibration  process  and  compared  to  Weingartner's  Er  value  [Ref.  6). 

The  refined  Er  curve  in  the  now  calibrated  echosounder  system  should 
be  used  to  develop  an  absolute  CT2  plot  of  an  air  mass.  An  established 
method  such  as  a  tower  mounted  vertical  array  of  thermocouples  would  be 


69 


used  simultaneously  to  generate  a  comparative  absolute  CT2  plot  of  the 
same  air  mass. 

C.  SUMMARY 

Realizing  the  potential  of  the  computer  controlled  echosounder  for 
analyzing  lower  atmospheric  turbulence  required  a  calibration  of  the 
echosounder.  This  thesis  described  theory  and  software  for  performing 
calculations  crucial  to  a  backscattered  cross  section  calibration  of  the 
echosounder.  The  backscattered  cross  section  calibration  was  evaluated 
in  this  thesis. 

Er's  values  indicate  the  evaluated  calibration  process  requires  some 

refinement.  They  also  indicate  the  calibration  method  possesses 
sufficient  merit  to  warrant  further  development  and  that  the  software 
performs  as  intended. 
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ANECHOIC  CHAMBER  DESCRIPTION 

The  anechoic  chamber  is  located  in  room  019  of  building  232  at  the 
Naval  Postgraduate  School.  Acoustic  research  on  sound  sources,  sound 
receivers  and  sound  scatterers  may  be  conducted  in  the  chamber  with  a 
minimum  of  interference  from  wall  reflection  and  external  noise. 

Wall  reflection  is  minimized  by  102  centimeter  deep  wedges  made  of 
P.F.  612  fiberglass.  The  wedges  are  attached  to  the  walls,  ceiling  and 
floor  of  the  chamber.  Approximately  142  cubic  meters  of  fiberglass  is 
used  to  absorb  99%  of  incident  sound  energy  for  frequencies  greater  than 
100  Hertz.  External  noise  isolation  is  accomplished  by  separating  the 
chamber's  inner  concrete  block  sides  from  the  outer  12  inch  concrete 
walls  with  a  two  inch  thick  lining  of  fiberglass  and  cork. 

A  usable  region  of  approximately  8.2  meters  by  4.3  meters  by  3.4 
meters  is  available.  This  region  is  floored  by  a  grid  of  225  wire  cables, 
with  each  cable  placed  in  a  tension  of  150  to  200  pounds  per  square  inch.  * 

*This  Appendix  was  compiled  from  the  description  posted  by  the 
entrance  to  the  anechoic  chamber  instrumentation  and  control  room. 


APEEND1X.-B 

DIFFERENTIAL  SCATTERING  CROSS  SECTION  PROGRAM 


Examples  of  the  output  of  this  program  are  Figures  4  and  5,  located  in 
Chapter  III. 
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170 
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190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 


REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
ACOUSTICALLY  HARO  SPHERE  SCATTERING  OF  PLANE  ACOUSTIC  WAVES 

THIS  PROGRAM  CA  CULATES  AND  GRAPHS  THE  DIFFERENTIAL  SCATTERING 
CROSS  SECTION  AND/OR  SCATTERING  MODULUS  FOR  AN  ACOUSTICALLY  HARD 
SPHERE  OF  DIAMETER  0  ENSONIFIED  BY  AN  INCIDENT  PLANE  WAUE  OF 
FREQUENCY  F. 


I++++  VARIABLE  DECLARATIONS  AND  DEFINITIONS  ++++++++++++++++++++++++++++ 


REAL 

C 

I 

j 

1  SPEEO  OF  SOUND. 

i 

REAL 

Freq 

!  FREQUENCY  OF  THE  EMITTER ! 

I 

REAL 

A 

!  RADIUS  OF  THE  SPHERE. 

1 

REAL 

L 

i 

!  ORDER  OF  THE  FUNCTION  CALCULATED 

1  IN  A  SUBROUTINE.  IT  IS  USED 
!  AS  A  LOOP  INDEX  AND  TO  DESIGNATE 
!  THE  ARRAY  ELEMENT  FOR  THE  FUNCTION 
1  OF  CORRESPONDING  ORDER. 

1 

REAL 

Lnax 

!  MAX.  SIGNIFICANT  ORDER  OF  THE  SUM 
!  DETERMINING  THE  DIFFERENTIAL 

1  SCATTERING  CROSS  SECTION  (DSCS) 

1  AND  SCATTERING  MODULUS  ( SM  ) . 

1 

REAL 

K 

1  THE  WAUE  NUMBER,  2*PI*FREQ/C. 

1 

REAL 

Ka 

! 

i  K • A ,  THE  ARGUMENT  OF  THE  REGULAR 

1  AND  IRREGULAR  SPHERICAL  BESSEL 
!  FUNCTIONS  AND  THEIR  DERIVATIVES. 

| 
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410  REAL  J( 51  ) 

420 

430 

440 

450 

460  REAL  D j ( 5 1  ) 

470 

480 

490  REAL  Y<  S 1  ) 

500 

510 

520 

530 


I  ARRAY  FOR  THE  REGULAR  SPHERICAL 
1  BESSEL  FUNCTIONS  <  RSB )  OF 
!  ARGUMENT  KA  AND  ORDER  0  THROUGH 
1  LMAX . 

1  ARRAY  FOR  THE  OERIUATIUES  OF  THE 
!  RSB  (THE  DRSB). 

!  ARRAY  FOR  THE  IRREGULAR  SPHERICAL 
I  BESSEL  FUNCTIONS  (ISB)  OF 
1  ARGUMENT  KA  ANO  ORDER  0  THROUGH 
1  LMAX. 

I 


540 

550 

560 

570 

580 

590 

600 

610 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

860 

870 

880 

890 

900 


REAL  Dy(51  ) 


ARRAY  FOR  THE  OERIUATIUES  OF  THE 
ISB  (THE  DISB  ) . 


REAL  P< 51  ) 


ARRAY  FOR  THE  LEGENDRE 
POLYNOMIALS  OF  ARGUMENT  COS(PHI). 


REAL  Sm( 360  ) 


THE  ARRAY  REPRESENTING  THE 
SCATTERING  MODULUS  ( SM  ) . 


REAL  Os( 360  ) 


THE  ARRAY  REPRESENTING  THE 
1  DIFFERENTIAL  SCATTERING  CROSS 
!  SECTION  NORMALIZED  BY  THE  CROSS 
!  SECTIONAL  AREA  OF  THE  SPHERE 
1  ( DSCS  ) . 

!  NOTE  DS( PHI  )=4*SM( PHI >“2/KAA2  . 


!  ++++  NOTE  ++  +  + 

1  THE  ARRAY  DIMENSIONS  (30+21)  NEED 
I  TO  BE  INCREASED  TO  ACCOMODATE 
1  KA>30. 


I 


H++++++++ +++++++++++ 


I +++++++++++++++++ +++++++++H 
! 

!++++  MAIN  PROGRAM  +++++++++++++++++++++++++++++++++++++++++++++++++++++ 


CALL  I  rti  t  (  C  ,Freq  ,  A  ,Lmax  ,Ka  ) 


INPUT  AND  CALCULATE  REQUIRED 
PARAMETERS. 


CALL  Difsctr(Lmax  ,  l\a  ,  u  j  ( *  )  ,  J  ( *  )  ,Dy( •  )  ,Y( •  )  ,P( *  >  ,Sn( *  >  ,Ds( *  )  ) 

1  CALCULATE  THE  SMC)  AND  DSC). 


CALL  Ou tp t ( C ,Freq  ,A  ,Ka  , S m( »  )  ,  Os  ( *  )  ) 


OUTPUT  THE  SMC)  AND/OR  DSC). 


i 

i 
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910  !+  +  +  +  MAIN  PROGRAM  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++ 

920  ' 

930  PRINTER  IS  I 
940  PRINT 

950  PRINT  "END"  !  ADUISES  THE  USER  OF  THE  MAIN 

9G0  PRINT  !  PROGRAM’S  CONCLUSION. 

970  END 

980  ! 

990  !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

1000  ! 

1010  ! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
1020  ! 
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1030 

1040 

1 

| $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

1050 

SUB  Imt(  C  ,Freq  ,A  ,Lnax  ,Ka  ) 

1060 

REM  ++++++++++++++ ++++++++++++++++++++++++++++*+++++++++++++++++++++++ 

1070 

REM  THIS  MODULE  REQUESTS  THE  APPROPRIATE  INPUT  VARIABLES  TO  DETERMINE 

1080 

REM  C,  FREQ,  A,  LMAX  ,  AND  KA. 

1080 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

1  100 

1 

1110 

!++++  MODULE  UARIABLE  DECLARATIONS 

AND  DEFINITIONS  +++++++++++++++++++++ 

1120 

J 

1130 

1 

1  140 

REAL  Tercp 

\ 

AMBIENT  AIR  TEMPERATURE  IN  DE6  C. 

1  150 

j 

1  160 

REAL  Circ 

1 

CIRCUMFERENCE  OF  THE  TGT .  SPHERE. 

1  170 

1 

1  180 

REAL  Di am 

| 

DIAMETER  OF  THE  TARGET  SPHERE. 

1  190 

! 

1200 

! 

1210 

! +++++++++++++++++++++++++++++++++++++++++++++++++++++*++++++++++++++++++ 

1220 

1 

1230 

!++++  REQUESTING  AND  CALCULATING  SYSTEM  PARAMETERS  ++4+++++++++++++++++++ 

1240 

! 

1250 

( 

1260 

PRINTER  IS  1 

1 

INPUTS  ARE  REQUESTED  ON  THE  CRT. 

1270 

! 

1280 

ChoiceD0 

| 

DETERMINING  C 

1290 

WHILE  Choi ce<>  1  AND  ChoiceOZ 

1300 

PRINT 

1310 

PRINT  “DO  YOU  DESIRE  TO  INPUT" 

1320 

PRINT  *  ( 1  )  THE  SPEED  OF  “  i 

1330 

PRINT  "S0UND,“ 

1340 

PRINT  “  OR  (2)  THE  AMBIENT  AIR  " 

5 

1350 

PRINT  "TEMPERATURE?" 

1360 

PRINT 

1370 

INPUT  Choice 

1380 

END  WHILE 

1390 

1 

1400 

IF  Choice= 1  THEN 

1 

DIRECTLY  ENTER  C, 

1410 

PRINT  "PLEASE  ENTER  C  IN  M/S." 

1420 

INPUT  C 

1  430 

END  IF 

1  440 

) 

1450 

IF  Choice=2  THEN 

I 

OR  CALCULATE  C  FROM  TEMP. 

1460 

PRINT  "PLEASE  ENTER  TEMP  IN  "i 

1470 

PRINT  "DEGREES  CELSIUS." 

1  480 

INPUT  Temp 

1490 

Tenp,=  Tenp  +  273.  15 

1500 

C=20.05*SQR<  Tenp  ) 

1510 

END  IF 

1520 

1 

*«•  DETERMINING  FREQ  ••• 

1530 

PRINT 

{ 

DIRECTLY  ENTER  FREQ. 

1540 

PRINT  "PLEASE  ENTER  THE  ECHOSOUNDER 

i 

1550 

PRINT  "  FREQUENCY  IN  HZ." 

1560 

PRINT 

1570 

INPUT  Freq 

75 


1580 

1590 

Choice«0 

! 

1 

««•  DETERMINING  KA 

1600 

WHILE  Choi ce<>  1  ANO  Cholce<>2  AND 

Cho 

ice<>3 

1610 

PRINT 

1620 

PRINT  "DO  YOU  DESIRE  TO  INPUT  “ 

1630 

PRINT  ”  (  1  )THE  SPHERICAL  “ 

i 

1640 

PRINT  "8ESSEL  ARGUMENT  K»A,“ 

1650 

PRINT  "  (2)  THE  SPHERE’S  “ 

1 

1660 

PRINT  "CIRCUMFERENCE," 

1670 

PRINT  “  OR  (3)  THE  SPHERE’S  “ 

1 

1680 

PRINT  "DIAMETER?" 

1690 

PRINT 

1700 

INPUT  Choice 

1710 

END  WHILE 

1720 

1 

1730 

IF  Choice=l  THEN 

! 

DIRECTLY  ENTER  KA, 

1740 

PRINT  “PLEASE  ENTER  K*A.“ 

1750 

INPUT  Ka 

1760 

END  IF 

1770 

1 

1780 

IF  Choice=2  THEN 

1 

CALCULATE  KA  FROM  THE 

1790 

PRINT  "PLEASE  ENTER  THE  TARGET 

i  1 

TARGET  CIRCUMFERENCE, 

1800 

PRINT  "CIRCUMFERENCE  IN  CM." 

1810 

INPUT  Circ 

1820 

Ka=Circ* ( 1 0" ( -2 ) )«Freq/C 

1830 

END  IF 

1840 

1 

I860 

IF  Choice=3  THEN 

I 

OR  CALCULATE  KA  FROM  THE 

1860 

PRINT  "PLEASE  ENTER  THE  TARGET 

“  i  ! 

TARGET  DIAMETER. 

1870 

PRINT  "DIAMETER  IN  CM." 

1880 

INPUT  Dian 

1890 

Ka-Diam«< 10" (-2 ) )*PI«Freq/C 

1900 

END  IF 

1910 

! 

1920 

! 

•»«  DETERMINING  REMAINING 

1930 

! 

SYSTEM  PARAMETERS  *•» 

1940 

K = 2  *  P I «Freq/C 

i 

CALCULATE  THE  WAVE  NUMBER. 

1950 

A=Ka/K 

1 

CALCULATE  THE  SPHERE  RADIUS. 

I960 

Lnax= INT ( Ka  +  2  1  ) 

l 

MAX.  SUMMATION  ORDER  CALCULATED. 

1  970 

1 

1980 

l 

1990 

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

2000 

I 

2010 

!++++  MODULE  CONCLUSION  ++++++++++++++++*+++++++++++++++++++++++++++++++ 

2020 

I 

2030 

SUBENO 

2040 

! 

2050 

!+  +  -t-  +  +  +  +  +  +  +  ++  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  -f-f  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  4  +  4  +  ++  +  +  +  +  +  +  +  +  + 

2060 

( 

2070 

! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$5$ 

2080 

1 
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2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320 

2330 

2340 

2350 

2360 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2580 

2590 

2600 

2610 


! 

I $$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB  Di  f  sctr<  Lmax  ,Ka  ,0j  C  )  ,  JC  )  ,QyC  )  ,  YC  )  .PC  )  ,S(*iC  )  ,DsC  )  > 


REM 

REM 

REM 

REM 

REM 

! 

(  +  +  +  + 


+++++++++++++++++++++ ++++++++++++++++++++4-+++++++++++++-M-+++++++++ 
THIS  MOOULE  CALCULATES  THE  ARRAYS  SM( * >  AND  DS<*>  REPRESENTING  THE 
SCATTERING  MOOULUS  AND  THE  DIFFERENTIAL  SCATTERING  CROSS  SECTION 
NORMALIZED  BY  THE  CROSS  SECTIONAL  AREA  OF  THE  SPHERE,  RESPECTIVELY . 
+++++  +  +  +  +  +  +  +++++  +  +  +++++++  +  +++++  +  +  +++  +  +  +  +  +  ++++  +  +++++  +  +++  +  +  +  +  +  ■*-+++++ 


MODULE  VARIABLE  DECLARATIONS 


REAL  Fac 
REAL  Rsm 

REAL  Ian 

REAL  Fact 
REAL  Phi 


AND 

I 


DEF  INI  T I ONS  ++++4-+++++  +++  ++++++++ 


FACTOR  COMMON  TO  RSM  AND  ISM. 

THE  REAL  COMPONENT  OF  THE 
SCATTERING  MODULUS.  SM. 

THE  IMAGINARY  COMPONENT  OF  THE 
SM. 

FACTOR  COMMON  TO  SM(PHI)  &  DS(PHI). 

THE  (POLAR)  ANGLE  OFF  THE  AXIS 
OF  PROPAGATION  OF  THE  INCIDENT 
PLANE  WAVE,  UITH  ORIGIN  AT  THE 
SPHERE  CENTER.  COS(PHI)  IS  THE 
ARGUMENT  OF  THE  LEGENDRE 
POLYNOMIALS  PC)  OF  ORDER  0 
THROUGH  LMAX . 


h ++++++++++++++++++++++++++++++++++ 


!++++  CALCULATING  THE  SM  AND  DSCS  ARRAYS  +++++++++++++++++++++++++++++++ 


CALL  OrablLnax  ,Ka  ,Dj(  *  )  ,J(  *  )  ) 

CALL  Oisb(Lmax  ,Ka,QyC),YC>) 

FOR  Ph l =0  TO  360 

CALL  Leg( Lmax  ,Ka  ,Phi  ,P<  *  )  ) 


CALCULATE  THE  DERIVATIVES  OF  THE 
REGULAR  SPHERICAL  BESSEL  FUNCTIONS 
OF  ORDER  0  THROUGH  LMAX. 

CALCULATE  THE  DERIVATIVES  OF  THE 
IRREGULAR  SPHERICAL  BESSEL  FUNC’NS 
I  OF  ORDER  0  THROUGH  LMAX. 


! 


! 


CALCULATE  THE  ARRAY  PC  >. 


Rsn=0 
I  sm=*0 
FOR  L=0  TO  Lnax 

Fac  =  ( Oj ( L  )•< 2«L+1  )‘P(L  )  )/<Oj<L  )*Dj(L  >  +  Dy(L )*Dy(L ) ) 
Rsn=Rsra+(Fac*Oy(L  )  ) 

Isi»i=Isf';+(Fac«Dj(L  )  ) 

NEXT  L 


!  CALCULATING  THE  REAL  AND  IMAGINARY 
!  COMPONENTS  OF  THE  SM. 


Fact=R5m*RsM+Isn«Isfn 


2B20 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 


Sm( Phi  )=SQR( Fac t  ) 

Os<  Phi  )“4« ( Fact  )/(Ka*Ka) 

NEXT  Phi 


CALCULATING  THE  ARRAY  REPRESENTING 
THE  SM( PHI  ). 

CALCULATING  THE  ARRAY  REPRESENTING 
THE  DIFFERENTIAL  SCATTERING  CROSS 
SECTION  NORMALIZED  BY  THE  CROSS 
SECTIONAL  AREA  OF  THE  SPHERE. 


! +++++++++++++i 


I++++  MODULE  CONCLUSION  +++h 


H+++++++++ 


SUBEND 


!+++++++++++++++++++++++++++++++++++++++++++++ ++++++4 


!$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

! 
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2840 
2850 
2860 
2870 
2880 
2890 
2900 
2910 
2920 
2930 
2940 
2950 
2980 
2970 
2980 
29S0 
3000 
3010 
3020 
3030 
3040 
3050 
3060 
3070 
3080 
3090 
3100 
31  10 
3120 
3130 
3140 
3150 


I 

!$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB  Orsb( Lmax  ,Ka  ,0 j ( »  )  ,  J( » >  ) 

REM 
REM 
REM 
REM 
REM 
REM 
! 

!++++  CALCULATING  THE  DERIUATIUES  OF  THE  RSB  +++++++++++++++++++++++++++ 


+++++++++++++-M-++++++++++++++4-+++++++++++++++++++ +++++++++++++++++ 

THIS  MOOULE  CALCULATES  THE  DERIUATIUES  OF  THE  REGULAR  SPHERICAL 
BESSEL  FUNCTIONS  ( DRSB )  OF  ARGUMENT  KA  AND  ORDERS  0  THROUGH  LMAX . 
THE  CALCULATION  OF  THE  DER IUATIUE  INUOLUES  THE  RSB'S  OF  THE  SAME 
AND  PREUIOUS  ORDER. 

+++++++++++++++++++++ H 


CALL  Rsh< Lnax  ,Ka  , J< *  ) ) 

RAO 

Dj ( 0  )"( - J ( 0  >+COS<  Ka  )  )/Ka 
FOR  L“ 1  TO  Lwax 

D j ( L  )«-( ( L+ 1  )/Ka  >« J ( L  >  + J ( L- 1  ) 
NEXT  L 


I 


! 


CALCULATE  THE  RSB  UALUES  J<«). 


!  ANGLE  UNITS  FOR  DRSB  CALC’NS. 

!  CALCULATE  THE  INITIAL  DRSB  DJ<0>. 
!  CALCULATE  THE  REMAINING  DRSB. 


I  +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


I++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 
SUBEND 

! 

! 

I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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3160 

3170 

3160 

3190 

3200 

3210 

3220 

3230 

3240 

3250 

3260 

3270 

3280 

3290 

3300 

3310 

3320 

3330 

3340 

3350 

3360 

3370 

3380 

3390 

3400 

3410 

3420 

3430 

3440 

3450 

3460 

3470 

3480 

3490 

3500 

3510 

3520 

3530 

3540 

3550 

3560 

3570 

3580 

3590 

3600 

3610 

3620 

3630 

3640 

3650 

3660 

3670 

3680 

3690 

3700 


I 

! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$5$$$$$$$$$$$$$$$$ 
SUB  Rsb( Lmax ,Ka ,J<  «  ) ) 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM  THIS  MODULE  CALCULATES  THE  RESULAR  SPHERICAL  BESSEL  FUNCTIONS 
REM  ( RSB  )  OF  ARGUMENT  KA  AND  ORDERS  0  THROUGH  LMAX.  THE 
REM  CALCULATION  USES  THE  CONTINUED  FRACTION  APPROACH. 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++♦+ 
I 

I++++  MODULE  UARIABLE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 


REAL  R1 
REAL  Numr 

REAL  Denr 

REAL  Scexp 

REAL  Nf lag 
REAL  Dflag 
REAL  Bnm 

REAL  Bine 


I 

I 

!  THE  EUOLUING  RATIO  OF  J(0)/J(L>. 

1 

!  THE  EUOLUING  NUMERATOR  FACTOR  USED 
I  IN  CALCULATION  OF  RL. 
t 

!  THE  EUOLUING  DENOMINATOR  FACTOR 
I  USEO  IN  CALCULATION  OF  RL. 

I 

!  EUOLUING  SCALING  EXPONENT  FOR  UERY 
I  LARGE  OR  UERY  SMALL  UALUES  OF  RL. 

! 

!  0  FLAG  FOR  NUMR. 

1 

i  0  FLAG  FOR  DENR. 

( 

I  EUOLUING  TERM  USED  IN  THE 
t  CALCULATION  OF  NUMR  AND  DENR. 

I 

!  INCREMENT  USED  IN  THE  EUOLUTION 
!  OF  BNM. 

i 

1 


+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


++++  CALCULATING  THE  RSB  ++++++++++++++++++++++++++++++++++++++++++++++ 


RAD 


ANGLE  UNITS  FOR  CALC’NS. 


J(  0 )«<  SINIKa  )  )/Ka 


CALCULATE  THE  INITIAL  ••• 
•  *»  RSB ,  J<0>.  ••• 


FOR  L-l  TO  Lmax 
Rl-I  .0 
Numr-0 . 

Denr«0 . 
Scexp«0. 

Nf lag«0 . 

Df lag“0. 

8nm« I . 0/Ka 
Blnc“2 .0/Ka 


CALCULATING  THE  RATIO  ••• 
•  J<  0  >/J(L  ) . 

INITIALIZING  RL. 

INITIALIZING  NUMR. 
INITIALIZING  DENR. 
INITIALIZING  SCEXP. 
INITIALIZING  NFL AG. 
INITIALIZING  DFLAG. 
INITIALIZING  BNM. 

CALCULATE  BINC. 
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3710 

1  *  CALCULATE  THE  L  “NAKED” 

3720 

FOR  Ordcnt=1  TO  L 

!  «  NUMR  FACTORS  OF  J(0)/J(L>.  • 

3730 

Bnm=8nm+B l nc 

!  INCREMENT  BNM . 

3740 

IF  Nf lag=0  THEN 

!  CHECK  NUMR-0  FLAG  NOT  SET  j 

3750 

Nui»ir“Bnm-Numr 

1  UPDATE  NUMR  IF  NUMRO0. 

3760 

IF  NunrOQ  THEN 

3770 

Rl=Rl*Nunr 

!  UPDATE  RL  IF  UPDATED  NUMRO0. 

3760 

Nunr  = 1 . 0/Nunr 

!  PREPARE  NUMR  FOR  NEXT  EVOLUTION. 

3790 

!  CHECK  SCALIN6  OF  UPDATED  RL. 

3800 

IF  A8S(R1  )> 1 0A250 

THEN 

CALL 

Scale( R1  ,Scexp  ) 

3810 

ELSE 

3820 

Nf  lag= 1 

3830 

END  IF 

3840 

ELSE 

!  IF  NFLA6O0  THE  PREVIOUS  STEPS  ARE 

3850 

Nf lag=0 

!  DEFERRED  FOR  ONE  INCREMENT  OF  BNM. 

3860 

ENO  IF 

3870 

NEXT  Ordcnt 

3880 

1  •  CALCULATE  THE  REMAINING  NUMR  * 

3890 

!  »  AND  DENR  TERMS  OF  J(0)/J(L).  • 

3900 

!  NOTE  WHEN  NUMR-DENR  THEN 

3910 

1  RL  =  J(0  )/J(L  ). 

3920 

WHILE  NunrODenr  OR  Nflag=l  1 

OR  Of Iag= 1 

3930 

Bnm“Bnm+Binc 

1  INCREMENT  BNM. 

3940 

) 

3950 

IF  Nf lag=0  THEN 

1  CHECK  NUMR=0  FLAG  NOT  SET; 

3960 

Nunr=Bnm-NuMr 

!  UPDATE  NUMR  IF  NUMRO0. 

3970 

IF  Numr<>0  THEN 

3980 

R1*R1 »Nunr 

1  UPDATE  RL  IF  UPDATED  NUMRO0 . 

3990 

Nunr= 1 . 0/Numr 

1  PREPARE  NUMR  FOR  NEXT  EVOLUTION. 

4000 

!  CHECK  SCALING  OF  UPDATED  RL . 

4010 

IF  ABS<R1  )> 1 0A2S0 

THEN 

CALL 

Scale! R1  ,Scexp  > 

4020 

ELSE 

4030 

Nf  lag=  1 

4040 

ENO  IF 

4050 

ELSE 

!  IF  NFLAGO0  THE  PREVIOUS  STEPS  ARE 

4060 

Nf 1 ag=0 

!  DEFERRED  FOR  ONE  INCREMENT  OF  BNM. 

4070 

ENO  IF 

4080 

1 

4090 

IF  Of lag=0  THEN 

!  CHECK  DENR-0  FLAG  NOT  SET; 

4100 

Denr=Bnn-Denr 

!  UPDATE  DENR  IF  DENRO0. 

4110 

IF  Denr < >0  THEN 

4120 

Rl-RI/Oenr 

!  UPDATE  RL  IF  UPDATED  DENRO0. 

4130 

Denr*1 . 0/Oenr 

I  PREPARE  DENR  FOR  NEXT  EVOLUTION. 

4140 

1  CHECK  SCALING  OF  UPDATED  RL . 

4  150 

IF  ABSIRl  >>!CA2SC 

THEN 

CALL 

Sea  1  e ( R 1  , Scexp  ) 

4160 

ELSE 

4170 

Of lag=  1 

4  180 

ENO  IF 

4  1  90 

ELSE 

1  IF  DFLAGO0  THE  PREVIOUS  STEPS  ARE 

4200 

0  f  1  ag  =  0 

!  DEFERRED  FOR  ONE  INCREMENT  OF  BNM. 

4210 

END  IF 

4220 

! 

4230 

END  WHILE 

4  2  4C 

I 
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4250  I 

4260  J<L)-J<0>/(Rl»10"<Scexp>)  !  J ( L )  DETERMINED .  •«* 

4270  ! 

4280  NEXT  L 

4290  ! 

4300  ! 

4310  !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

4320  ! 

4330  !++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

4340  ! 

4350  SUBEND 

4360  ! 

4370  !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

4380  ! 

4390  !$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$«$$$ 

4400  ! 
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4410  ! 

4420  !$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$*$$$$$$$$$$$$$$$$$$$$$$$$$ 

4430  SUB  Scale! R1  ,Scexp  ) 

4440  REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
4450  REM  THIS  MODULE  SCALES  RL  TO  MAINTAIN  NUMERICAL  ACCURACY  IN  THE 
4460  REM  ONGOING  CALCULATION  OF  RL. 

4470  REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
4480  I 

4490  !+  +  <■+  SCALING  RL  +++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

4500  ! 

4510  ! 

4520  IF  R1 > 1 0*250  THEN 
4530  RI-R1 *  1 0“ ( -2S0  > 

4540  Scexp=Scexp+2S0 

4550  END  IF 

4560  IF  R1 < 1 0A  < -250  )  THEN 
4570  R1«R1*10A250 

4580  Scexp“Scexp-2S0 

4590  END  IF 

4600  ! 

4610  ! 

4620  ! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

4630  ! 

4640  !++++  MOOULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

4650  ! 

4660  SUBEND 

4670  I 

4 £80  ) +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

4690  ! 

4710  ! 
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4720 

4730 

4740 

4750 

4760 

4770 

4780 

4790 

4800 

4810 

4820 

4830 

4840 

4850 

4860 

4870 

4880 

4890 

4900 

4910 

4920 

4930 

4940 

4950 

4960 

4970 

4980 

4990 

5000 

5010 

5020 

5030 


! 

! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$5$ 

SUB  Disb! Lmax  ,Ka  ,Dy( * >  ,Y< *  )  ) 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM  THIS  MODULE  CALCULATES  THE  DERIVATIVES  OF  THE  IRREGULAR  SPHERICAL 
REM  BESSEL  FUNCTIONS  OF  ARGUMENT  KA  AND  ORDERS  0  THROUGH  LMAX.  THE 
REM  CALCULATION  OF  THE  DERIVATIVE  USES  THE  ISB'S  OF  THE  SAME  AND 
REM  PREVIOUS  ORDER. 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
! 

! ++++  CALCULATING  THE  DERIVATIVES  OF  THE  ISB  +++++++++++++++++++++++++++ 


CALL  IsblLmax  ,Ka  ,Y< *  )  ) 


! 

I  CALCULATE  THE  ISB  VALUES. 


RAD  !  ANGLE  UNITS  FOR  CALC’NS. 

Dy(0)»<-Y(0>+SIN(Ka))/Ka  I  THE  INITIAL  ISB  DERIVATIVE. 

FOR  L*=  I  TO  Lmax  I  THE  REMAINING  ISB  DERIVATIVES. 

Dy(  L  )«*-(  (  L+ 1  )/Ka  ) * Y<  L  )+Y(L-1  ) 

NEXT  L 

! 

I 

!+  +  +  +  +  +  +  +  +  +  ++++++  +  ++++++  +++  ++++++++++  +  +  -*■  +++  +  +++++++  +  +  +  +++  +  +++  +  +  +  +  +  +  +++  +  + 


!+  +  +  +  MODULE  CONCLUSION  +++++++++■+++++++++++++++++++++++++++++++++++++++ 

I 

SUBEND 

I 

+++++++++++++++++++++++++++++++T+++++++++++++++++++++++++++++++++++++++ 


!  $$$$$$£$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
I 
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5040 

5050 

5060 

5070 

5080 

5090 

5100 

5110 

5120 

5130 

5140 

5150 

5160 

5170 

5180 

5190 

5200 

5210 

5220 

5230 

5240 

5250 

5260 

5270 

5280 

5290 

5300 

5310 

5320 

5330 


! 

! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB  Isb(  Lna.x  .Ka  ,Y(»  ) ) 

REM  ++++++ ++  +  -M-++-M-+++++ ++++++++++ ++++++++++++++++++++++++++++++++++++ 

REM  THIS  MODULE  CALCULATES  THE  IRREGULAR  SPHERICAL  BESSEL  FUNCTIONS 
REM  OF  ARGUMENT  KA  AND  ORDER  L  THROUGH  LMAX .  THE  CALCULATION  IS 
REM  ACCOMPLISHED  BY  FORWARD  RECURSION. 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
! 

!++++  CALCULATING  THE  ISB  ++++++++++++++++++++++++++++++++++++++++++++++ 

I 


RAD 

Y ( 0  )=-CuS( Ka  )/Ka 

Y<1  )“<  Y( 0  >-SIN( Ka  )  )/Ka 

FOR  L-2  TO  Lmax 

Y(L)»((2*L-1  )/Ka  )«Y( L- 1  )- Y<  L-2 ) 
NEXT  L 


!  ANGLE  UNITS  FOR  CALC’NS. 

!  INITIAL  ISB’S  CALCULATED. 

I  REMAINING  ISB’S  CALC’D. 


++++++++++++++++++++++++++ H 


!•  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  + 


++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

SUBEND 

•  I 


!  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

! 
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5340 

5350 

5360 

5370 

5380 

5390 

5400 

5410 

5420 

5430 

5440 

5450 

5460 

S470 

5480 

S490 

5500 

5510 

5520 

5530 

5540 

5550 

S560 

5570 

5580 

5590 

5600 

5610 

5620 

5630 

5640 

5650 


I 

!$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB  Leg< Lnax  ,Ka  ,Phi  ,P( «  )  ) 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM  THIS  MODULE  CALCULATES  THE  LEGENDRE  POLYNOMIALS  OF  ARGUMENT 
REM  COS(PHI)  AND  ORDERS  0  THROUGH  LMAX .  CALCULATION  IS  ACCOMPLISHED 
REM  BY  FORUARO  RECURSION. 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

!++++  CALCULATING  THE  LEGENDRE  POLYNOMIALS  +++++++++++++++++++++++++++++ 


OEG 

X“COS( Phi  ) 


1  ANGLE  UNITS  FOR  CALC’NS. 

I  ARGUMENT  OF  THE  LEG.  POL’S. 


P(0)-l 
P<  1  >»X 

FOR  L=2  TO  Lnax 

P(L)=(2-< 1 /L  )  )»X*P( L- 
NEXT  L 


INITIAL  LEG.  POL’S.  CALC’D. 


!  REMAINING  LEG. 
1  )-( 1 -( 1/L))»P(L-2) 


POL’S.  CALC ’ 0 . 


! 

1 

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

l 

SUBEND 

1 

! ++++++++++++++++++++++++++++++++++++++++ +++++++++++++++++++++++++++++++ 
l 


I  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$ 

| 
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5660 
5670 
5680 
5G90 
5700 
5710 
5720 
5730 
5740 
5750 
5760 
5770 
5780 
5790 
5800 
58)0 
5820 
5830 
5840 
5850 
5860 
5870 
5880 
5890 
5900 
5910 
5920 
5930 
5940 
5950 
5960 
5970 
5980 
5990 
6000 
6010 
6020 
6030 
6040 
6050 
6060 
6070 
6080 
6090 
6)00 
6)10 
6)20 
6)30 
6  I  40 
6150 
6)60 
6  I  70 
6180 
6190 


I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB  Outpt  <  C  ,Freq  ,A ,Ka ,Sn( « ),0a(»)> 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM  THIS  MODULE  CONTROLS  THE  OUTPUT  OF  THE  CALCULATED  DSCS  AND  SM  IN 
REM  GRAPHICAL  PRESENTATIONS. 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

!++++  MODULE  UARIABLE  DECLARATIONS  AND  DEFINITIONS  +  +  i 


REAL  Desire 
REAL  Rpt 


INTERACTIVE  LOOP  CONTROL. 

POLARPLT  SCALING  CONTROL i  A 
REPETITION  INDEX. 


++++++++++++++++++++++++++++++++4++++++++++H 

++++  DETERMINING  THE  DESIRED  OUTPUT  +++++++++++++++++++++++++++++++++++ 


AND  0esire<>3 


Oes ire-0 

WHILE  DesireOI  AND  0esire<>2 
PRINT  "DO  YOU  DESIRE" 

PRINT  “  (I)  A  PLOT  OF  THE  DSCS 

PRINT  “  (2)  A  PLOT  OF  THE  SM . “ 

PRINT  "  OR  (3)  BOTH?" 

INPUT  Desire 
END  WHILE 
PRINTER  IS  701 

WHILE  Desire<4 

IF  DesireOI  OR  Oesire-3  THEN 
FOR  Rpt- I  TO  2 

CALL  Polarpl t (Ds( »  )  ,Rpt  ) 

PRINT 
NEXT  Rpt 
END  IF 

IF  Desire=2  THEN 
CALL  Recp 1 t  (  Sn  ( *  )  ) 

PRINT 
ENO  IF 

PRINT  “  " ,"C  =" iCi "M/S. “ 

PRINT  "  " , “F  =" i Freq i “ HZ .  ” 

PRINT  ”  "  ,"A  ="  A s  "M.  " 

PRINT  "  " , " X • A  ="  i  K  a  i  "  .  " 

PRINT  "  “/'NORMALIZED  "i 
PRINT  "BCKSCTR  =  " iOs<  1 80  ) 

PRINT  CHR$(  1 2  ) 


I 

!  PLOTTING  THE  DSCS. 

! 

!  PLOTTING  THE  SM. 

! 

!  ECHO  SIGNIFICANT  CALCULATED  UALUES 
!  AND  INPUT  PARAMETERS. 


!  PRINTER  FORM  FEED. 

I 

IF  Oesire=l  OR  0esire=2  THEN  Oesire-4 
IF  Desire=3  THEN  Oesire=2 

END  WHILE 
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++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


6200 
6210 
6220 

6230  !++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

6240  I 

6250  SUGENO 

6260  ! 

6270  ! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

6280  ! 

6290  !*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$** 

6300  ! 
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6310 

6320 

6330 

6340 

6350 

6360 

6370 

6380 

6390 

6400 

6410 

6420 

6430 

6440 

6450 

6460 

6470 

6480 

6490 

6500 

6510 

6520 

6530 

6540 

6550 

6560 

6570 

6580 

6590 

6600 

6610 

6620 

6630 

6640 

6650 

6660 

6670 

6680 

6690 

6700 

6710 

6720 

6730 

6740 


! 

I $$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB  Polarp 1 1 < Ds( «  )  ,Rpt  ) 

REM  ++++++++++++++++++++++++++ +++++++++++++++++++++++++++++++f++++++++ 

REM  THIS  MOOUIE  MAPS  THE  CALCULATED,  NORMALIZED  DSCS  ARRAY  DS( *  >  ON 

REM  LINEAR  POLAR  PLOTS.  THE  FIRST  PLOT  INCLUDES  THE  ENTIRE  DSCS, 

REM  WHILE  THE  SECOND  PLOT  SHOWS  THE  SCATTERING  DETAIL  AROUND  THE 
REM  NORMALIZED  RANGE.  PLOT  SELECTION  IS  GOUERNED  BY  A  REPETITION 
REM  INDEX,  Rpt,  AND  IS  CONTROLLED  IN  THE  CALLING  MODULE  (OUTPT). 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

!++++  MODULE  UARI ABLE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 


REAL 

Xguroax 

•  THE  MAXIMUM  ABSCISSA  UALUE  IN 
!  GRAPHICS  DISPLAY  UNITS. 

REAL 

Ygunax 

! 

!  THE  MAXIMUM  ORDINATE  UALUE  IN 
!  GRAPHICS  DISPLAY  UNITS. 

| 

REAL 

Gumax 

!  THE  MINIMUM  OF  XGUMAX  AND  YGUMAX  i 
!  USED  TO  DETERMINE  PLOTTING  AREA. 

| 

REAL 

Dtnx 

1  THE  MAX.  UALUE  OF  THE  ARRAY  OS<»>. 

| 

REAL 

Endpt 

!  THE  MAX.  RADIAL  SCALE  UALUE  IN 

!  USER  DEFINED  UNITS. 

1 

REAL 

Ring 

! 

!  GRID  RING  SPACING  IN  USER  UNITS. 

REAL 

Tick 

!  RADIAL  AXIS  TICK  MARK  SPACINGS  IN 
!  USER  DEFINED  UNITS. 

! 

1 

! 4 444 +444-4 4444444+444+4444 4+44444+44 4444 4 4+4 444+44+4444+4+4+4+4+444444++ 

i 

++++  INITIALIZING  THE  PLOTTER  ++++++++++++++t++++++++++++++++++++++++++ 


GINIT 

PLOTTER  IS  CRT  ."INTERNAL" 

GRAPHICS  ON 
GCLEAR 

! 

! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
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I++++  GRAPH  CONSTRUCTION  AND  LABELLING  +++++++++++++++++++++++++++++++++ 


Xgumax“100*MAX( I  .RATIO  ) 
Ygumax” 1 00»MAX( I  ,  I /RAT  10  ) 
6umax“MIN<  Xgumax  ,  Ygumax  ) 


6750  I 
6760 
6770 
6780 
6790 
6800 
6810 
6820 
6830  CSIZE  3 
6840  LINE  TYPE  I 
6850  MOUE  Gumax  ,Gumax/2 
6860  LORG  8 
6870  LABEL  “0“ 

6880  MOUE  Gumax/2  .Gumax 

6890  LORG  6 

6900  LABEL  “PI/2" 

6910  MOUE  0 ,6umax/2 
6920  LORG  2 
6930  LABEL  “PI" 

6940  MOUE  Gumax/2  .0 

6950  LORG  4 

6960  LABEL  “3»PI/2" 

6970 

6980  UIEWPORT  5  .Gumax-5  .5  ,Gumax~5 

6990 

7000 

7010 

7020  Dtmax“0 

7030  IF  Rpt-1  THEN 

7040  FOR  1=0  TO  360 

7050  IF  Dtmx<Ds< I  )  THEN  Otmx“Ds< I  ) 

7060  NEXT  I 

7070  END  IF 

7080  Endpt  =  INT(Dtmx+l  ) 

7090  IF  Endpt  <2  THEN  Endpt-2 
7100  SHOW  -Endpt  .Endpt  , -Endpt  .Endpt 
71  10 

7120  LINE  TYPE  4 

7130  Ring=INT<Endpt/5+l ) 

7140  Tick=.2*Ring 

7150  AXES  Tick. Tick. 0.0. 5. 5. 2 

7160 

7170  LINE  TYPE  3 
7180  OEG 
7190  LORG  7 

7200  FOR  R=1  TO  Endpt/Ring 

7210  LINE  TYPE  1 

7220 
7230 
7240 
7250 
7260 
7270 
7280 
7290 


!  ••  LABELLING  MAJOR  RADIALS  •• 


•»  DEFINES  THE  TOTAL  PL0TTIN6  ** 
»»  AREA  IN  6RAPH  UNITS.  ** 

•«  SCALING  THE  GRAPH  TO  *« 

**  USER  UNITS.  •* 

DETERMINING  MAX.  RADIUS  UALUE 
REQ'D  FOR  THE  PLOT. 


!  ISOTROPICALLY  SCALES  THE  GRAPH. 

I 

!  *•  MAJOR  RADIALS  CONSTRUCTION  •* 


I 

I  ••  RANGE  RING  CONSTRUCTION  •• 


MOUE  0  ,R*R i ng 

LABEL  R  *R  j  ng  !  LABELLING  THE  RANGE. 

LINE  TYPE  3 

FOR  Angl e=0  TO  360  !  DRAWING  THE  RING. 

PLOT  R*Ring*COS(Angle),R*Ring*SIN(Angle) 

NEXT  Angle 
PENUP 
NEXT  R 
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7300 

73)0 

7320 

7330 

7340 

73S0 

7360 

7370 

7380 

7390 

7400 

7410 

7420 

7430 

7440 

7450 

7460 

7470 

7480 

7490 

7500 

7510 

7520 

7530 

7540 

7550 

7560 

7570 

7580 

7590 

7600 

76)0 

7620 

7630 

764C 

7650 

7660 


FOR  Angle-30  TO  150  STEP  30  I  30  DE6REE  RADIALS  •* 

IF  Angle<>90  THEN  I  *»  CONSTRUCTION.  •» 

PLOT  Endp  t  «COS<  Angle  )  ,Endpt  »S IN( Angl e  )  .  —  2 
DRAW  Endpt  »C0S< 180+Angle  )  .Endpt«SIN( 180+Angle  ) 

ENO  IF 
NEXT  Angle 
PENUP 


I 

14444  PLOTTING  THE  CALCULATED  DSCS  4444444 


4444444444 


LINE  TYPE  1 
DEG 

FOR  1-0  TO  360  I  PLOTTING  RECTAN6ULAR  COORDINATES 

PLOT  Os( I  )«COS< I  )  ,Ds< I  )«SIN( I )  !  IN  POLAR  FORMAT. 

NEXT  I 


44444444444444444+ 


+4444444444+444444444444444444 


4444  PRINTING  THE  PLOT  44444444444444444444444444+444444444444444444444 


DUMP  GRAPHICS  #701 


I 

! 


! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++4++++++++ 
I 

14444  MODULE  CONCLUSION  4+4444444444444444444444444444444444444.444444444 


GRAPHICS  OFF 
SUBEND 

I 


I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

! 
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SUB  Recp 1 1 ( Sn < •  )  ) 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM  THIS  MODULE  MAPS  THE  CALCULATED  SCATTERING  MODULUS  ARRAY  SM(»>  ON 
REM  A  SEMI -LOG  PLOT. 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
! 

!++++  MODULE  UARIABLE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 


REAL  Xgumax 


!  THE  MAXIMUM  ABSCISSA  UALUE  IN 
I  GRAPHICS  DISPLAY  UNITS. 


REAL  Ygunax 


REAL  Dtmx 


REAL  Dtnn 


REAL  Mxy 


!  THE  MAXIMUM  ORDINATE  UALUE  IN 
!  GRAPHICS  DISPLAY  UNITS. 

1 

!  THE  MAX.  UALUE  OF  THE  ARRAY  SM<«>. 

! 

I  THE  MIN.  UALUE  OF  THE  ARRAY  SM<*>. 

I 

!  THE  MAX.  ORDINATE  SCALE  UALUE  IN 
!  USER  DEFINED  UNITS. 


REAL  Mny 


•  THE  MIN.  ORDINATE  SCALE  UALUE  IN 
!  USER  DEFINED  UNITS. 


^  +  +  +  +  +  4-  +  +  *M-  +  +  +  +  ++++++++++++++++++++ 


!++++  INITIALIZING  THE  PLOTTER  +H 


h++++++++++++++++++++ 


GINIT 

PLOTTER  IS  CRT  INTERNAL" 

GRAPHICS  ON 
GCLEAR 

I 

I +  +  +  +  +  +  +  +  +  +  +++++++  +  +  +  +++++  +  ++++++  ++++  +  +  +  +  +  +++  +  +  ++  +  +  +  +  ++++++++++++++  +++  +  + 
) 

!++++  GRAPH  CONSTRUCTION  AND  LABELLING  +++++++++++++++++++++++++++++++++ 

i 

! 

Xgumax=100*MAX( I  .RATIO  ) 

Ygumax=t00»MAX( 1  ,1/RATIO) 

1 

CSIZE  4 


LORG  4 

MOUE  Xgumax/2,0 

LABEL  "PHI  (DEGREES)’ 


I  «*  LABELLING  THE  ABSCISSA  AXIS  •  • 


I 


81  70 

8180  OEG 

8190  LOIR  90 

8200  LORG  S 

8210  MOVE  0  .Ygurcax/Z 

8220  LABEL  'SCTRG  MODULUS" 

8230  LOIR  0 
8240 

8250  VIEWPORT  1 S  ,Xgumax-5 , 1 0  .  Ygumax-5 

82G0 

8270 

8280  Dtmx»0 
8290  Dtein*Sm(  180  > 

8300  FOR  I»0  TO  180 
8310  IF  Dtfix<Sn<  I  )  THEN  Dinx»Sn<  I  > 
8320  IF  Dtmn>Sm< I  )  THEN  Otnn=Sn( I) 
8330  NEXT  I 

8340  Mxy=INT(  LGT ( Dtfix  )  )  +  I 
8350  Mny*INT ( LGT( Dt mn  )  ) 

83G0 

8370  WINDOW  0,1 80. Mny.Mxy 
8380 

8390  AXES  30,1  ,0,Mny  ,1  ,1,1.5 
8400 

8410  LINE  TYPE  3 
8420  FOR  I=Mny  TO  Mxy 
8430  FOR  J=1  TO  9 

8440  MOVE  0  ,I+LGT( J  ) 

8450  PLOT  180  ,I+L6T(  J  )  ..-1 

8460  PENUP 

8470  NEXT  J 

8480  NEXT  I 

8490 

8500  FOR  1*60  TO  180  STEP  60 

8510  MOVE  I  ,Mny 

8520  PLOT  I  .Mxy  ,- 1 

8530  PENUP 

8540  NEXT  I 

8550 

8560  CLIP  OFF 
8570  LINE  TYPE  1 
8580 

8590  LORG  6 

8600  FOR  I>0  TO  180  STEP  30 

8610  MOVE  I ,Mny 

8620  LABEL  I 

8630  NEXT  I 

8640 

8650  LORG  8 
8660  FOR  I-Mny  TO  Mxy 
8670  MOVE  0,1 

8680  LABEL  1C~I 

8690  MOVE  0 , 1 +LGT ( 5  ) 

8700  LABEL  < 10A<  1  +  1  )  )/2 

8710  NEXT  I 


!  •«  LABELLING  THE  ORDINATE  AXIS  » 


!  **  GRAPH  CONSTRUCTION  *» 

I  DEFINES  THE  TOTAL  PLOTTING  AREA 
!  IN  GRAPH  UNITS. 

!  ••  ORDINATE  USER  UNITS  *• 

I  »•  DETERMINED. 


!  ANISOTROPICALLY  SCALES  THE  GRAPH 
I  TO  USER  UNITS. 

!  PLOTS  THE  AXES. 

! 

I  *  PLOTS  THE  HORIZONTAL 
1  *  LOGARITHMIC  6RID  LINES  I.A.W. 

!  *  THE  VERTICAL  SCALE. 


t 

!  *  PLOTS  THE  VERTICAL  GRID  LINES 
I  *  I.A.W.  THE  HORIZONTAL  SCALE. 


I 


I 

I  «•  LABELLING  THE  ABSCISSA  SCALE 


! 

I  **  LABELLING  THE  ORDINATE  SCALE 


!  POWER  OF  10  LABELS. 

1  HALF  THE  NEXT  POWER  OF  10  LABELS 


8720 

8730 

8740 

8750 

8760 

8770 

8780 

8790 

8800 

8810 

8820 

8830 

8840 

8850 

8860 

8870 

8880 

8890 

8900 

8910 

8920 

8930 

8940 

8950 

8960 

8970 

8980 

8990 

9000 

9010 


CLIP  ON 

i 


++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++T 


++++  PLOTTING  THE  CALCULATED  SCATTERING  MODULUS  +++++++++++++++++++++++ 


LINE  TYPE  1 
FOR  1-0  TO  180 

PLOT  I  ,LGT ( Sm( 1  +  180') 
NEXT  I 


!  THIS  PLOT  IS  FOR  THE  PLANE  WAUE 
!  INCIDENT  FROM  THE  RIGHT  <0  DEG’S >. 

I  USE  SM( I  )  IN  PLACE  OF  SM< 1  +  180) 

1  FOR  A  PLANE  WAUE  INCIDENT  FROM  THE 
!  LEFT  (180  DEG'S ). 


! ++++++++++++++++++++++++ +H 


!+  +++  PRINTING  THE  PLOT  +++++++++++++++++++++++++++++++4-++++++++++++++++ 


DUMP  GRAPHICS  #701 


I 


!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

) 

I+  +  +  +  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

GRAPHICS  OFF 
SUBEND 


I 


! +++++++++++H 


M-  +  +  +  -I 


I  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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BACKSCATTERED  CROSS  SECTION  PROGRAM 


Examples  of  the  output  of  this  program  are  Figures  6  and  7,  located  in 
Chapter  III. 


10  !$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$*$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

20  ! 

30  REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

40  REM  ACOUSTICALLY  HARD  SPHERE  BACKSCATTERING  OF  PLANE  ACOUSTIC  WAUES 

50  REM 

60  REM  THIS  PROGRAM  CALCULATES  AND  GRAPHS  THE  NORMALIZED.  BACKSCATTERED 

70  REM  CROSS  SECTION  FOR  ACOUSTICALLY  HARD  SPHERES  OF  A  SELECTED  RANGE 

80  REM  OF  DIAMETERS  0  ENSONIFIED  BY  AN  INCIDENT  PLANE  WAUE. 

90  REM 

100  REM 

110  REM 

1 20  REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

130  ! 

140  I++++  MAIN  PROGRAM  UARI ABLE  DECLARATION  AND  DEFINITION  +++++++++++++++++ 

150  ! 

160 

170  REAL  Nbs( 301  > 

180 
190 
200 

210  REAL  K 
220 

230  REAL  A 
240 

250  REAL  Ka 
260 
270 
280 

290  REAL  Kanin 
300 

310  REAL  Kanax 
320 

330  REAL  Stp 
340 
350 

360  1 ++++++ +++++++++++++++++++++++++ ++++++++++++++++++++++++++++++++++++++++ 

370  I 


THE  ARRAY  REPRESENTING  THE 
NORMALIZED,  BACKSCATTERED  CROSS 
SECTION  ( NBCS  ) . 

THE  WAVE  NUMBER,  2«PI*FREQ/C. 

RADIUS  OF  THE  SPHERE. 

K*A ,  THE  ARGUMENT  OF  THE  REGULAR 
AND  IRREGULAR  SPHERICAL  BESSEL 
FUNCTIONS  AND  THEIR  OERIUATIUES. 

INTEGER  OF  THE  MIN.  DALUE  OF  KA. 

INTEGER  OF  THE  MAX.  UALUE  OF  KA. 

KA  INCREMENT. 
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380  !++++  MAIN  PROGRAM  +++++++++++++++++++++++++++++++++++++++++++++++++++++ 

390  1 

400  CALL  Bcksc tr< Nbs( * >  .Kami n  ,Kamax  ,S tp  )  I  CALCULATE  THE  NBS(»>. 

410  * 

420  CALL  Recplt(Nbs( *  )  .Kanin  ,Kamax  ,Stp)  !  OUTPUT  THE  NBS<«>. 

430  ! 

440  i +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
450  ! 

460  !++++  MAIN  PROGRAM  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++ 

470  ! 

480  PRINTER  IS  70! 

490  PRINT  CHRK 12  ) 

500  END 

510  ! 

530  ! 

540  ! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$*$$$$$$$$ 

550  ! 
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560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 
770 
780 
790 
800 
810 
820 
830 
840 
850 
860 
870 
880 
890 
900 
910 
920 
930 
940 
950 
960 
970 
980 
990 
1000 
1010 
1020 
1030 
1040 
1  050 
1060 
1070 
1080 
1090 


I 

! $$$$$$$$$$$$$$$$$$$$$$S$$$$$$$$$$$$$$$$$$$$s$$$$$$$$$$$$$$S$$$$$*$S$S$S 

SUB  Bcksctr(Nba<  *  )  .Kanin  .Kamax  ,Stp  > 

REM  ++++++++++++++++++++++++++++++++++++++++ +++++++++++++++++++t++++++ 
REM  THIS  MODULE  CALCULATES  THE  ARRAY  N8S<«>  REPRESENTING  THE 
REM  BACKSCATTERED  CROSS  SECTION  NORMALIZED  BY  THE  CROSS  SECTIONAL  AREA 
REM  OF  THE  SPHERE. 

REM  ++++++ ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


!++++  MODULE  VARIABLE  DECLARATIONS 

REAL  L 

REAL  Lnax 

REAL  J(S1 ) 

REAL  Dj ( 5 1  ) 

« 

REAL  Y( 5  I  ) 

REAL  Dy ( 5 1  ) 

REAL  P( 5 1  ) 

REAL  Rsn 

REAL  Isn 

REAL  Fac 


AND  DEFINITIONS  +++++++++++++++++++++ 


ORDER  OF  THE  FUNCTION  CALCULATED 
IN  A  SUBROUTINE.  IT  IS  USED 
AS  A  LOOP  INDEX  AND  TO  DESIGNATE 
THE  ARRAY  ELEMENT  FOR  THE  FUNCTION 
OF  CORRESPONDING  ORDER. 

MAX.  SIGNIFICANT  ORDER  OF  THE  SUM 
DETERMINING  THE  DIFFERENTIAL 
SCATTERING  CROSS  SECTION  ( DSCS  ) 

AND  SCATTERING  MODULUS  (  SM > . 

ARRAY  FOR  THE  REGULAR  SPHERICAL 
BESSEL  FUNCTIONS  ( RSB  )  OF 
ARGUMENT  KA  AND  ORDER  0  THROUGH 
LMAX . 

ARRAY  FOR  THE  DERIVATIVES  OF  THE 
RSB  (THE  DRSB  ) . 

ARRAY  FOR  THE  IRREGULAR  SPHERICAL 
BESSEL  FUNCTIONS  ( ISB  )  OF 
AR6UMENT  KA  AND  ORDER  0  THROUGH 
LMAX. 

ARRAY  FOR  THE  DERIVATIVES  OF  THE 
ISB  (THE  DISB  ) . 

ARRAY  FOR  THE  LEGENORE 
POLYNOMIALS  OF  ARGUMENT  -I. 

THE  REAL  COMPONENT  OF  THE 
SCATTERING  MODULUS. 

THE  IMAGINARY  COMPONENT  OF  THE 
SCATTERING  MODULUS. 

FACTOR  COMMON  10  RSM  AND  ISM. 

+  +  +  +  NOTE  ++++ 

THE  ARRAY  DIMENSIONS  (30+21)  NEED 
TO  BE  INCREASED  TO  ACCOMODATE 
KAMAX>30 . 
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nee 
i  tie 
1120 
1  130 
1  140 
1  150 
I  160 
1  170 
1  180 
1190 
1200 
1210 
1220 
1230 
1240 
1 250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1  470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 


! 

l 

1++++++4+4+++++++++++++++4+4+++4+++++44+4+4++44+444+44+4444+44444+444444 

! 

!++++  CALCULATING  THE  NBC5  ARRAY  ++++4++++++++++++++++++++++++++++++++++ 

l 

l 

PRINTED  IS  1 

INPUT  Kamin.Kamax  !  DETERMINE  DESIRED  RANGE  OF  KA 

Kanin«INT(Kamln  )  !  UALUES. 

Kamax* INT( Kanax  ) 

Stp= .  1 

! 

FOR  Ka”Kamin+Stp  TO  Kamax  STEP  Stp 
Lmax“INT ( Ka+2  1  1 

CALL  Drsb(Lnax  ,Ka  ,Dj( •  )  ,J( *  )  1  !  CALCULATE  THE  DRSB  OF  ARGUMENT  KA 

!  ANO  ORDERS  0  THROUGH  LMAX . 

CALL  Disb(Lnax  .Ka  ,Oy( *  )  ,Y! •  )  )  !  CALCULATE  THE  DISB  OF  ARGUMENT  KA 

!  AND  ORDERS  0  THROUGH  LMAX. 

CALL  Leg180(Lmax  ,Ka  ,P< *  )  1  !  CALCULATE  P(»>  OF  ARGUMENT  -I  AND 

!  ORDERS  0  THROUGH  LMAX. 

Rsm»0  !  CALCULATING  THE  REAL  AND  IMAGINARY 

I  sn*=0  !  COMPONENTS  OF  THE  SCATTERING 

FOR  L-0  TO  Lmax  I  MODULUS. 

Fac*=(  Dj  (  L  )*  <  2*L  +  1  )*P<  L  )  )/(Dj(L  )*Dj(L  >  +  Dy(L  )*Dy<L  >  ) 

Rsm“Rsci+(  Fac»Dy ( L  ) ) 

Isrv=Ism4(  Fac*0  j  (  L  )  ) 

NEXT  L 

!  CALCULATING  THE  ARRAY  REPRESENTING 
!  THE  NBCS . 

Nbs(  <Ka-Kanln)/Stp  >“4«(Rsn»Rsn4l5(<(»Isn)/(Ka»Ka) 

I 

NEXT  Ka 

I 

! 

I  44  +  444444  +  44444444444444444444444  +  4444444444444444444444444444444444444 
| 

144+4  MOOULE  CONCLUSION  44+4++44+44+444444+4444+44+++444+4444++4444444+4 

I 

SUBEND 

! 

144++44+4444+4444++44+44+444+444++++444444+4+44++4++++444+++4+4+++4+4++4 

I 
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1 550  ! 

1 560  ! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 

1570  SUB  Drsb(Lma>.  ,Ka  ,0j<  •  )  ,J<  *  ) ) 

1 580  REfi  ++++++++++♦+++++++++ 4++++++++++++++++++ +++++++++++++++++++ ++++++++ 

1590  REM  THIS  MODULE  CALCULATES  THE  DERIVATIVES  OF  THE  RE6ULAR  SPHERICAL 

1600  REM  BESSEL  FUNCTIONS  ( DRSB  )  OF  ARGUMENT  KA  AND  ORDERS  0  THROUGH  LMAX 

16)0  REM  THE  CALCULATION  OF  THE  DERIVATIVE  INVOLVES  THE  RSB’S  OF  THE  SAME 

1620  REM  rtNO  PREVIOUS  ORDER. 

1630  REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

1640  I 

1650  14444  CALCULATING  THE  DERIVATIVES  OF  THE  RSB  +++++++++++++++++++++++++++ 

1660  ! 

1670  I 

1680  CALL  Rsb( Lmax ,Ka  ,  J <  «  )  >  !  CALCULATE  THE  RSB  VALUES  J <  *  > . 

1690  ! 

1700  RAD  i  ANGLE  UNITS  FOR  DRSB  CALC’NS. 

1710  Dj<0)=(-J(0)+COS(Ka))/Ka  !  CALCULATE  THE  INITIAL  DRSB  DJ<0). 

1720  FOR  L“ 1  TO  Lmax  !  CALCULATE  THE  REMAINING  DRSB. 

1730  Dj  <  L  )■-( ( L+ 1  )/Ka  >* J ( L )+J<  L-1  ) 

1740  NEXT  L 

1750  I 

1760  I 

1 770  ! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

1780  ! 

1790  !++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

1800  I 

1810  SUBEND 

1820  ! 

I  830  I 44444444444444444444444444444444444444444444444444444444444444444444444 
1840  ! 

1 850  ! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$93$ 

1860  ! 
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1870 
1880 
1880 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2330 
2340 
2350 
2360 
2370 
2380 
2390 
2400 
24  10 


I 

SUB  Rsb( Lmax  ,Ka  ,  J( *  )) 

REM  ++++++++++  +++■*-+++++++++  +  +  +  +++  +  +  ++++  ++++  +  +  ++++++++++++++  +++  +  +  +++  +  +  + 
REM  THIS  MODULE  CALCULATES  THE  RE6ULAR  SPHERICAL  BESSEL  FUNCTIONS 
REM  <  RSB  )  OF  ARGUMENT  KA  AND  ORDERS  0  THROUGH  LMAX .  THE 
REM  CALCULATION  USES  THE  CONTINUED  FRACTION  APPROACH. 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


!++++  MODULE  VARIABLE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 


REAL 

R1 

! 

1 

1  THE  EVOLVING  RATIO  OF  J(0)/J(L>. 

REAL 

Nuror 

1 

!  THE  EVOLVING  NUMERATOR  FACTOR  USED 
1  IN  CALCULATION  OF  RL. 

| 

REAL 

Denr 

!  THE  EVOLVING  DENOMINATOR  FACTOR 
!  USED  IN  CALCULATION  OF  RL. 

| 

REAL 

Scexp 

1 

!  EVOLVING  SCALING  EXPONENT  FOR  VERY 

I  LARGE  OR  VERY  SMALL  VALUES  OF  RL. 

| 

REAL 

Nf  lag 

!  0  FLA6  FOR  NUMR. 

j 

REAL 

Of  lag 

!  0  FLAG  FOR  DENR. 

1 

REAL 

Bni*i 

\ 

1  EVOLVING  TERM  USED  IN  THE 
!  CALCULATION  OF  NUMR  AND  DENR. 

| 

REAL 

Bine 

!  INCREMENT  USED  IN  THE  EVOLUTION 

1  OF  BNM. 

! 

1 

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


++++  CALCULATING  THE  RSB  +++++++++++++++++++++++++•+++++++++++++++++++++ 


RAD 


ANGLE  UNITS  FOR  CALC’NS. 


J(0  >  =  <  SIN(Ka  )  )/Ka 


««•  CALCULATE  THE  INITIAL 
»•«  RSB,  J(0>.  *•* 


FOR  L=1  TO  Lmax 
Rl-1  .0 
Numr=0 . 

Denr*=0 . 
Scexp=B  . 

Nf lag=0. 

Df lag»0. 

Bnm= 1 . 0/Ka 
P  :ic«2 . 0/Ka 


•••  CALCULATING  THE  RATIO  •  •• 
«»«J(0)/J(L). 

INITIALIZING  RL. 

INITIALIZING  NUMR. 
INITIALIZING  DENR. 
INITIALIZING  SCEXP. 
INITIALIZING  NFLAG. 
INITIALIZING  DFLAG . 
INITIALIZING  BNM. 

CALCULATE  BINC. 
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2420 

1  •  CALCULATE  THE  L  'NAKEO “ 

2430 

FOR  Ordcnt-1  TO  L 

1  «  NUMR  FACTORS  OF  J<0>/J(L>.  • 

2440 

Bnm-Bnm+Binc 

1  INCREMENT  BNM. 

2450 

IF  Nflag-0  THEN 

!  CHECK  NUMR-0  FLAG  NOT  SET; 

2460 

Numr-Bnm-Nunr 

1  UPDATE  NUMR  IF  NUMRO0 . 

2470 

IF  Nunr<>0  THEN 

2480 

Rl-Rl«Numr 

1  UPDATE  RL  IF  UPDATED  NUMRO0. 

2480 

Numr-1  .0/Nufnr 

I  PREPARE  NUMR  FOR  NEXT  EUOLUTION. 

2500 

1  CHECK  SCALING  OF  UPOATED  RL . 

2510 

IF  ABS! R1 )> 1 0*250 

THEN 

CALL 

Seal e( Ri  .Scexp > 

2520 

ELSE 

2530 

Nf lag-l 

2540 

END  IF 

2550 

ELSE 

!  IF  NFLAGO0  THE  PREUIOUS  STEPS  ARE 

2560 

Nf  lag-0 

I  OEFERREO  FOR  ONE  INCREMENT  OF  BNM. 

2570 

END  IF 

2580 

NEXT  Ordcnt 

2590 

1  *  CALCULATE  THE  REMAINING  NUMR  • 

2600 

1  •  ANO  DENR  TERMS  OF  J!0)/J!L>.  * 

2610 

(  NOTE  WHEN  NUMR-DENR  THEN 

2620 

1  RL-J ( 0  >/ J( L  ) . 

2630 

WHILE  NunrODenr  OR  Nflag-I  I 

DR  Dflag-1 

2640 

Bnei-Bnit+Binc 

1  INCREMENT  BNM. 

2650 

1 

2660 

IF  Nflag-0  THEN 

!  CHECK  NUMR-0  FLAG  NOT  SETi 

2670 

Nunr-Brim-Nunr 

!  UPDATE  NUMR  IF  NUMRO0 . 

2680 

IF  Nunr<>0  THEN 

2690 

Rl-Rl •Nunr 

!  UPDATE  RL  IF  UPDATED  NUMRO0. 

2700 

Nueir-I .0/Nunr 

1  PREPARE  NUMR  FOR  NEXT  EUOLUTION. 

2710 

1  CHECK  SCALING  OF  UPDATED  RL . 

2720 

IF  ABS!R1 >>10*250 

THEN 

CALL 

Scale! Rl , Scexp ) 

2730 

ELSE 

2740 

Nf lag-1 

2750 

END  IF 

2760 

ELSE 

1  IF  NFLAGO0  THE  PREUIOUS  STEPS  ARE 

2770 

Nf lag-0 

1  DEFERRED  FOR  ONE  INCREMENT  OF  BNM. 

2780 

END  IF 

2790 

! 

2800 

IF  Of  lag-0  THEN 

1  CHECK  DENR-0  FLAG  NOT  SET i 

2810 

Denr-Bnrc-'Oenr 

I  UPDATE  DENR  IF  DENRO0. 

2820 

IF  Denr<>0  THEN 

2830 

Rl-Rl/Denr 

1  UPDATE  RL  IF  UPDATED  DENRO0. 

2840 

Denr-1 .0/Oenr 

!  PREPARE  DENR  FOR  NEXT  EUOLUTION. 

2850 

!  CHECK  SCALING  OF  UPDATED  RL . 

2860 

IF  ABS! R 1  )> 10*250 

THEN 

CALL 

Sea  1 e ( R 1  .Scexp  > 

2870 

ELSE 

2880 

Of  1 ag= 1 

2890 

END  IF 

2900 

ELSE 

!  IF  DFLAGO0  THE  PREUIOUS  STEPS  ARE 

2910 

Of  lag-0 

!  DEFERRED  FOR  ONE  INCREMENT  OF  BNM. 

2920 

END  IF 

2930 

! 

2940 

END  WHILE 
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2950 

2950 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 


J(L  >-=J<0  )/<Rl*10''(Scexp  ) ) 


i 

I  J < L  >  DETERMINED.  *•• 


NEXT  L 

I 

j 

!+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

j 


!++++  MODULE  CONCLUSION 
SUBEND 


I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
1 
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3120 

3130 

3140 

3IS0 

3160 

3170 

3180 

3190 

3200 

3210 

3220 

3230 

3240 

3250 

3260 

3270 

3280 

3290 

3300 

3310 

3320 

3330 

3340 

3350 

3360 

3370 

3380 

3390 

3400 

3410 

3420 


I 

$•$$$$$$$$$$$$$$$ 

SU8  Scale(Rl  ,Scexp  ) 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM  THIS  MODULE  SCALES  RL  TO  MAINTAIN  NUMERICAL  ACCURACY  IN  THE 
REM  ONGOING  CALCULATION  OF  RL. 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
! 

I++++  SCALING  RL  +++++++++++4-f++++++++++++++++++++++++++++++++++++++++++ 

I 

! 

IF  Rl>10'250  THEN 
R1“R1 *  1 0A( -250  > 

ScexpijScexp+250 
END  IF 

IF  Rl<10',(-250  >  THEN 
R1»R1 • 1 0"250 
Scexp“Scexp-2S0 
END  IF 

! 

I 

! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
! 

!++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

! 

SUBEND 

1 

! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

I $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$*$$$$$$$$$$$$$$$$$$ 

I 
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SUB  DlablLmax ,Ka  ,0y( •  )  ,Y<  * ) ) 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM  THIS  MODULE  CALCULATES  THE  DERIVATIVES  OF  THE  IRREGULAR  SPHERICAL 
REM  BESSEL  FUNCTIONS  OF  ARGUMENT  KA  AND  ORDERS  0  THROUGH  LMAX .  THE 
REM  CALCULATION  OF  THE  DERIVATIVE  USES  THE  ISB'S  OF  THE  SAME  AND 
REM  PREVIOUS  ORDER. 

REM  ++++++++++++++++++++++♦+++++++++++++++++++++++++++++++++++++++♦+++ 


++++  CALCULATING  THE  DERIVATIVES  OF  THE  ISB  +++++++++++++++++++++++++++ 

I 
I 

CALL  Isb( Lmax ,Ka , Y( •  ) )  I  CALCULATE  THE  ISB  VALUES. 

I 

AD  I  ANGLE  UNITS  FOR  CALC’NS. 

y(0)“(-Y(0)+SIN(Ka))/Ka  I  THE  INITIAL  ISB  DERIVATIVE. 

OR  L-l  TO  Lmax  I  THE  REMAINING  ISB  DERIVATIVES. 


RAD 

Dy(0)«(-Y(0)+SIN(Ka))/Ka 
FOR  L“1  TO  Lmax 

Dy ( L  )“-( ( L+ 1  )/Ka  )«Y(L  )  +  Y( L-l ) 
NEXT  L 


! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
! 

!++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

SUBEND 

I 

I++++++++++++++++++++++4++++++++++++++++++++++++++++4+++++++++++++++++++ 


! $$$S$$$$$$$$$$$$$$$$$SS$$SS$$$$SS$$$$$$$$$S$$$$$$$$$$$$$$$$$$SS$$$$$$£$ 


3750  * 

3760  l**S$$**$$*$S$SSSS$$$SS**$$$$$$$$$SSS$$$$$$S$S$$$$S«$S$$$$$$$S*$$S*S#$S* 

3770  SUB  IabfLmax  ,Ka  ,Y( »  ) ) 

3780  REM  ++++++++++++++++*+++++++++++++++++++++++++++++++++++++++++++++++++ 

3790  REM  THIS  MODULE  CALCULATES  THE  IRREGULAR  SPHERICAL  BESSEL  FUNCTIONS 

3800  REM  OF  ARGUMENT  KA  AND  OROER  L  THROUGH  LMAX .  I,.'  CALCULATION  IS 

3810  REM  ACCOMPLISHED  BV  FORWARD  RECURSION. 

3820  REM  ++++++++++++++++++++++++++++++ ++++++++++++++++++++++++++ ++++++++++ 

3830  I 

3840  !++++  CALCULATING  THE  ISB  +++++++++++++++++++++++++++++++++++++++++++♦++ 

3850  ! 

3860  ! 

3870  RAO  !  ANGLE  UNITS  FOR  CALC ’NS. 

3880  Y(0)»-COS(Ka>/Ka  !  INITIAL  ISB’S  CALCULATED. 

3890  Y( 1  >“<  Y( 0  >-SIN( Ka  )  >/Ka 

3900  FOR  L-2  TO  Lnax  !  REMAINING  ISB’S  CALC’D. 

3910  Y(  L  )»<  (  2«L-1  )/Ka  ) •»  Y <  L—  1  )-Y<  L-2  > 

3920  NEXT  L 

3930  ! 

3940  ! 

3950  ! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

3960  ! 

3970  !++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

3980  ! 

3990  SUBEND 

4000  ! 

4010  I +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
4020  ! 

4040  i 
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4050  ! 

4060  ! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$*$$$$$$$$$$$$$$$$$$$$$$$$* 

4070  SUB  Leg  1 80( Lmax  ,Ka  ,P( «  )  > 

4080  REM  +++++++++++++++++++++++++++ t+++++++++f +++++++++ ++++++++++ +++++++++ 

4090  REM  THIS  MODULE  CALCULATES  THE  LEGENDRE  POLYNOMIALS  OF  ARGUMENT  -1  AND 

4100  REM  ORDER  0  THROUGH  LMAX .  CALCULATION  IS  ACCOMPLISHED  BY  FORWARD 

4(10  REM  RECURSION. 

4 1 20  REM  +++++++++++++■+++++++++++++++++++++++++++++++++++++++++++++++++++++ 

4130  I 

4140  !+  +  +  +  CALCULATING  THE  LEGENDRE  POLYNOMIALS  +++++++++++++++++++++++++++++ 

4150  ! 

41G0  X  — 1  1  ARGUMENT  OF  THE  LEG.  POL’S. 

4170  ! 

4180  P(0)-1  !  INITIAL  LEG.  POL’S.  CALC’D. 

4190  P< 1  )“X 

4200  FOR  L=2  TO  Lnax  !  REMAINING  LEG.  POL’S.  CALC’D. 

4210  P<L  )«<2-< 1/L  >  >*X*P(L-1  )-<)-( 1/L  )  >»P(L-2  > 

4220  NEXT  L 

4230  ! 

4240  ! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

4250  ! 

4260  !++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

4270  I 

4280  SUBEND 

4290  ! 

4300  ! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

4310  ! 

4320  ! $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$£$$$&$$$$$$$$$ 

4330  I 
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4340 

4350 

4360 

4370 

4380 

4390 

4400 

4410 

4420 

4430 

4440 

4450 

4460 

4470 

4480 

4490 

4500 

4510 

4520 

4530 

4540 

4550 

4560 

4570 

4580 

4590 

4600 

4610 

4620 

4630 

4640 

4650 

4660 

4670 

4680 

4690 

4700 

4710 

4720 

4730 

4740 

4750 

4760 

4770 

4780 

4790 

4800 

4810 

4820 

4830 

4840 

4850 

4860 

4870 

4880 


SUB  Recpl t<  Nbs( *  )  , Kan In  .Kanax  ,Stp  ) 

REM  +++++■+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM  THIS  MODULE  MAPS  THE  CALCULATED  SCATTERING  MODULUS  ARRAY  SMC)  ON 
REM  A  SEMI -LOG  PLOT. 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

!++++  MODULE  VARIABLE  DECLARATIONS  AND  DEFINITIONS  444444444444444444444 


1 

! 


REAL 

Xgumax 

!  THE  MAXIMUM  ABSCISSA  VALUE 
!  GRAPHICS  DISPLAY  UNITS. 

( 

REAL 

Ygunax 

1  THE  MAXIMUM  ORDINATE  VALUE 
1  GRAPHICS  DISPLAY  UNITS. 

44444444444444444444444444i 


14444  INITIALIZING  THE  PLOTTER  44444444444444444444444444444444444444444 

l 

GINIT 

PLOTTER  IS  CRT  .“INTERNAL" 

GRAPHICS  ON 
GCLEAR 


I 


! 


14444  GRAPH  CONSTRUCTION  AND  LABELLING  44444444444444444444444444444444 

i 

I 

Xgunax»100«MAX<  1  .RATIO) 

Ygumax” 1 00* MAX ( 1  ,1/RATIO) 

j 

CSIZE  4 
LINE  TYPE  1 

I 

LORG  4  !  ««  LABELLING  THE  ABSCISSA  AXIS  •* 

MOVE  Xgumax/2  .0 
LABEL  "  K  «  A " 

! 

OEG  !  ..  LABELLING  THE  ORDINATE  AXIS  ** 

LOIR  90 
LORG  6 

MOVE  0  ,Ygunax/2 

LABEL  “NORM’D  BCKSCTR6  CROSS  SEC’N" 

LDIR  0 


VIEWPORT  10  .Xgunax -5 , 10  ,Ygunax-5 


*•  GRAPH  CONSTRUCTION  «• 

DEFINES  THE  TOTAL  PLOTTING  AREA 
IN  GRAPH  UNITS. 


WINDOW  Kanin .Kanax  ,0 , 1 . 2 
LINE  TYPE  4 

GRID  Stp»10,.l,Kanln,0,61,10,.2 
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4890  I 

4900  CLIP  OFF 
4910  CSIZE  3 
4920  LINE  TYPE  1 

4930  ! 

4940  LORG  6  !  •*  LABELLING  THE  ABSCISSA  SCALE  •« 

4950  FOR  I=Kamtn  TO  Kamax  STEP  5tp*20 

49G0  MOVE  I ,0 

4970  LABEL  I 

4980  NEXT  I 

4990  ! 

5000  LORG  8  !  *«  LABELLING  THE  ORDINATE  SCALE  •« 

5010  FOR  1-0  TO  I  STEP  ,5 

5020  MOVE  Kamin  ,1 

5030  LABEL  I 

5040  NEXT  I 

5050  I 

5060  CLIP  ON 

5070  ! 

5080  ! 

5090  !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

5100  i 

5110  !++++  PLOTTING  THE  CALCULATED  NBCS  +++++++++++++++++++++++++++++++++++++ 

5120  •  1 

5130  FOR  I-Kamin+Stp  TO  Kamax  STEP  Stp 
5140  PLOT  I  ,Nbs(< I-Kanin  )/Stp > 

5150  NEXT  I 

5160  ! 

5170  i ++++++++++++++++++++++++++++++++++++++++++++++++++++++++ +++++++++++++++ 
5180  ! 

5190  !++++  PRINTING  THE  PLOT  ++++++++++++++++++++++++++++++++++++++++++++++++ 

5200  ! 

5210  DUMP  GRAPHICS  4701 

5220  I 

5230  I +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
5240  ! 

5250  !++++  MODULE  CONCLUSION  ++++++++++++++++++++++++++++++++++++++++++++++++ 

5260  I 

5270  GRAPHICS  OFF 
5280  SUBEND 

5290  ! 

5300  I +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
5310  ! 

S320  !$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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ACOUSTIC  ARRAY  EFFICIENCY  PROGRAM 


This  program  was  used  to  analyze  data  and  determine  the  calibration 
product  Er  Et.  The  output  of  this  program  was  also  used  to  construct 

Figures  13,  14,  and  15. 

The  power  supply  mentioned  in  the  Gain30  and  Gain35  subprograms 
was  the  HP33 1 4A  function  generator.  All  data  contained  in  the  data 
subprograms  are  repeated  in  Appendix  E.  Data  units  are  available  in 
Appendix  E  with  their  respective  data. 


te 

20 

30 

40 

50 

G0 

70 

80 

90 

100 

130 

140 

150 

IG0 

170 

180 

190 

200 

210 

220 

230 

240 

250 


I 

REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
! 

!++++  VARIABLE  DECLARATIONS  AND  DEFINITIONS 


++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
ACOUSTIC  ARRAY  EFFICIENCY 

THIS  PROGRAM  CALCULATES  THE  EFFICIENCY  AND  CENTERLINE  GAIN  OF 
AN  ACOUSTIC  ARRAY  BASED  ON  THE  RETURN  FROM  AN  ACOUSTICALLY  HARD 
SPHERE  OF  DIAMETER  D  ENSONIFIED  BY  AN  INCIDENT  PLANE  WAVE  OF 
FREQUENCY  F. 


H  +  +  +  + 


REAL  Vmsr 


MEAN  SQUARE  VOLTAGE  RECEIVED  FROM 
THE  TARGET  SPHERE  AND  MEASURED 
AFTER  PRE-AMPLIFICATION. 


REAL  Vnsn 


MEAN  SQUARE  VOLTAGE  RECEIVED  IN 
THE  ABSCENCE  OF  A  TARGET  AND 
MEASURED  AFTER  PRE-AMP. 
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MEAN  SQUARE  VOLTAGE  TRANSMITTED. 


260 

270 

REAL 

Vnst 

| 

1 

MEAN  SQUARE  VOLTAGE  TRANSMITTED. 

280 

290 

300 

310 

320 

REAL 

R 

I 

1 

1 

I 

1 

THE  RANGE  FROM  THE  PLANE  OF  THE 
ARRAY  SPEAKER  DIAPHRAGMS  TO  THE 
CENTER  OF  THE  TARGET  SPHERE  AND  TO 
THE  SPECTRUM  ANALYZER  MICROPHONE. 

330 

340 

350 

360 

REAL 

Et 

j 

| 

I 

j 

THE  EFFICIENCY  OF  CONVERSION 

FROM  TRANSMITTED  ELECTRICAL  POWER 

TO  TRANSMITTED  ACOUSTICAL  POUER. 

370 

380 

390 

REAL 

Go 

! 

| 

| 

THE  CENTERLINE  GAIN  OF  THE  ACOUSTIC 
ARRAY. 

400 

410 

REAL 

Etgo 

| 

j 

ET«GO 

420 

430 

440 

REAL 

At 

? 

I 

j 

THE  ATTENUATION  COEFFICIENT, 
"ALPHA". 

450 

460 

REAL 

C  .Cmoist 

i 

j 

THE  SPEED  OF  SOUND  IN  HUMID  AIR. 

470 

480 

490 

REAL 

Ef  req 

f 

| 

THE  SELECTED  FREQUENCY  OF  THE 
ACOUSTIC  ARRAY. 

500 

510 

REAL 

K 

I 

l 

THE  WAVE  NUMBER.  2*PI*EFREQ/C . 

520 

530 

REAL 

A 

i 

1 

THE  RADIUS  OF  THE  TARGET  SPHERE. 

540 

550 

560 

570 

REAL 

Ka 

! 

! 

f 

| 

K*A ,  THE  ARGUMENT  OF  THE  REGULAR 

AND  IRREGULAR  SPHERICAL  BESSEL 
FUNCTIONS  AND  THEIR  DERIVATIVES. 

580 

590 

REAL 

Aa 

! 

! 

THE  APERATURE  AREA  OF  THE  ARRAY. 

600 

610 

REAL 

Gpre 

i 

j 

THE  PRE-AMP’S  ELECTRICAL  GAIN. 

620 

630 

640 

650 

REAL 

Nbs 

l 

| 

i 

f 

THE  VARIABLE  REPRESENTING  THE 
NORMALIZED,  BACKSCATTERED  CROSS 
SECTION  ( NBCS  ) . 

660 

670 

680 

REAL 

Atgt 

i 

i 

1 

CROSS  SECTIONAL  AREA  OF  THE 

TARGET  SPHERE. 

690 

700 

710 

720 

730 

740 

REAL 

Sasqb 

; 

i 

i 

1 

i 

1 

SURFACE  AREA  OF  AN  IMAGINARY 

SPHERE  CENTERED  ON  THE  ARRAY 

AND  HAVING  A  RADIUS  EQUAL  TO 

THE  RANGE  TO  THE  CENTER  OF  THE 
TARGET  SPHERE. 

750 

760 

770 

780 

REAL 

Prat 

i 

i 

! 

/ 

THE  RATIO  OF  RECEIVED  TARGET 
ELECTRICAL  POWER  TO  THE 

TRANSMITTED  ELECTRICAL  POWER. 

THE  EFFICIENCY  OF  CONVERSION 
FROM  RECEIVED  ACOUSTICAL  POWER 
TO  RECEIVED  ELECTRICAL  POWER. 


790 
800 
810 
820 
830 
840 
850 
860 
870 
880 
890 
900 
910 
920 
930 
940 
950 
960 
970 
980 
990 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1  100 
1110 
I  120 
1  130 
1  140 
1150 
1160 
1  170 


REAL  E r 


++++  MAIN  PR06RAM 


! 

I 


CALL  Ini t  <  Vnsr , Oman  ,Vnst  ,R,Etgo,At  ,A,Ka,Aa ,Gpre ) 

!  INPUT  AND  CALCULATE  REQUIRED 
!  PARAMETERS. 

i 

CALL  Bcksctr(Ka.Nbs)  1  CALCULATE  THE  NBCS. 

PRINT  *NBS«-;Nbs 

i 

Atgt-PI*A«A 

PRINT  -AT6T-- jAtgt 

Sasph»4»PI»R»R 

Prat“<  Vmsr-Vitsn  )/(  6pre*Gpre*Vmst  ) 

Er“(Prat )»<Sasph»Sasph/Etgo  >*<EXP<2*At»R  >/(Aa»Atgt  >  >/Nbs 

I 

PRINT  "Er-'iEr 
PRINT  “Et6o“" sEtgo 
PRINT  *Go-‘ lEtgo/Er 
PRINT  "EtEr-“ sEr*Er 

i 

i 


++++++++-M-+++++++++ 


++++  MAIN  PROGRAM  CONCLUSION 


i 


PRINT  CHRS( 12  ) 

PRINTER  IS  I 
PRINT 

PRINT  "END"  '  ADVISES  THE  USER  OF  THE  MAIN 

PRINT  1  PROGRAM’S  CONCLUSION. 

END 


1180 

1190  !+++++++++++++++++++++++++* 

1200  I 
>210  1 $$$$$$$ 

1220  ! 

1230  ' 


M-  +  +  +  +  +  ++  +  +  -M-+  +  +  +  -M 


I  !  ! 


1240 
1250 
I  260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1 690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 


SUB  Ini t  <  Vnsr .Unsn.Vmst  ,R  ,E tgo  ,At ,A  ,Ka , Aa ,6pre  ) 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM  THIS  MOOULE  DETERMINES  THE  VALUES  OF  THE  PARAMETERS  REQUIRED  FOR 
REM  THE  CALCULATION  OF  Er . 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

I 

!++++  MODULE  VARIABLE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 


REAL  Ch 

REAL  Vr( 60 ) 

REAL  J  ,Jt  , Jn , J tn 

REAL  Tc 
REAL  Patn 

REAL  Rh 

REAL  Tgtln 

REAL  Dia 

REAL  Nuspkr 

REAL  Diaspkr 

REAL  Vt  <  60  ) 

REAL  Nol se<  60  > 

REAL  Tnoi9e(60) 

REAL  Ten 
REAL  Rhn 


A  CONTROL  VARIABLE  USED  FOR 
CONTROLLING  SELECTION  OF  THE 
CORRECT  GAIN  DATA  SETS. 

THE  ARRAY  REPRESENTING  THE 
MEASURED  TARGET  RETURN  VOLTAGES. 

THE  NUMBER  OF  ARRAY  ELEMENTS  MINUS 
ONE  AND  A  LOOP  CONTROL  VARIABLE. 

THE  AMBIENT  TEMPERATURE  IN  OEG.  C. 

THE  ATMOSPHERIC  PRESSURE  MEASURED 
IN  MILLIBARS. 

THE  RELATIVE  HUMIDITY  MEASURED  IN 
PERCENT. 

THE  MEASURED  TIME  FOR  THE  ACOUSTIC 
SIGNAL  TO  MAKE  THE  ROUND  TRIP 
BETWEEN  THE  ARRAY  AND  THE  TARGET. 

THE  MEASURED  DIAMETER  OF  THE 
TARGET  SPHERE. 

THE  NUMBER  OF  SPEAKERS  IN  THE 
ARRAY. 

THE  AVERAGE  DIAMETER  OF  AN  ARRAY 
SPEAKER. 

THE  ARRAY  REPRESENTING  THE 
MEASURED  TRANSMISSION  VOLTAGES. 

THE  ARRAY  REPRESENTING  THE 
MEASURED  NOISE  RETURN  VOLTAGES. 

THE  ARRAY  REPRESENTING  THE 
MEASURED  TRANSMISSION  VOLTAGES  FOR 
NOISE  DETERMINATION. 

Tc  FOR  NOISE  DETERMINATION. 


Rh  FOR  NOISE  DETERMINATION. 


1790 

1800 

REAL 

Patnn 

!  Patn  FOR  NOISE  DETERMINATION. 

I 

1810 

1820 

REAL 

Efreqn 

1  E  f reqn^E f req  . 

I 

1830 

1840 

REAL 

Cn 

!  C  FOR  NOISE  DETERMINATION. 

| 

1850 

I860 

REAL 

Atn 

!  At  FOR  NOISE  DETERMINATION. 

j 

1870 

1880 

REAL 

Umstn 

!  Vmst  FOR  NOISE  DETERMINATION. 

! 

1890 

1900 

1910 

1920 

REAL 

Rspkr 

!  THE  AVERAGE  RADIUS  OF  AN  ARRAY 
!  SPEAKER. 

I 

j 

1 930  i  ++++++++++++++++ +++  +  +  +++++++++++++■•*■+++++ ++++++++  +++++++++++ ++++++++  +++++ 
1940  ! 

1950  !++++  DETERMINING  SYSTEM  PARAMETERS  +++++++++++++++++++++++++++++++++++++ 

I960  ! 

1970  ! 

1980  PRINTER  IS  1 

1990  PRINT  "WHICH  TARGET  DATA  SET 

2000  PRINT  "DO  YOU  WISH  ANALYZED?" 

2010  PRINTER  IS  701 

2020  INPUT  Sphere 

2030  PRINT  "TARGET  NO.  "iSphere 

2040  PRINT 

2050  ! 

2060  IF  Sphere3 1  THEN 

2070  Ch*=3 . 0 

2080  CALL  Gaino< Ch  ,Etgo  ) 

2090  CALL  Sphere 1< Ur («),J1Tc, Pa tn.Rh.Tgttn.Dia.Efreq.Nuspkr, Diasp kr.Gpre) 

2100  CALL  Trans30( J t  ,Ut ( •  )  > 

2110  CALL  Noise30<Jn,JtntNoise(»),Tnoise(#),Tcn,Rhn,Patnn .Efreqn ) 

2120  END  IF 

2130  ! 

2140  IF  Sphere=2  THEN 

2150  Ch=3 . 5 

2160  CALL  Gaino( Ch ,EtQO  ) 

2170  CALL  Sphere2( Ur( *  )  ,J  ,Tc  ,Patn,Rh,Tgttm,Dia  ,E f req  .Nuspkr .Diaspkr ,6pre ) 

2180  CALL  Trans35( J t  ,Ut < •  )  ) 

2  I  90  CALL  No ise35(Jn,Jtn,  No lse(«)tTnoise<«),  Ten,  Rhn.Patnn, Efreqn) 

2200  END  IF 

2210  ! 

2220  IF  Sphere-3  THEN 

2230  Ch=3 . 5 

2240  CALL  Ga i no ( Ch  ,E t qo  ) 

2250  CALL  Sphere3( Ur( »  )  ,  J  ,Tc  .Patn  .Rh.Tgttn  ,Di  a  ,E  f req  .Nuspkr , Diaspkr .Gpre ) 

2260  CALL  Trans 35 <  Jt ,Vt( •  )  ) 

2270  CALL  Noise35(Jn,Jtn,Noise(»),Tnoise(»  >, Ten, Rhn.Patnn, Efreqn) 

2280  END  IF 
2290 


I  13 


I 


2300 

2310 

2320 

2330 

2340 

23S0 

2360 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2580 

2590 

2600 

2610 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 


IF  Sphere-4  THEN 
Ch-3.5 

CALL  6aino( Ch ,£tgo ) 

CALL  Sphere4<  Vrl  *  ) ,J ,Tc ,Pat«  ,Rh  ,Tgt  tm  ,Dia  .Efreq .Nuspkr .Oiaspkr ,Gpre  > 
CALL  Trans35< J t ,Vt < *  ) ) 

CALL  Noise35( Jn.Jtn.NoiselO.TnoisefO.Tcn.Rhn.Patmn.Efreqn) 

END  IF 

i 

CALL  Voltsmsl Vr< •  ) ,  J  ,Vmsr  ) 

PRINT  "VMSR»“ ;Vmsr 

CALL  Atten<  Ten ,Rhn , Pat nn.Cn.Atn, Efreq n> 

CALL  Vo It sms ( Tnoisel *  > , Jtn .Vmstn ) 

CALL  Voltsms(Noise( ♦ ) , Jn  ,Vmsn  ) 

CALL  Voltsmsl Vt< * >  ,Jt  ,Vmst  ) 

PRINT  'VMST-'iVnst 

! 

A-Dia/2 
PRINT  "A»" ) A 

CALL  Range ( Tgttm,A,Tc ,Rh,Patm,C,At  .Efreq.R) 

PRINT  “C=*,-C 

PRINT  “AT=“ ; At 

PRINT  " ATTEN=" ; EXP( At *R  ) 

PRINT  “R=* iR 

Vmsn-Vmsn»£XP<  2*R*(  Atn-At ))»Vmst/Vmstn 

PRINT  "VMSN-" jVmsn 

Ka=2*PI*Efreq*A/C 

PRINT  “KA” * . Ka 

Rspkr=0iaspkr/2 

Aa=Nuspkr*PI *Rspkr»Rspkr 

PRINT  “AA=“;Aa 


+  +  ^ 

++++  MODULE  CONCLUSION 
SUBEND 


! ++++++++++++++++H 


I-+++  +  +  +  +++  +  H 


I 


2720 

2730 

2740 

27S0 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2890 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 

3120 

3130 

3140 

3150 


!$$$$$S$$$$$$$$$$$$S$$$$$$$$$$$$SSS$SS$$S$S$$$$SS$S$$?$S$$SS$$$SS$$$S$$S 

SUB  6alno(Ch  ,Etgo  ) 

REM  THIS  MODULE  DETERMINES  THE  VALUE  OF  Et6o. 

!++++  MODULE  VARIABLE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 


ROOT  MEAN  SQUARE  VQLTA6E  MEASURED 
IN  SAIN  DETERMINATION. 

MEASURED  INTENSITY  AT  R; 
Ptgt/Atgt. 

ACOUSTIC  ARRAY  IMPEDANCE. 


REAL  Urns 


REAL  Ir 


REAL  Za 


++++  DETERMINATION  OF  EtGo 


IF  Ch-3.0  THEN  CALL  6ain30< Urns  ,Tc .Rh  ,R .Patn . Ir ,Ef req .Za > 

IF  Ch=3 . 5  THEN  CALL  Gain3S< Urns  .Tc  ,Rh  ,R  ,Patn ,Ir .Efreq ,Za > 

!  OBTAINING  UALUES  FROM  APPROPRIATE 
!  DATA. 

t 

CALL  Atten(Tc .Rh. Patn, C, At  .Efreq  )  !  CALCULATING  AT. 

Pt “Urns* Urns /Za 

Etgo-Ir*4»PI*R*R*EXP< At*R  )/Pt  !  CALCULATING  EtGo. 


I 

SUBEND 

i 

! +++++++++++++++++++++++++++++++++++♦+++++++++++++++++++++++++++++++++++ 

! 

!XS$SSSS$$$SSS$SS$$$SSS$SS$SSSSSS$$SSSSS$S$$SSS$$SSSSSS$S$SSSSSSSSSS$SSS 


3160 

3170 

3180 

3190 

3200 

3210 

3220 

3230 

3240 

3250 

32G0 

3270 

3280 

3290 

3300 

3310 

3320 

3330 

3340 

3350 

3360 

3370 

3380 

3390 

3400 

3410 

3420 

3430 

3440 

3450 

34G0 

3470 

3480 

3490 

3500 

3510 

3520 

3530 

3540 

3550 

3560 

3570 

3580 

3590 

3600 

3610 

3620 

3630 

3640 


SUB  Voltsms( Vol t  <  * ) , J ,Vns ) 

REM  THIS  MODULE  DETERMINES  THE  MEAN  SQUARE  VALUE  OF  THE  INPUT  VOLTASE 
REM  ARRAY. 

REM 
| 

!++++  MODULE  VARIABLE  DECLARATIONS  AND  DEFINITIONS 


(REAL  VOLT ( • ) 

REAL  Volts 
REAL  VavQ 
REAL  V 

IREAL  Vns 

++++++++^ 

++++  DETERMINATION  OF  VMS 


Volt  s=0 
FOR  1=0  TO  J 

Volts-Volts+Volt( I  ) 
NEXT  I 

Vavg=Volts/<  J+1  ) 

FOR  1=0  TO  J 
V=Vol t ( I  )-Vavg 
Vms  =  ( V*V  )+Vns 
NEXT  I 

Vns=Vms/<  J+1  ) 


+++++++++++++++ +H 

++++  MODULE  CONCLUSION  +++h 
SUBENO 


1  THE  ARRAY  REPRESENTING  THE  INPUT 
VOLTAGE  ARRAY. 

THE  SUM  OF  VOLT < • > * S  ELEMENTS. 

THE  AVERAGE  VALUE  OF  VOLT(»). 

THE  DIFFERENCE  BETWEEN  VOLT( I  ) 
AND  VAVG. 

1  THE  MEAN  SQUARE  VALUE  OF  VOLTO. 


K++++++++++++++ ++++++++++++++++++ 


I-  +  +  +  +  +  +  +  H 


►  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  ++  +  +  +  +  +  +  -*■  +  +  +  +  +  +  +  +  +  +  +  ++  +  +  +  +  +  +  +  +  +  +  + 


3650  i 
3660  • 


3G70 

3680 

3690 

3700 

3710 

3720 

3730 

3740 

3750 

3760 

3770 

3780 

3790 

3800 

3810 

3820 

3830 

3840 

3850 

3860 

3870 

3880 

3890 

3900 

3910 

3920 


i 


SUB  Rangel Tgttn, A, Tc.Rh,  Pa  tn.C.At.Efreq.R) 

REM  THIS  MODULE  DETERMINES  THE  AVERA6E  PATH  ATTENUATION  AND  THE  RANGE 
REM  BETWEEN  THE  ARRAY  AND  THE  TARGET  CENTER. 

REM 
! 

!++++  DETERMINATION  OF  R  AND  AT 


Tgt  tm=Tgt  tm/ 1 000 

CALL  AttenITc ,Rh,Patm,C.At  .Efreq  ! 
R»Tgt  tm»C/2+A 


AT  AND  C  DETERMINED. 
R  DETERMINED. 


!++++  MODULE  CONCLUSION 
SUBEND 


I 


117 


3930 

3940 

3950 

3960 

3970 

3980 

3990 

4000 

4010 

4020 

4030 

4040 

4050 

4060 

4070 

4080 

4090 

4100 

4110 

4120 

4130 

4140 

4150 

4160 

4170 

4180 

4190 

4200 

4210 

4220 

4230 

4240 

4250 

4260 

4270 

4280 

4290 

4300 

4310 

4320 

4330 


SUB  At  ten ( Tc ,Rb , Pa tn .Crooi s t , At .Efrep ) 

REM  THIS  MODULE  CALCULATES  THE  ATTENUATION  USING  RELATIVE  HUMIDITY 
REM  IN  PERCENT,  ATMOSPHERIC  PRESSURE  IN  MILLIBARS,  ANO  TEMPERATURE 
REM  IN  DEGREES  CENTIGRADE. 

!++++  MODULE  VARIA8LE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 


REAL  Tk 


REAL  Es 


REAL  Pratio 


REAL  Cdry 


REAL  H 


REAL  Past 


REAL  Tast 


REAL  Fn 


REAL  Anax 


REAL  Frat 


REAL  F2 


REAL  Acl  ,Amol 


THE  AMBIENT  TEMPERATURE  IN  DEG.  K. 

THE  SATURATION  VAPOR  PRESSURE. 

THE  RATIO  OF  WATER  PRESSURE  IN 
MILLIBARS  TO  ATMOSPHERIC  PRESSURE 
IN  MILLIBARS. 

THE  SPEED  OF  SOUND  IN  DRY  AIR. 


P  ASTERISK. 


T  ASTERISK. 


MAX.  FRED. 


MAX  EXPECTED  ATTENUATION. 


EFREQ/FM 


FRAT*FRAT 


ATTENUATION  EXPONENT  COMPONENTS. 


>+++++++++ 


4340 

4350 

4360 

4370 

4380 

4390 

4400 

4410 

4420 

4430 

4440 

4450 

4460 

4470 

4480 

4490 

4500 

4510 

4520 

4530 

4540 

4550 

4560 

4570 

4S80 

4590 

4600 

4610 

4620 

4630 

4640 

4650 

4660 

4670 

4680 

4690 

4700 

4710 


!  ++++  CALCULATING  THE  ATTENUATION 


Tk-Tc+273.15 


Rh«Rh/100 

Es^lO'M  9.4-2353/Tk  ) 
Pratio«Es*Rh/(Patm-Es*( 1+Rh)1 

Cdry-20 . 05*SQR<  Tk  ) 

Cnoj  st «Cdry*< 1  +  . 1 4«Pra 1 10  1 


CONVERSION  OF  TEMP  IN  CELSIUS  TO 
TEMP  IN  KELVIN . 

CONVERSION  OF  RH  IN  PCT.  TO  PRATI0 . 


CALCULATION  OF  THE  SPEED  OF  SOUND. 


H»Pratio*100 
Past-Patm/1014 
Tast“< 1 . 8*Tc+492 )/519 

Fn=( 1 0+6600 *H+44400»H*H  1»Past/(Tast“.8 ) 

Amax-.0078*Fm*<TastA<-2.5>  >»EXP<7.77»( 1-1/Tast  )) 

Frat-Efreq/Fn 

F2=Frat*Frat 

Amol-An?x«SOR(  .  1  8* .  1  B«F2+<  2»F2/<1 +F2  1  T2  1/304.8 
Acl  =  1 .74*I0“< -10  >*Efreq*Efreq 
At  =  <  Amol+Acl )»LOG< 101/10 

1 


++++  MODULE  CONCLUSION 
SUBEND 


i  ■ 
i 

I 


4720 

4730 

4740 

4750 

4760 

4770 

4780 

4790 

4800 

4810 

4820 

4830 

4840 

4850 

4860 

4870 

4880 

4890 

4900 

4910 

4  920 

4930 

4940 

4950 

4960 

4970 

4980 

4990 

5000 

5010 

5020 

5030 

5040 

5050 

5060 

5070 

5080 

5090 

5100 

51  10 

5120 

5130 

5140 

5150 

5160 

5170 

5180 

5190 

5200 


I 

!*$$$*$* ****$$$$$$$$$$$$$$$*$«$ 

SUB  Bcksctr< Ka  ,Nbs  ) 

REM  + ++++++++++++++++++++++++++++++++++++++++++++++++++++++ +++++++++4.+ 

REM  THIS  MODULE  CALCULATES  THE  BACKSCATTERED  CROSS  SECTION  NORMALIZED 
REM  BY  THE  CROSS  SECTIONAL  AREA  OF  THE  SPHERE. 

REM  +++++++++++++++++++++++++++++++++++++++++++++ ++++++ +++++++++++++++ 

! 

I++++  MODULE  VARIABLE  DECLARATIONS  AND  DEFINITIONS  +++++++++++++++++++++ 


! 

! 


REAL 

L 

!  ORDER  OF  THE  FUNCTION  CALCULATED 

1  IN  A  SUBROUTINE.  IT  IS  USED 
!  AS  A  LOOP  INDEX  AND  TO  DESIGNATE 

1  THE  ARRAY  ELEMENT  FOR  THE  FUNCTION 
!  OF  CORRESPONDING  ORDER. 

I 

REAL 

Lnax 

•  MAX.  SIGNIFICANT  ORDER  OF  THE  SUM 
)  DETERMINING  THE  NORMALIZED 
!  BACKSCATTERING  CROSS  SECTION 
!  ( NBCS  ) . 

| 

REAL 

J<  5 1  ) 

!  ARRAY  FOR  THE  REGULAR  SPHERICAL 

1  BESSEL  FUNCTIONS  (RSB)  OF 
!  ARGUMENT  KA  AND  ORDER  0  THROUGH 
!  LMAX . 

| 

REAL 

D  j  (  5  1  ) 

!  ARRAY  FOR  THE  DERIVATIVES  OF  THE 
!  RSB  (THE  DRSB  ) . 

I 

REAL 

Y(S1  ) 

!  ARRAY  FOR  THE  IRREGULAR  SPHERICAL 
!  BESSEL  FUNCTIONS  (ISB)  OF 

1  ARGUMENT  KA  AND  ORDER  0  THROUGH 

1  LMAX. 

| 

REAL 

Dy<  51  ) 

!  ARRAY  FOR  THE  DERIVATIVES  OF  THE 
!  ISB  (THE  DISB  ) . 

| 

REAL 

P(  51  ) 

!  ARRAY  FOR  THE  LEGENDRE 
!  POLYNOMIALS  OF  ARGUMENT  COSU80). 

| 

REAL 

Fac 

1  FACTOR  COMMON  TO  RSM  AND  ISM. 

1 

REAL 

Rsm 

!  THE  REAL  COMPONENT  OF  THE 
!  SCATTERING  MODULUS,  SM. 

| 

REAL 

I  SCI 

!  THE  IMAGINARY  COMPONENT  OF  THE 
!  SM. 

| 

120 


++++  NOTE  +++  + 

THE  ARRAY  DIMENSIONS  (30+21)  NEED 
TO  BE  INCREASED  TO  ACCOMODATE 
KA>30. 


++++  CALCULATING  THE  N8CS  +++++++++++++++++++++++++++++++++++++++++++++ 


Li»iax“INT ( Ka+2 1 ) 


CALL  Orsb< Lmax ,Ka ,D j < »  )  , J< * > ) 


CALL  Di sb<  Lmax  ,Ka  ,DyC  )  #YC  ) ) 


CALL  Leg! 80 ( Lmax ,Ka  ,P( *  ) ) 


CALCULATE  THE  DERIVATIVES  OF  THE 
REGULAR  SPHERICAL  BESSEL  FUNCTIONS 
OF  ORDER  0  THR0U6H  LMAX . 

CALCULATE  THE  DERIVATIVES  OF  THE 
IRREGULAR  SPHERICAL  BESSEL  FUNC’NS 
OF  ORDER  0  THR0U6H  LMAX. 

CALCULATE  THE  ARRAY  PC). 


Rsm=0  !  CALCULATING  TH 

lsm=0  !  COMPONENTS  OF 

FOR  L=0  TO  Lmax 

Fac=< Dj < L  >*(2*L  +  I  )*P(L  )  )/(Dj(L  >*Dj(L )+Dy(L )*Dy(L ) ) 
Rsn“Rsm+(  Fac*Dy  (  L ) ) 

Ism=*Ism+(Fac*Dj(L)) 

NEXT  L 


CALCULATING  THE  REAL  AND  IMAGINARY 
COMPONENTS  OF  THE  SM. 


Nbs“4*(Rsm*R5m+Ism*Ism)/(Ka*Ka> 


CALCULATING  THE  NBCS. 


! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
I 

! 

SUBEND 

I 

! 


12! 


5630 

5640 

5650 

5660 

5670 

5680 

5690 

5700 

5710 

5720 

5730 

5740 

5750 

5760 

S770 

5780 

5790 

5800 

5810 

5820 

5830 

5840 

5850 

5860 

5870 

5880 

5890 

5900 

5910 

5920 

5930 

5940 


SUB  Drsb( Lmax ,Ka  ,D j ( *  ) , J< *  ) ) 


REM 

REM 

REM 

REM 

REM 

REM 

! 


THIS  MODULE  CALCULATES  THE  DERIVATIVES  OF  THE  REGULAR  SPHERICAL 
BESSEL  FUNCTIONS  ( DRSB )  OF  ARGUMENT  KA  AND  ORDERS  0  THROUGH  LMAX. 
THE  CALCULATION  OF- .THE  DERIVATIVE  INVOLVES  THE  RSB’S  OF  THE  SAME 
ANO  PREVIOUS  ORDER. 


!++++  CALCULATING  THE  DERIVATIVES  OF  THE  RSB 

I 


CALL  Rsb<  Lmax  ,Ka  ,  J<  *  )  > 

RAD 

0 j  <  0  )“< -J( 0  )+COS<  Ka  )  >/Ka 
FOR  L-1  TO  Lmax 

D j < L )«-<  ( L+ 1 )/Ka)*J<L)+J(L-l  ) 
NEXT  L 


CALCULATE  THE  RSB  VALUES  JO. 

ANGLE  UNITS  FOR  DRSB  CALC’NS. 
CALCULATE  THE  INITIAL  DRSB  OJ<0). 
CALCULATE  THE  REMAINING  DRSB. 


i 


++++  MODULE  CONCLUSION 
SUBEND 
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5950 
5960 
5970 
5980 
5990 
6000 
6010 
6020 
6030 
6040 
6050 
6060 
6070 
6080 
6090 
6100 
61  10 
6120 
6130 
6140 
6150 
6160 
6170 
6180 
6190 
6200 
6210 
6220 
6230 
6240 
6250 
6260 
6270 
6280 
6290 
6300 
63)0 
6320 
6330 
6340 
6350 
6360 
6370 
6380 


! 

SUB  RsbILmax ,Ka , J( *  ) ) 

REM  THIS  MODULE  CALCULATES  THE  RE6ULAR  SPHERICAL  BESSEL  FUNCTIONS 
REM  ( RSB  )  OF  ARGUMENT  KA  AND  ORDERS  0  THROUGH  LMAX .  THE 
REM  CALCULATION  USES  THE  CONTINUED  FRACTION  APPROACH. 


!++++  MODULE  VARIABLE  DECLARATIONS  AND  DEFINITIONS 


REAL 

R1 

1 

I 

!  THE  EVOLVING  RATIO  OF  J(0)/J(L). 

i 

REAL 

Numr 

!  THE  EVOLVING  NUMERATOR  FACTOR  USED 
!  IN  CALCULATION  OF  RL . 

i 

REAL 

Oenr 

!  THE  EVOLVING  DENOMINATOR  FACTOR 
!  USED  IN  CALCULATION  OF  RL. 

j 

REAL 

Scexp 

!  EVOLVING  SCALING  EXPONENT  FOR  VERY 
•  LARGE  OR  VERY  SMALL  VALUES  OF  RL. 

I 

REAL 

Nf  lag 

!  0  FLAG  FOR  NUMR. 

t 

REAL 

Df  lag 

!  0  FLAG  FOR  DENR. 

i 

REAL 

Brin 

!  EVOLVING  TERM  USED  IN  THE 
!  CALCULATION  OF  NUMR  AND  DENR. 

I 

REAL 

Bine 

!  INCREMENT  USED  IN  THE  EVOLUTION 
!  OF  BNM. 

! 

I 

RAO 


ANGLE  UNITS  FOR  CALC’NS. 


J ( 0  >*( SIN( Ka  )  >/Ka 


***  CALCULATE  THE  INITIAL 

RSB ,  J( 0  ) .  *«• 
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6390 

6400 

6410 

6420 

6430 

6440 

6450 

6460 

6470 

6460 

6490 

6500 

6510 

6520 

6530 

6540 

6550 

6560 

6570 

6580 

6590 

6600 

6610 

6620 

6630 

6640 

6650 

6660 

6670 

6680 

6690 

6700 

6710 

6720 

6730 

6740 

6750 

6760 

6770 

6780 

6790 

6800 

6810 

6820 

6830 

6840 

6850 

6860 

6870 

6880 

6890 

6900 

6910 

6920 


FOR  L=  1  TO  Lr-iax 
R 1  - 1 .0 
Nunr=0 . 

Denr=0 . 
Scexp=0 . 

Nf lag«0 . 

Of lag“0. 

Bnn“ 1 . 0/Ka 
Binc=2 . 0/Ka 


FOR  Ordcnt=1  TO  L 
Bnrt“Bnn+Binc 
IF  Nf lag=0  THEN 
Nunr'Bnm-Numr 
IF  Nunr<>0  THEN 
RI=R1 »Nunr 
Nunr»  1  . 0/Nufnr 

IF  ABS< R1  )> 1 0"250  THEN  CALL 
ELSE 
Nf lag»1 
END  IF 
ELSE 
Nf lag*0 
END  IF 
NEXT  Ordcnt 


WHILE  NunrOOenr  OR  Nflag=l  OR 
8nn=Bnn+Binc 


!  CALCULATING  THE  RATIO 

!  J( 0  )/ J( L  ) . 

!  INITIALIZING  RL . 

!  INITIALIZING  NUMR . 

!  INITIALIZING  DENR. 

!  INITIALIZING  SCEXP. 

1  INITIALIZING  NFLAG. 

!  INITIALIZING  OFLAG. 

!  INITIALIZING  BNM. 

I  CALCULATE  BINC. 

I 

!  *  CALCULATE  THE  L  “NAKED"  * 

!  ♦  NUMR  FACTORS  OF  J(0)/J(L).  * 

I  INCREMENT  BNM. 

!  CHECK  NUMR=0  FLAG  NOT  SETi 
!  UPDATE  NUMR  IF  NUMRO0. 

i  UPDATE  RL  IF  UPDATED  NUMRO0. 

I  PREPARE  NUMR  FOR  NEXT  EVOLUTION. 
!  CHECK  SCALING  OF  UPDATED  RL. 
Scale( R1  .Scexp  ) 


I  IF  NFLAGO0  THE  PREVIOUS  STEPS  ARE 
!  DEFERRED  FOR  ONE  INCREMENT  OF  BNM. 


!  *  CALCULATE  THE  REMAINING  NUMR  • 
!  •  AND  DENR  TERMS  OF  J(0)/J(L>.  • 
!  NOTE  WHEN  NUMR-DENR  THEN 
!  RL»J(0)/J(L>. 

ag=  1 

!  INCREMENT  BNM. 


IF  Nf lag*=0  THEN 
Nunr=Bnn-Nunr 
IF  Nunr<>0  THEN 
Rl'Rl »Nunr 
Nunr» 1 . 0/Nunr 

IF  ABS( R1  )> 1 0 "250  THEN  CALL 
ELSE 


I 

I  CHECK  NUMR=0  FLAG  NOT  SETs 
!  UPDATE  NUMR  IF  NUMRO0. 

!  UPDATE  RL  IF  UPDATED  NUMRO0. 

!  PREPARE  NUMR  FOR  NEXT  EVOLUTION. 
!  CHECK  SCALING  OF  UPDATED  RL . 
Scalef  R I  .Scexp  ) 


Nf lag= 1 
END  IF 
ELSE 

Nf lag»0 
ENO  IF 

IF  Df lag=0  THEN 
Oenr=Bnm-Denr 
IF  Denr <>0  THEN 
R1 "R 1 /Denr 
Oenr- 1 . 0/Oenr 


l  IF  NFLAGO0  THE  PREVIOUS  STEPS  ARE 
1  DEFERRED  FOR  ONE  INCREMENT  OF  BNM. 


I  CHECK  DENR=0  FLAG  NOT  SET; 
i  UPDATE  DENR  IF  DENRO0. 

i  UPDATE  RL  IF  UPDATED  DENRO0. 

I  PREPARE  DENR  FOR  NEXT  EVOLUTION. 


» 
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£930 
6940 
6950 
6860 
6970 
6980 
6990 
7000 
7010 
7020 
7030 
7040 
7050 
7060 
7070 
7080 
7090 
7100 
71  10 
7120 
7130 
7140 
7150 
7160 
7170 
7180 
71  90 


!  CHECK  SCALING  OF  UPDATED  RL. 
IF  ABS<R1 >>10'2S0  THEN  CALL  Scale* R1  .Scexp  ) 

ELSE 


!  IF  DFLAGO0  THE  PREVIOUS  STEPS  ARE 
'  DEFERRED  FOR  ONE  INCREMENT  OF  BNM. 


**♦  J < L  >  DETERMINED.  *♦* 


Of  lag'! 
END  IF 
ELSE 
Of lag“0 
END  IF 

END  WHILE 


J(L  )«J< 0 )/< Rl» 1 0A( Scexp  ) ) 
NEXT  L 

+++  +  +  •4 

++++  MODULE  CONCLUSION 
SUBEND 


7200 

72) 0 
7220 
7230 
7240 
7250 
72B0 
7270 
7280 
7290 
7300 

73) 0 
7320 
7330 
7340 
7350 
7360 
7370 
7380 
7390 
7400 
7410 
7420 
7430 
7440 
7450 
7460 
7470 
7480 
7490 
7500 


I 


I$SSSSSS$SS*$S$**$SS$$SS$$S$$S$$S$SSSS$$$$$$$$$$$$$$$S$$$$$$$S$$S$$$S$$$ 

SUB  Scale( R1  .Scexp  ) 


REM 

REM 

REM 

REM 


THIS  MODULE  SCALES  RL  TO  MAINTAIN  NUMERICAL  ACCURACY  IN  THE 
ONGOINS  CALCULATION  OF  RL. 


i 


'++++  SCALING  RL 


IF  R 1 > 1 0*250  THEN 
R1»R1*I0*(-2S0> 
Scexp**Scexp+250 
END  IF 

IF  R1 < 1 0* ( -250  )  THEN 
R1-R1»10*2S0 
Scexp«Scexp-250 
END  IF 


M +++++++ 


!++++  MODULE  CONCLUSION 
SUBENO 


! 
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7510 

7520 

7530 

7540 

7550 

75G0 

7570 

7580 

7590 

7B00 

7610 

7620 

7630 

7640 

7650 

7660 

7670 

7680 

7690 

7700 

7710 

7720 

7730 

7740 

7750 

7760 

7770 

7780 

7790 

7800 

7810 

7820 


! 


SUB  Disb( Lnax  ,Ka,Dy(«),Y<»>) 

REM  THIS  MODULE  CALCULATES  THE  DERIUATIUES  OF  THE  IRRE6ULAR  SPHERICAL 
REM  BESSEL  FUNCTIONS  OF  ARGUMENT  KA  AND  ORDERS  0  THROUGH  LMAX .  THE 
REM  CALCULATION  OF  THE  DERIVATIVE  USES  THE  ISO’S  OF  THE  SAME  AND 
REM  PREVIOUS  ORDER. 

REM 
! 

!++++  CALCULATING  THE  DERIVATIVES  OF  THE  ISB 


CALL  Isb( Lnax  ,Ka  ,  Y(  *  ) ) 

RAD 

Dy(0)=<-Y<0)+SIN(Ka))/Ka 
FOR  L-1  TO  Lnax 

Dy( L  )«*-( ( L+ 1  )/Ka  )*Y(L)+Y(L-1  ) 
NEXT  L 


CALCULATE  THE  ISB  VALUES, 

ANGLE  UNITS  FOR  CALC’NS. 

THE  INITIAL  ISB  DERIVATIVE. 

THE  REMAINING  ISB  DERIVATIVES. 


i 


+++++++++++ 
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7830 
7840 
7850 
7860 
7870 
7880 
7890 
7900 
7910 
7920 
7930 
7940 
7950 
7960 
7970 
7980 
7990 
8000 
8010 
8020 
8030 
8040 
8050 
8060 
8070 
8080 
8090 
8100 
81  10 
8120 


i 


!  $$$$$$$$ $$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$ 


SUB 

REM 

REM 

REM 

REM 

REM 


Isb(  Leiax  ,Ka  ,Y<  »  )  ) 


h++++++++++++4 


h  +  +  + 


THIS  MODULE  CALCULATES  THE  IRRE6ULAR  SPHERICAL  BESSEL  FUNCTIONS 
OF  ARGUMENT  KA  AND  ORDER  L  THROUGH  LMAX .  THE  CALCULATION  IS 
ACCOMPLISHED  BY  FORWARD  RECURSION. 


!++++  CALCULATING  THE  ISB 


RAD 

Y<0  )—  COSIKa  )/Ka 
Y ( 1  )=< Y(0)-SIN(Ka  )  )/Ka 
FOR  L=2  TO  Lnax 

Y<  L  )■( < 2*L-1  )/Ka  )*Y< L-1  )-Y( L-2  ) 
NEXT  L 


!  ANGLE  UNITS  FOR  CALC’NS. 

!  INITIAL  ISB’S  CALCULATED. 

!  REMAINING  ISB’S  CALC’D. 


i 


i 


K+++++H 


++++  MODULE  CONCLUSION 
SUBEND 


K  +  +  T  +  '4‘  +  -f  +  +  4 


b++++++++++ 


$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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8130 

8140 

8150 

8160 

8170 

8180 

8190 

8200 

8210 

8220 

8230 

8240 

8250 

8260 

8270 

8280 

8290 

8300 

8310 

8320 

8330 

8340 

8350 

8360 

8370 

8380 

8390 

8400 

8410 


SUB  Leg  1 80<  Lmax  ,Ka  ,P< ♦ > ) 

REM  THIS  MODULE  CALCULATES  THE  LESENDRE  POLYNOMIALS  OF  ARSUMENT  -1  AND 
REM 
REM 
REM 
! 

!++++  CALCULATING  THE  LEGENDRE  POLYNOMIALS 


THIS  MODULE  CALCULATES  THE  LESENDRE  POLYNOMIALS  OF  ARSUMENT  -1 
ORDER  0  THR0U6H  LMAX .  CALCULATION  IS  ACCOMPLISHED  BY  FORWARD 
RECURSION. 


X=- I  1  ARGUMENT  OF  THE  LEG.  POL’S. 

I 

P(0)=1  !  INITIAL  LES.  POL’S.  CALC’D. 

P< 1 )=X 

FOR  L=2  TO  Lnax  !  REMAINING  LEG.  POL’S.  CALC’D. 

P(L)=(2-( 1/L)>*X*P(L-1  >-( 1-< 1/L))*P<L-2> 

NEXT  L 


i 

! 


129 


130 


SUB  Gain35( Vrms  ,Tc,Rh,R,Patm,Ir,EfreQ,Za> 

REM  THIS  MODULE  CONTAINS  THE  DATA  USED  FOR  DETERMINATION  OF  Et6o  WITH 
REM  3.5V  INPUT  FROM  THE  POWER  SUPPLY. 

i 


RESTORE  S35 

READ  Vrns  .Tc.Rh.R.Patm.Ir ,E  freq ,Za 
635:  DATA  .523 ,21 .0  ,50.0,5.30 . 1021 .5 .4. 467E-3 .5000 . 15. 8 


++++  MODULE  CONCLUSION 


SUBEND 
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8940 

8950 

8960 

8970 

8980 

8990 

9000 

9010 

9020 

9030 

9040 

9050 

9060 

9070 

9080 

9090 

9100 

9110 

9120 

9130 

9140 

9150 

9160 

9170 

9180 

9190 

9200 

9210 

9220 

9230 

9240 

9250 

9260 

9270 


SUB  Sphere) <0t< •  )  ,J ,Tc ,Patn  ,Rh  ,  Tgttm  ,Dia .Efreq .Nuspkr.Diaspkr ,Gpre  ) 

REM  THIS  MODULE  CONTAINS  THE  DATA  USED  FOR  DETERMINATION  OF  Er  WITH 
REM  TARGET  NO.  I  EMPLOYED. 

REM 
! 

!++++  LOADING  DATA 


J=39 

RESTORE  SI 
FOR  1-0  TO  J 
READ  Dt<I) 

NEXT  I 

READ  Tc,Patm,Rh,Tgttn,Dia  ,Efreq  .Nuspkr.Diaspkr.Gpre 
Sis  DATA  9. 54,9. 4,8. 32.6. 42. 3. 88. .96. -2. 07, -4. 9, -7. 24, -8. 87 

DATA  -9.64,-9.45.-8.32,-6.41  ,-3 . 83  .- . 89  ,2 . 1 4 ,4 . 94 .7 . 26  .8 . 88 
DATA  9.6 ,9.38,8.22 ,6.26  ,3. 67,. 73, -2. 3, -5. 11  ,-7.4  ,-8.96 
DATA  -9.64, -9. 39,-8. 2, -6. 21  ,-3.61  ,-.65,2. 3S  ,5. 12,7.38,8.91 
DATA  20. 2,1016. 8, 52. 8. 29. 91  .  .2546 .5000 ,19,. 0762 ,11094 

I 

i 


++++  MODULE  CONCLUSION 


!■  +  ++  +  +  +  +  +  +  +  +  +  ■( 


I 


SUBEND 
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3280 

9290 

9300 

93) 0 
9320 
9330 
9340 
9350 
9360 
9370 
9380 
9390 
9400 

94) 0 
9420 
9430 
9440 
9450 
9460 
9470 
9480 
9490 
9500 
9510 
9520 
9530 
9540 
9550 
9560 
9570 
9580 
9590 
9500 
96)0 


SUB  Sphere2(0t(*),J,Tc, Patm, Rh,Tgttm,Dia,Efreq,Nuspkr,Diaspkr,Gp re) 

REM  THIS  MOOULE  CONTAINS  THE  DATA  USED  FOR  DETERMINATION  OF  Er  WITH 
REM  TARGET  NO.  2  EMPLOYED. 

! 

!++++  LOADING  DATA 


J-39 

RESTORE  S2 
FOR  1=0  TO  J 
READ  Dt  < I > 

NEXT  I 

READ  Tc,Patf>i,Rh,Tgttn,0ia,Efreq,Nuspkr,Diaspkr,Gpre 
S2:  DATA  4 . 72 ,4  .  S6 .4 . 1 6  ,3 . 23  .2 . 0 1  . . 560 . 940  ,-2 . 35  .-3 . 52  ,-4 . 36 
DATA  -4.77 .-4.70 ,-4. 1 8 .-3.25  .-) . 99  .-.560  . . 960  ,2 .34  ,3 . 5 1  .4.34 
DATA  4.74 ,4.66 ,4. 13 .3. t 9,1 .95  .  .500  .-) .00  .-2.40  ,-3.57,-4.38 
DATA  -4.78,-4.68 ,-4. 15  ,-3.20 .-) .92  ,-.480,1 .02  ,2.39 ,3.54 ,4.34 
DATA  2 1 . 0 , 1 009 . 4 ,49 . 8 , 30 . 55  , . ! 0 1 3 , 5000 , 1  9 , . 0762 ,11094 


++++  MODULE  CONCLUSION 


SUBEND 


! 
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SUB  Sphere3<  Dt( *  )  ,J  ,Tc  ,Patm  ,Rh  ,  Tgttm  ,Dia .Efreq ,Nuspkr,Oiaspkr,Gpre) 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM  THIS  MODULE  CONTAINS  THE  DATA  USED  FOR  DETERMINATION  OF  Er  UITH 
REM  TARGET  NO.  3  EMPLOYED. 

! 


J-39 

RESTORE  S3 
FOR  I“0  TO  J 
READ  Dt  < I > 

NEXT  I 

READ  Tc.  Pa tm.Rh, Tgttm, Dia  .Efreq  .Nuspkr.Diaspkr.Gpre 
S3:  DATA  3 .6 1  ,3 . 40 .2 .86 .2 .04 . 1 .01 1 20 .-1 . 25 .-2 . 25  .-3. 04  .-3. 54 
DATA  -3.68,-3.46  ,-2.9)  .-2.07.-! .04.. 1 10.) .25.2.24,3.00,3.48 
DATA  3. 63, 3. 39,2. 84. 2. 02.. 870. -.180.-1 .30,-2.31  .-3.08.-3.56 
DATA  -3.69,-3.46  ,-2.88  ,-2.04  .-.990 ,. 160,1 .29,2.26 ,3.05 .3.50 
DATA  21 .3, 1009.6  ,49.8  ,30.59,-07634 ,5000 , 1 9 , . 0762 , 1 1094 


| ++++++++++++4++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
! 

! 

SUBEND 
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9960 
9970 
9980 
9990 
10000 
10010 
10020 
10030 
10040 
10050 
10060 
10070 
10080 
10090 
10100 
101  10 
10120 
10130 
10140 
10150 
10160 
10170 
10180 
10190 
10200 
10210 
10220 
10230 
10240 
10250 
10260 
10270 
10280 
10290 


! 

j  $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$ 

SUB  Sphere4<  Dt( *  )  ,J  ,Tc  ,Patm  ,Rh,Tgttm,Dia  .Efreq  .Nuspkr.Diaspkr.Gpre) 

REH  THIS  MODULE  CONTAINS  THE  OATA  USED  FOR  DETERMINATION  OF  E r  WITH 
REM  TARGET  NO.  4  EMPLOYED. 

REM  ++++++++++++++++H 
| 

!++++  LOADING  DATA 


J=39 

RESTORE  S4 
FOR  1=0  TO  J 
READ  Dt  < I  ) 

NEXT  I 

READ  Tc .Patm.Rh ,Tgt  tro  ,Dia .Efreq .Nuspkr.Diaspkr.Gpre 


S4 :  DATA  2.28,2. 10,1 .68 ,1 , 
DATA  -2.34 ,-2.12,-1.71 
DATA  2.29,2.10,1 .68,1 . 


10.. 420, -.320,-1 .03,-1 .64,-2.07,-2.31 
,-! . 12, -.440,. 320, 1 .02,1 .62,2.06 ,2.28 
09,.  400,-.  340,-1..  04,-1  .66,-2. 10,-2.32 
DATA  -2.3S.-2. 12  ,-1 .70  ,-1 . 1 1  ,-.400  ,  .340 , 1 .04 , 1 .64  ,2.08,2.28 
DATA  21 .4,1009.9,49.6  ,30 . 64  ,  .06240  ,5000 , 1 9 , . 0762 , 1 1094 

! 

I 


!++++  MODULE  CONCLUSION 
SUBEND 


i 


++++++++++++++++++++++T-M 
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10300 
10310 
10320 
10330 
10340 
10350 
10360 
10370 
10380 
10390 
10400 
10410 
10420 
10430 
10440 
10450 
10460 
10470 
10480 
10490 
10500 
10510 
10520 
10530 
10540 
10550 
1  0560 
10570 
10580 
10590 
10600 
10610 


SUB  Trans30<  J  ,Dt  <  *  )) 

REM  +++++++++++++++++++++++++  + <-+++++■►+++++++++++++++++++++++++++++++++ 
REM  THIS  MODULE  CONTAINS  DATA  USED  IN  DETERMINATION  OF  Er  WITH 
REM  3.0U  APPLIED. 

REM 


I 


!++++  LOAOING  DATA 


J=39 

RESTORE  T30 
FOR  1*0  TO  J 
READ  Dt< I > 

NEXT  I 

T30:  DATA  .690 , . 680  . . 590  , . 440  . . 320  , . 1 60  . 020  . 230  . 4 1 0  .550 

DATA  -.610  ,-.6 10. -.520  .-.360, -.240  ,-.080 . . 1 10. .320,. 490 ..630 
DATA  .690  .  .680  ..590  . . 430  , . 320  , . 1 50  .-.040  .-.240  .-.420  ,-.560 
DATA  - . 620 .- . 600 ,- . 5 1 0  . 380 .- . 240  .- . 070 . . 1 20  , . 330  , . 500  , . 640 


++++  MODULE  CONCLUSION 
SUBEND 


!  ++++++++++++H 


!$$$$ $$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ 
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10620 

10630 

10640 

10650 

10660 

10670 

10680 

10690 

10700 

10710 

10720 

10730 

10740 

10750 

10760 

10770 

10780 

10790 

10800 

10810 

10820 

10830 

10840 

10850 

10860 

10870 

10880 

10890 

10900 

10910 

10920 

10930 


I 

! $$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$ 

SUB  Trans3S<  J  ,Dt(»  ) ) 

REM  THIS  MODULE  CONTAINS  DATA  USED  IN  DETERMINATION  OF  Er  WITH 
REM  3.SU  APPLIED. 

REM 

! 

!++++  LOAOING  DATA 

l 

! 

J-39 

RESTORE  T35 
FOR  1=0  TO  J 
READ  Dt  < I  ) 

NEXT  I 

T35:  DATA  . 830  .  .770  , . 640  . . 490  ,  .290  . . 080  . 1 80  400 .6 1 0 .730 

DATA  - . 770  ,  - . 720  . 580  . 430  , - . 240 , - . 020  , . 240  . . 460  , . 660 . . 780 

DATA  .830, .770,. 630. .470, .290, .060, -.190, -.420, -.600, -.720 
DATA  -.770, -.7 10. -.570, -.420, -.230, -.010,. 240.. 480,. 670,. 790 

! 

i 

»■++++++++++++++++++++ 


++++  MODULE  CONCLUSION 

! 

SUBEND 

! 

i 
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10940 
10950 
10950 
10970 
10980 
10990 
11000 
11010 
1(020 
11030 
1  1040 
1 1050 
11060 
11070 
1 1080 
11090 
11100 
11)10 
11120 
11130 
11140 
11150 
1  1  160 
11170 
11180 
11190 
11200 
1  1210 
11220 
11230 
11240 
11250 
11260 
11270 
11280 
1  1290 
1  1300 
11310 
1  1320 
1  1330 
1  1340 
1  1350 
1  1360 
1  1  370 
11380 


i 


! $$$$$$$$$$$$$$$$$$$$$$$$$$$ 

SUB  Noise30<  J,Jt,Dt<*),Tdt<*),Tcn,Rhn,Patein,Efreqn> 

REM  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
REM  THIS  MODULE  CONTAINS  DATA  USED  IN  DETERMINATION  OF  Er  WITH 
REM  3.0V  APPLIED. 


!++++  LOADING  DATA 


J»59 

RESTORE  N30 
FOR  1=0  TO  J 
READ  Ot(I  ) 

NEXT  I 
J  t=39 

RESTORE  Tn30 
FOR  1=0  TO  Jt 
READ  Tdt ( I  ) 

NEXT  I 

READ  Ten  ,Rhn  .Patnn  .Efreqn 

N30:  OATA  0  .0  .  02  .03  .  04  .  06  .  08  .  1  .  1  1  .  1  2 

DATA  -.12, -.11  ,-.1  ,-.08, -.07, -.05, -.04, -.01  ,0,0 
DATA  0.0 ,-.02, -.03, -.05, -.07, -.08,-.!  ,-.11 ,-.12 
DATA  -.12.-. 12,-. 1 1  ,-.08, -.07, -.04, -.03, -.02 ,0,0 
OATA  .01  ,0.0,-. 03,-. 04  .-.07  ,-.09,-. 11  ,-. I  2,-.  12 
DATA  -.12, -.12, -.11  ,-.08, -.07, -.04, -.03, 0,0,. 01 
Tn30:  DATA  .69 . .68 . .59 , . 44  ,  .32  , .  16  ,-.02  ,-.23  .- .41  . 55 

DATA  -.61  .-.61  ,-.52, -'.36  ,-.24, -.08..  11  ,.32..  49..  63 
DATA  .69,. 68  ,. 59,. 43.. 32.. 15, -.04, -.24, -.42 ,-.56 
DATA  -.62, -.6, -.51  ,-.38  .-.24  ,-.07  . . 1 2  ,  .33  ,  .5 , .64 
DATA  20.8,51 .7.1016.7,5000 


i 


H++++++++ ++++++++++++++ 


I++++  MOOULE  CONCLUSION 
SUBEND 


i 


' +++++++ 


i 


+++++H 
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11390 
11400 
11410 
11420 
I  1430 
1  1440 
1  1  4S0 
11460 
11470 
11480 
11490 
1  1500 
11510 
1  1520 
1  1530 
11540 
1  1550 
1  1560 
1  1570 
1  1580 
1  1590 
1  1600 
11610 
1  1620 
1  1630 
1  1640 
11650 
11660 
11670 
11680 
11690 
1  1700 
1  1710 
11720 
1  1730 
1  1740 
1  1750 
11760 
1  1770 
11780 
1  1790 
1  1800 
11810 


SUB  Noise35<J,Jt,Dt(»),Tdt<»>,Tcn,Rhn,Patmn ,Ef reqn ) 

REM  THIS  MODULE  CONTAINS  DATA  USED  IN  DETERMINATION  OF  Er  WITH 
REM  3.5V  APPLIED. 

REM 
! 

!++++  LOADING  DATA 

! 

j 

J-39 

RESTORE  N35 
FOR  1=0  TO  J 
READ  Ot(I) 

NEXT  I 
Jt=39 

RESTORE  Tn3S 
FOR  1-0  TO  Jt 
READ  Tdt  < I > 

NEXT  I 

READ  Ten  ,Rhn  ,Patmn ,Ef reqn 

N3S:  DATA  -.06 ,-.08 .-.08 ,-.08  .-.07. -.06  ,-.04  ,-.02  ,0  , . 

DATA  . 05  , . 06  , . 06  , . 06  . . 04  , . 03  . . 02  ,0  , - . 03 . - . 05 
DATA  -.07. -.08. -.09, -.08  .-.08. -.06  ,-.04  ,-.02 ,.01 
DATA  .04 ,.07,. 06 ,. 07,. 05.. 04  ,.01  .-.01  .-.04  .-.06 
Tn35:  DATA  . 74 . . 83 . . 82 , . 72 . . 53  ,  .38  , . 1 8  ,- . 06 .- .3 1  ,- .53 
DATA  -. 68, -.77, -.75, -.64, -.46,-. 3, -.1  ,.14, .4, .6 
DATA  .75, .83, .81  ,  .7  ,  .5  ,  .36  , .  16  ,- .08  ,-.33  ,- .54 
DATA  -.69, -.76  .-.74, -.63  ,-.46  ,-.28, -.08,. 16 ,.41  , 

DATA  21 .3  ,51 .3 ,1016.6  ,5000 

! 

i 


.03 


.62 


++++  MODULE  CONCLUSION 
SUBEND 


+++++++++++ ++++-H 


H-  +  +  + 
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ACOUSTIC  ARRAY  CALIBRATION  DATA 


This  appendix  contains  data  used  to  determine  the  calibration 
product  Er  Et.  Data  units  are  contained  in  square  brackets. 

For  listings  of  the  same  type  of  data  the  units  and  uncertainty  are 
explicitly  given  for  the  first  datum.  Remaining  data  in  the  listing  have  the 
same  units  and  uncertainty  as  the  first  datum. 

1 .  TWO  WAY  PROP  AG  AT  ION  PATH  DATA 

Pre-amplif  jer  gain  Gfl  and  array  impedance  1  were  previously 

determined  by  Moxcey.  G8was  determined  to  be  1 1,094.  Z  was  determined 
to  be  15.8  [01.  [Ref.  8] 

The  remaining  data  are  presented  in  the  following  tables. 

^supply  is  vo,ta9e  supplied  by  the  HP33 14A  function  generator  to 


the  pre-amplifier  (see  Figure  8). 


Table  5a.  Two  Way  Propagation  Path  Data  for  Sphere  One, 
Collected  30  November  1987;  Run  *2:  Data  Re-measured 
Because  First  Run  Vr  (with  3.5  [V]  Supply)  Was  Clipped. 


Diameter  =  2a  =  .2546  +  0.0006  [m] 

Tc  =  20.2  ±0.1  [°  c] 

Rh  *  52.8  ±  0. 1  [%) 

P=  1016.8  ±  0.2  (mb] 

No.  pulses  =  20 
Vsupply'^M 

tr  =  30.31  [mS]  -  400  [jjlS]  =  29.91  i  0.1  [mS] 


t 

t 

Vt 

2.25  ±  0.01  [mS] 

690  ±  5  [mV] 

2.45 

690 

2.26 

680 

2.46 

680 

2.27 

590 

2.47 

590 

2.28 

440 

2.48 

430 

2.29 

320 

2.49 

320 

2.30 

160 

2.50 

150 

2.31 

-  20 

2.51 

-40 

2.32 

-230 

2.52 

-240 

2.33 

-410 

2.53 

-420 

2.34 

-550 

2.54 

-560 

2.35 

-610 

2.55 

-620 

2.36 

-610 

2.56 

-600 

2.37 

-520 

2.57 

-510 

2.38 

-360 

2.58 

-380 

2.39 

-240 

2.59 

-240 

2.40 

-  80 

2.60 

i 

o 

2.41 

1 10 

2.61 

120 

2.42 

320 

2.62 

330 

2.43 

490 

2.63 

500 

2.44 

630 

2.64 

640 
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Table  5b:  Two  Way  Propagation  Path  Data  for  Sphere  One  Continued, 
Collected  30  November  1987;  Run  *2.  Data  Re-measured 
Because  First  Run  Vr  (with  3.5  [V]  Supply)  Was  Clipped. 


t 

Vr 

t 

Vr 

33.68  ±0.0 1  [mS] 

9.54  ±  0.03  [V] 

33.88 

9.60 

33.69 

9.40 

33.89 

9.38 

33.70 

8.32 

33.90 

8.22 

33.71 

6.42 

33.91 

6.26 

33.72 

3.88 

33.92 

3.67 

33.73 

0.960 

33.93 

0.730 

33.74 

-2.07 

33.94 

-2.30 

33.75 

-4.90 

33.95 

-5.1 1 

33.76 

-7.24 

33.96 

-7.40 

33.77 

-8.87 

33.97 

-8.96 

33.78 

-9.64 

33.98 

-9.64 

33.79 

-9.45 

33.99 

-9.39 

33.80 

-8.32 

34.00 

-8.20 

33.81 

-6.41 

34.01 

-6.21 

33.82 

-3.83 

34.02 

-3.61 

33.83 

-0.890 

34.03 

-0.650 

33.84 

2.14 

34.04 

2.36 

33.85 

4.94 

34.05 

5.12 

33.86 

7.26 

34.06 

7.38 

33.87 

8.88 

34.07 

8.91 
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Table  6a.  Two  Way  Propagation  Path  Data  for  Sphere  Two, 
Collected  28  November  1987. 


Diameter  =  2a  =  1 0. 1 3  ±  0.0 1  [cm) 

Tc  =  21.0  ±0.1  1°  c] 

Rh  =  49.8  ±0.1  [%] 

P=  1009.4  +  0.2  [mb] 

No.  pulses  =  20 

Vsupply  =  3-51V] 

tr  =  30.95  [mS]  -  400  [pS]  =  30.55  +  0. 1  [mS] 


t 

Vt 

t 

Vt 

2.26  +  0.01  [mS] 

830  +  5  [mV] 

2.46 

830 

2.27 

770 

247 

770 

2.28 

640 

2.48 

630 

2.29 

490 

2.49 

470 

2.30 

290 

2.50 

290 

2.31 

80 

2.51 

60 

2.32 

-180 

2.52 

-190 

2.33 

-400 

2.53 

-420 

2.34 

-610 

2.54 

-600 

2.35 

-730 

2.55 

-720 

2.36 

-770 

2.56 

-770 

2.37 

-720 

2.57 

-710 

2.38 

-580 

2.58 

-570 

2.39 

-430 

2.59 

-420 

2.40 

-240 

2.60 

-230 

2.41 

-  20 

2.61 

-  10 

2.42 

240 

2.62 

240 

2.43 

460 

2.63 

480 

2.44 

660 

2.64 

670 

2.45 

780 

2.65 

790 
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Table  6b.  Two  Way  Propagation  Path  Data  for  Sphere  Two  Continued, 
Collected  28  November  1987. 


t 

vr 

t 

Vr 

34.13  iO.OUmS] 

472  ±  0.01  [V] 

3433 

474 

34.14 

466 

3434 

4.66 

34.15 

4.16 

34.35 

4.13 

34.16 

3.23 

3436 

3.19 

34.17 

2.0! 

3437 

1.95 

34.18 

0.560 

34.38 

0.500 

34.19 

-0.940 

3439 

-1.00 

34.20 

-2.35 

3440 

-2.40 

34.21 

-3.52 

3441 

-3.57 

34.22 

-436 

3442 

-4.38 

34.23 

-4.77 

3443 

-4.78 

34.24 

-4.70 

34.44 

-4.68 

3425 

-418 

34.45 

-4.15 

34.26 

-3.25 

3446 

-3.20 

34.27 

-1.99 

3447 

-1.92 

3428 

-0.560 

34.48 

-0.480 

3429 

0.960 

3449 

1.02 

34.30 

2.34 

34.50 

2.39 

34.31 

3.51 

34.51 

3.54 

3432 

434 

3452 

4.34 

M4 


Table  7:  Two  Way  Propagation  Path  Data  for  Sphere  Three, 
Collected  28  November  1987. 


Diameter  *  2a  =  7.634 ±  0.003  (cm) 

Tc  *  21 .3  i  0.1  1°  c] 

Rh  =  49.8  ±  0. 1  IX] 

P=  1009.6  ±  0.2  [mb] 

No.  pulses  =  20 

Vsupply  *  3  5  tV] 

tr  -  30.99  [mS]  -  400  [yS]  -  30.59  *  0. 1  tmSl 


Vt  data  is  the  same  as  sphere  two's  (see  Table  6a). 


t 

Vr 

t 

Vr 

33.98  ±0.0 1[mS] 

3.61  ±0.03  [V] 

34.18 

3.63 

33.99 

3.40 

34.19 

3.39 

34.00 

2.86 

34.20 

2.84 

34.01 

2.04 

34.21 

2.02 

34.02 

1.01 

34.22 

0.970 

34.03 

-0.120 

34.23 

-0.180 

34.04 

-1.25 

34.24 

-1.30 

34.05 

-2.25 

3425 

-2.31 

34.06 

-3.04 

34.26 

-3.08 

34.07 

-3.54 

34.27 

-3.56 

34.08 

-3.68 

34.28 

-3.69 

3409 

-3.46 

3429 

-3.46 

34.10 

-2.91 

34.30 

-2.89 

34.1 1 

-2.07 

34.31 

-2.04 

34.12 

-1.04 

34.32 

-0.990 

34.13 

0.1 10 

3433 

0.160 

34.14 

1.25 

34.34 

1.29 

34.15 

2.24 

34.35 

2.26 

34.16 

3.00 

34.36 

3.05 

3417 

3.48 

3437 

3.50 
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Table  8:  Two  Way  Propagation  Path  Data  for  Sphere  Four, 
Collected  28  November  1987. 


Diameter  -  2a  *  6.240  ±  0.008  {cm] 

Tc  =  21.4  ±0.1  {°  c] 

Rh  =  49.6  ±  0.1  [%] 

P=  1009.9  ±  0.2  [mb] 

No.  pulses  =  20 

^supply  *  3  5  IV1 

tr  -  3 1 .04  [mS]  -  400  [pS]  =  30.64  ±  0. 1  [mS] 

Vt  data  is  the  same  as  sphere  two's  (see  Table  6a). 


t 

Vr 

t 

vr 

33.83  ±0.0![mS] 

2.28  ±  0.03  [V] 

34.03 

2.29 

33.84 

2.10 

34.04 

2.10 

33.85 

1.68 

34.05 

1.68 

33.86 

1.10 

34.06 

1.09 

33.87 

0.420 

34.07 

0.400 

33.88 

-0.320 

34.08 

-0.340 

33.89 

-1.03 

34.09 

-1.04 

33.90 

-1.64 

34.10 

-1.66 

33.91 

-2.07 

34.11 

-2.10 

33.92 

-2.31 

34.12 

-2.32 

33.93 

-2.34 

34.13 

-2.35 

33.94 

-7.12 

34.14 

-2.12 

33.95 

-1.71 

34.15 

-1.70 

33.96 

-1.12 

34.16 

-1.11 

33.97 

-0.440 

34.17 

-0.400 

33.98 

0.320 

34.18 

0.340 

33.99 

1.02 

34.19 

1.04 

34.00 

1.62 

34.20 

1.64 

34.01 

2.06 

34.21 

2.08 

34.02 

2.28 

34.22 

2.28 
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Table  9:  Noise  Data  for  3.0  [V]  Supply  Voltage; 
Collected  30  November  1987. 


Tc  -  20.8  ±  0. 11°  c]  Rh  -  51.7  ±  0.1  [%]  P  *  1016.7  ±  0.2  [mb] 

V»ippIy  =  3  0 IV1  No.  pulses  =  20 

V,  data  is  the  same  as  sphere  one  s  (see  Table  5a). 


t 

t 

Vn 

33.73  lO.OllmS] 

00  t  5  (mV] 

34.03 

-120 

33.74 

00 

34.04 

-120 

33.75 

-  20 

34.05 

-1 10 

33.76 

-  30 

34.06 

-  80 

33.77 

-  40 

34.07 

-  70 

33.78 

-  60 

34.08 

-  40 

33.79 

-  80 

34.09 

-  30 

33.80 

-100 

34.10 

-  20 

33.81 

-1  10 

34.1 1 

00 

33.82 

-120 

34.12 

00 

33.83 

-120 

34.13 

10 

33.84 

-110 

34.14 

00 

33.85 

-100 

34.15 

00 

33.86 

-  80 

34.16 

-  30 

33.87 

-  70 

34.17 

-  40 

33.88 

-  50 

34.18 

-  70 

33.89 

-  40 

34.19 

-  90 

33.90 

-  10 

34.20 

-1 10 

33.91 

00 

34.21 

-120 

33.92 

00 

34.22 

-120 

33.93 

00 

34.23 

-120 

33.94 

00 

34.24 

-120 

33.95 

-  20 

34.25 

-1 10 

33.96 

-  30 

34.26 

-  80 

33.97 

-  50 

34.27 

-  70 

33.98 

-  70 

34.28 

-  40 

33.99 

-  80 

34.29 

-  30 

34.00 

-100 

34.30 

00 

34.01 

-110 

34.31 

00 

34.02 

-120 

34.32 

10 

1 47 


Table  10a:  Noise  Data  for  3.5  [V]  Supply  Voltage; 
Collected  30  November  1987. 


Tc  =  21.3  ±0.1  [°c] 

Rh  s  5 1 .3  ±  0.1  [%] 

P=  1016.6  ±  0.2  [mb] 
No.  pulses  =  20 


t 

Vt 

t 

Vt 

2.24  ±0.0  HmS] 

740  ±  5  [mV] 

2.44 

750 

2.25 

830 

2.45 

830 

2.26 

820 

2.46 

810 

2.27 

720 

2.47 

700 

2.28 

530 

2.48 

500 

2.29 

380 

2.49 

360 

2.30 

180 

2.50 

160 

2.31 

-  60 

2.51 

-  80 

232 

-310 

2.52 

-330 

2.33 

-530 

2.53 

-540 

2.34 

-680 

2.54 

-690 

2.35 

-770 

2.55 

-760 

2.36 

-750 

2.56 

-740 

2.37 

-640 

2.57 

-630 

2.38 

-460 

2.58 

-460 

2.39 

-300 

2.59 

-280 

2.40 

-100 

2.60 

-  80 

2.41 

140 

2.6! 

160 

2.42 

400 

2.62 

410 

2.43 

600 

2.63 

620 
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Table  10b:  Noise  Data  for  3.5  [V]  Supply  Voltage  Continued, 
Collected  30  November  1987. 


t 

Vn 

t 

V 

n 

33.80  ±0.0  UmS] 

-60  +  5  [mV] 

34.00 

-  70 

33.81 

-  80 

34.01 

-  80 

33.82 

-  80 

34.02 

-  90 

33.83 

-  80 

34.03 

-  80 

33.84 

-  70 

34.04 

-  80 

33.85 

-  60 

34.05 

-  60 

33.86 

-  40 

34.06 

-  40 

33.87 

-  20 

34.07 

-  20 

33.88 

00 

34.08 

10 

33.89 

20 

34.09 

30 

33.90 

50 

34.10 

40 

33.91 

60 

34.1 1 

70 

33.92 

60 

34.12 

60 

33.93 

60 

34.13 

70 

33.94 

40 

34.14 

50 

33.95 

30 

34.15 

40 

33.96 

20 

34.16 

10 

33.97 

00 

34.17 

-  10 

33.98 

-  30 

34.18 

-  40 

33.99 

-  50 

34.19 

-  60 
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2.  ONE  WAY  PROPAGATION  PATH  DATA 

One  way  propagation  path  data  are  presented  in  the  following  tables. 
Tc  and  Rh  are  estimated.  The  estimates  were  based  on  hygrothermographic 

comparisons  and  28  and  30  November  data. 


Table  1 1:  One  Way  Propagation  Path  Data  for  3.0  IV], 
5000  [Hz]  Continuous  Wave  Supply  from  the  HP33 1 4A 
Function  Generator;  Collected  7  December  1987. 


(Vrms)t  =  0.246  i  0.001  [V] 

T.  *  2 1 .0  [°  c] 

Rh  ~  50.0  [%] 

P  =  1021.5  ±  0.2  [mb] 

R  =  5.30  [m]  (calculated  from  two  way  propagation  path  data  for 
spheres  two,  three  and  four) 

I  (referenced  to  10"12  [W/m2]) 

r 

94.7  i  0.1  [dB] 

94.5 

94.7 
95.0 

94.8 

95.1 

94.2 
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Table  12:  One  Way  Propagation  Path  Data  for  3.5  [V], 

5000  [Hz]  Continuous  Wave  Supply  from  the  HP33 1 4A 
Function  Generator;  Collected  7  December  1987. 

(Vrms)t  =  0523±00°l  IV] 

Tc  ~  21.0  [°  cj 
Rh  *  50.0  [%} 

P=  1021.5  i  0.2  [mb] 

R  =  5.30  [m]  (calculated  from  two  way  propagation  path  data  for 
spheres  two,  three  and  four) 

lp  (referenced  to  10"12  [W/m2]) 

97.4  +  0. 1  [dB] 

96.4 

95.1 
97.0 

96.2 
97.1 

96.4 

95.4 

96.3 

97.4 
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